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Abstract 

Factors  influencing  Army  accessions  were  studied  using 
the  Army  Research  Institute's  (ARI)  Military  Entrance  Proces¬ 
sing  Station  (MEPS)  data  base,  a  US  Army  Recruiting  Command 

(USAREC)  data  base,  and  a  Defense  Manpower  Data  Center  ( DMDC ) 

\ 

data  base.  "  Six  types  of  analysis  were  performed:  regres¬ 
sion  analysis,  data  base  profile  analysis  through  SPSS  Fre¬ 
quencies  and  Crosstabs,  factor  analysis,  discriminant  anal¬ 
ysis,  time  series  analysis,  and  analysis  of  leading  economic 
indicators.  The  USAREC  Enlistment  Prediction  Model  was  im¬ 
proved:  the  MEPS  file  was  reorganized,  analyzed  with  respect 
to  District  Recruiting  Command  (DRC),  and  merged  with  the 
USAREC  file:  factor  and  discriminant  analysis  procedures 
were  developed  and  employed  to  examine  the  merged  file  by 
DRC:  Box-Jenkins  time  series  analysis  was  applied  to  the 
DMDC  data;  and  procedures  were  implemented  using  leading 
economic  indicators.  Computer  programs  were  written  by  the 
authors  to  accomplish ^the^foregoing  analyses.  Procedures 
developed  in  the  study  should  aid  USAREC  and  ARI  researchers 
and  managers  to  better  understand  the  variables  influencing 
recruiting,  and  enhance  their  ability  to  effectively  allo¬ 
cate  scarce  resources  by  Region  and  DRC.  Ultimately,  the 
information  derived  from  these  analyses  could  aid  ARI  and 
USAREC  planners  in  their  formulation  of  recruiting  policy. 


FACTORS  INFLUENCING  ARMY  ACCESSIONS 


/  I .  Introduction 

The  subject  of  this  thesis  is  accessions  and  the  factors 

influencing  a  potential  recruit's  decision  to  join  the  Army. 

The  topic  of  accessions  has  been  a  great  concern  for  Army 

planners  since  the  birth  of  the  All-Volunteer  Force.  Cur- 

» 

rently  the  United  States  Army  Recruiting  Command  (USAREC)  is 
experiencing  a  "boom"  in  recruitments,  which  can  be  attrib¬ 
uted  to  the  zealous  efforts  of  USAREC  recruiters  and  to  the 
current  economic  conditions.  However,  once  the  economic 
situation  reverses,  there  will  be  more  pressure  put  on*  the 
internal  system  to  meet  enlistment  quotas.  Whether  operating 
in  "boom"  or  "bust"  years.  Army  planners  must  be  able  to  ef¬ 
fectively  allocate  both  dollars  and  manpower  within  USAREC, 

and  the  Army  as  a  whole,  to  make  the  most  productive  use  of 

\ 

these  limited  resources. 

This  thesis  is  organized  in  the  following  manner: 

Chapter  II  Presents  general  background  information 

on  the  topic  of  accessions.  Additionally, 
it  discusses  current  related  research  that 
influenced  this  study,  and  highlights  the 
importance  of  this  effort. 

Chapter  III  Discusses  the  methodology  employed  in 
Chapters  IV,  V,  VI,  and  VII. 

Chapter  IV  Presents  regression  analysis  performed  on 
the  USAREC  data  base. 


p 


p 


I  ft 


Chapter  V 


Chapter  VI 


Chapter  VII 


Chapter  VIII 


Discusses  the  Military  Entrance  Processing 
Station  (MEPS)  applicant  file  and  the  key 
problems  encountered  in  analyzing  this 
MEPS  file.  Also  included  is  a  discussion 
of  the  merging  of  key  elements  in  the 
USAREC  data  base  with  the  MEPS  file. 
Presents  discriminant  and  factor  analysis 
performed  on  the  merged  file.  The  factor 
analysis  section  demonstrates  a  technique 
for  analyzing  the  critical  factors  affect¬ 
ing  recruiting  at  the  regional  and  the 
DRC  level.  Next,  procedures  are  developed 
for  performing  discriminant  analysis  on 
the  merged  MEPS  and  USAREC  data  to  deter¬ 
mine  the  key  variables  that  influence 
personnel  to  join  the  Army  over  the  other 
services. 

Presents  time  series  analysis  performed 
on  nationwide  accession  rates.  It  pre¬ 
sents  the  development  of  three  leading 
indicator  models  for  accession  rates. 
Highlights  the  major  accomplishments  of 
this  thesis  and  discusses  the  results  and 
conclusions  from  the  analytic  chapters. 
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II.  BACKGROUND 


General  Topic 

The  debate  over  the  All-Volunteer  Force  (AVF )  encouraged 
the  creation  of  a  presidential  commission,  named  the  Gates 
Commission,  which  was  tasked  with  developing  a  comprehensive 
plan  for  eliminating  the  draft  and  moving  toward  an  all  vol¬ 
unteer  armed  force.  In  1970  the  Gates  Commission  delivered 
its  final  report  to  the  President  and  to  Congress  with  the 
recommendation  that  first-term  military  pay  be  raised  to  a 
level  that  could  compete  with  current  levels  of  pay  received 
by  18  to  21  year  old  male  high  school  graduates.  This  would 
allow  the  Armed  Forces  to  attract  enough  qualified  volunteers 
without  the  need  for  the  draft  (Ref  4). 

The  movement  away  from  the  draft  toward  the  AVF  was  par¬ 
tially  due  to  the  increasing  resistance  to  the  Vietnam  War, 
but  was  also  due  to  a  natural  evolution  that  was  dictated  by 
the  economic  forces  of  the  human  resource  market.  Increasing 
numbers  of  men  were  reaching  draft  eligible  age  each  year, 
while  the  "winding  down"  of  the  war  moved  the  Armed  Services 
toward  decreasing  their  force  size.  These  coincident  fac¬ 
tors  operating  together  meant  that  a  smaller  percentage  of 
the  eligible  population  would  actually  serve  in  the  Armed 
Forces.  Based  on  these  conditions  and  the  Gates  Commission 
recommendations,  Congress  raised  first-term  pay  in  1971, 
which  coincided  with  plans  for  termination  of  the  draft  by 
1973.  Secretary  of  Defense  Melvin  Laird  in  fact  ended  draft 
calls  in  January  1973,  six  months  prior  to  the  expiration  of 


the  draft  authority.  The  nation  thus  became  committed  to 
the  AVF,  and  all  new  recruits  were  strictly  volunteers. 

Initial  reports  to  Congress  (Ref  5:5928)  indicated  that 
the  services  had  been  able  to  successfully  recruit  a  social- 
ly  representative  mix  of  qualified  personnel  to  meet  defense 
requirements.  Despite  the  rosy  picture  painted  to  Congress, 
the  Army  experienced  recruiting  shortfalls:  failing  to  meet 
requirements  for  gross  number  of  enlistments,  failing  to 
attract  and  enlist  sufficient  numbers  of  recruits  into  the 
combat  arms  specialties,  and  failing  to  recruit  enough  "high¬ 
er  category”  personnel  to  fill  the  Army's  increasing  number 
of  technical  positions  (Ref  4). 

Manning  the  Force  became  one  of  the  key  issues  faced  by 
the  Chief  of  Staff  of  the  Army.  This  was  evidenced  by  the 
intensification  and  expansion  of  resources,  both  in  dollars 
and  in  manpower,  allocated  to  the  U.S.  Army  Recruiting  Com¬ 
mand  (USAREC).  Greater  and  greater  pressures  were  placed 
on  recruiters  to  meet  their  mission  requirements  in  the  face 
of  continued  shortfalls,  so  much  so  that  grievous  violations 
of  recruitment  standards  occurred  leading  to  the  conviction 
and  imprisonment  of  several  recruiting  personnel  including 
officers.  Today,  only  officers  who  are  in  the  top  third  of 
their  year  group  are  selected  for  assignment  to  USAREC  and 
the  enlisted  recruiters  are  only  chosen  after  a  very  exten¬ 
sive  screening  process.  Despite  the  increased  emphasis, 
advertising,  dollars,  and  manpower,  USAREC  continued  to  ex¬ 
perience  shortfalls  until  the  economic  downturn.  An 


important,  intuitively  obvious  factor  that  influenced  enlist¬ 


ment  rates  is  the  condition  of  the  marketplace* 

Current  Related  Research 

The  topic  of  accessions  is  of  critical  importance  to 
Army  planners.  Thus  there  are  a  myriad  of  studies  in  this 
field*  An  entire  thesis  could  be  directed  to  a  literature 
review.  However,  since  this  thesis  is  an  application  effort, 
the  main  concern  is  that  no  current  research  is  being  dupli¬ 
cated.  With  this  in  mind,  the  Deputy  Chief  of  Staff  for 
Personnel's  analysis  shop.  The  Army  Research  Institute  (ARI) 
and  the  United  States  Army  Recruiting  Command  (USAREC)  were 
contacted  with  this  proposed  study  topic  (accessions).  There 
are  many  ongoing  studies  of  which  a  few  are  discussed  below 
and  both  ARI  and  USAREC  enlisted  our  services  in  studying 
accessions. 

Several  models  have  been  developed  which  show  the  rela¬ 
tionships  between  recruiting  variables  and  mission  accomp¬ 
lishment  (actual  accession  of  personnel  into  military  ser¬ 
vice).  USAREC  uses  an  automated  predictive  model  along  with 
professional  judgement  and  institutional  memory  in  the  allo¬ 
cation  of  resources  and  mission  requirements  among  its  five 
Region  Recruiting  Commands  (RRC's)  and  its  57  District  Re¬ 
cruiting  Commands  (DRC's).  This  model,  the  Enlistment 
Prediction  Model  (EPM),  and  regression  analysis  are  used  to 
produce  the  various  missions  that  are  assigned  to  the  RRC's. 
Next,  market  recruitment  potentials  are  converted  to  per¬ 
centages  of  the  command  and  specific  regional  missions  are 


assigned  to  the  DRC's.  Program  Analysis  and  Evaluation 
(PA&E)  Directorate  personnel  use  professional  judgement  to 
refine  the  mission  which  is  further  adjusted  after  negotia¬ 
tion  between  the  USAREC  Commander  and  the  DRC  Commanders. 
Current  research  on  the  relationships  between  the  various 
variables  and  the  recruiting  mission  has  not  yet  provided 
adequate  models  for  the  allocation  of  recruiters,  advertis¬ 
ing  dollars,  recruiter  aides,  etc.,  in  determining  where 
these  resources  can  be  most  productive. 

Two  recruiting  resource  allocation  models  have  been 
developed  by  USAREC:  the  EPM,  discussed  above,  and  the 
Veteran's  Educational  Assistance  Program  (VEAP)  Incentives 
Model.  The  VEAP  Model  uses  unemployment  rates  and  number  of 
active  recruiters  in  the  system  as  two  of  its  four  explana¬ 
tory  variables.  There  is  only  limited  data  on  VEAP  contracts 
(only  that  available  since  the  end  of  the  G.I.  Bill),  so  that 
simultaneous  equations  are  used  to  bade  predictions  on  term 
of  enlistment  and  the  number  of  accessions.  Once  more  data 
becomes  available,  a  regression  model  will  be  developed  with 
the  basic  VEAP  Model. 

Two  recruiting  resource  allocation  models  have  been 
developed  by  contract:  the  Claritas  Geodemographic  Model, 
and  the  N.W.  Ayer  Advertising  Effectiveness  Model  (Ayer 
Model ) . 

The  Claritas  model  uses  over  500  demographic  variables 
that  are  combined  into  34  factors.  These  factors,  combined 
with  measures  of  recruiter  strength,  proportion  of 


population  serving  in  the  Armed  Forces ,  are  used  in  a  linear 
regression  model  to  predict  the  annual  recruit  penetration 
into  the  population  of  available  males  by  mental  category. 
The  Claritas  model  was  done  under  contract  for  the  Army,  so 
that  the  specifics  of  the  model  are  known  only  to  the  con¬ 
tractor.  Claritas  uses  r  squared  as  a  measure  of  effective¬ 
ness  in  order  to  determine  where  the  Army  should  concentrate 
its  recruiting  efforts  in  order  to  maximize  accessions. 

Another  study  sponsored  by  USAREC,  was  the  evaluation 
of  the  effectiveness  of  the  Army's  National  Advertising  Ex¬ 
penditures  (the  Ayer  Model).  The  object  of  the  modelling 
project  was  to  measure  the  impact  of  dollars  spent  for  re¬ 
cruiting  on  the  quantity  and  quality  of  recruits.  A  second¬ 
ary  objective  was  to  measure  and  compare  the  influences  of 
other  major  factors  such  as  pay,  size  of  the  recruiting 
force,  type  and  extent  of  recruiting  mission,  and  the  amount 
of  money  spent  on  advertising  in  order  to  determine  the  most 
cost  effective  way  of  increasing  the  quality  of  accessions. 
The  model  assumed  that  any  omitted  variables  were  not  domi¬ 
nating  influences  and  that  past  (1976-1980)  relationships 
will  persist,  staying  relatively  close  to  their  observed 
ranges.  The  main  findings  of  this  study  were  that  dollars 
spent  on  advertising  were  the  most  cost  effective  way  of 
increasing  the  quantity  and  quality  of  accessions.  It 
should  be  noted  that  N.W.  Ayer  International  is  the  adver¬ 
tising  agency  that  handles  a  majority  of  the  Army  adver- 


The  Optimal  Recruiter  Allocation  Model  (ORAM)  was 
developed  to  show  the  relationship  between  the  many  depen¬ 
dent  variables  and  the  District  Recruiting  Command's  mis¬ 
sion.  (Ref  23).  Its  findings  suggest  that  the  optimal  allo¬ 
cation  of  recruiters  would  move  recruiters  from  the  Northern 
and  Western  Regions  and  place  them  into  the  two  southern  re¬ 
gions.  This  seems  to  tie  in  with  the  1980  census  data  which 

* 

indicated  population  shifts  into  the  HSun  Belt",  and  with 
the  traditional  conservatism  associated  with  the  South  and 
its  support  of  the  military.  ORAM  uses  past  performance  to 
allocate  recruiters  to  areas  where  recruiting  has  been  suc¬ 
cessful. 

Two  models  developed  for  the  Navy  are  Morey's  Budget 
Allocation  and  Enlistment  Prediction  Model,  and  Goldberg's 
Navy  Enlisted  Supply  Study.  All  of  the  models  discussed 
above  are  summarized  in  Ref  16. 

\ 

Dale  and  Gilroy  (Ref  5)  present  an  application  of  time 
series  analysis  used  to  predict  accession  rates.  This  was 
a  working  paper  which  the  authors  permitted  us  to  use  as  a 
source  in  this  thesis.  An  objective  of  the  paper  was  to 
study  the  relationship  between  economic  variables  and  acces¬ 
sion  rates.  Unemployment  rates  were  one  of  the  key  indepen¬ 
dent  variables  in  this  study  and  the  dependent  variable  was 
male  high  school  graduate  mental  category  I-IIIA  accessions. 
Regression  analysis  was  used  lagging  unemployment  rates  to 
determine  the  correlation  between  accessions  and  the  lagged 
unemployment  rates. 


The  Dale  and  Gilroy  study  was  important  to  this  thesis 
effort  because  it  is  a  current  working  project  and  it  stud¬ 
ied  one  of  the  variables  important  to  this  effort,  unemploy¬ 
ment  rates. 

Problem  Statement  and  Research  Objectives 

The  thrust  of  this  thesis  is  identified  as  follows: 

1.  Study  accessions  in  an  attempt  to  understand  what 
factors  influence  them. 

2.  Build  a  discriminant  function  to  predict  whether  a 
potential  accession  will  join  the  Army  or  other  services. 

3.  Enhance  the  USAREC  Enlistment  Prediction  Model. 

4.  Conduct  time  series  analysis  on  accessions  using 
Box  and  Jenkins  techniques. 

5.  Developed  leading  indicator  models  to  predict  ac¬ 
cessions  based  on  the  following  input  variables:  unemploy¬ 
ment  rates,  the  producer  price  index  and  the  prime  rate. 

6.  Provide  to  ARI  and  USAREC  usable  computer  programs 


to  accomplish  the  above 


III.  METHODOLOGY 


Introduction 

This  chapter  presents  the  methodology  employed  in  the 
following  four  chapters. 

The  Data 

Initially  this  study  started  as  an  analysis  of  a  ten 
percent  random  sample  of  the  Military  Entrance  Processing 
Station  (MEPS)  Applicants  file  which  was  received  from  the 
Army  Research  Institute  (ARI).  This  file  consists  of  up  to 
256  variables  of  information  about  individuals  who  show  an 
interest  in  the  armed  services.  This  interest  is  manifested 
by  the  applicant  arriving  at  a  MEPS  to  begin  enlistment  pro¬ 
cessing.  Approximately  58.18  percent  of  the  applicants  turn 
into  accessions.  There  are  48,520  applicants  or  cases  in 
the  file,  so  that,  in  order  to  facilitate  some  form  of  a 
cost  effective  analysis,  the  number  of  variables  per  case 
was  reduced  to  what  was  considered  a  manageable  and  minimum 
number  of  essential  variables.  These  included  such  variables 
as  date  of  birth,  sex,  service  processed  for,  etc. 

Additionally,  the  United  States  Army  Recruiting  Command 
(USAREC)  was  visited  and  another  data  base  was  obtained, 
which  will  be  referred  to  as  the  USAREC  file.  The  USAREC 
file  is  input  into  an  SPSS  regression  routine  to  predict  var¬ 
ious  levels  of  accessions.  This  model  is  the  initial  phase 
in  the  quarterly  process  of  assigning  the  recruiting  mission 
to  the  fifty-seven  District  Recruiting  Commands  (DRC's). 


The  focal  point  of  this  phase  of  the  analysis  was  to  enhance 
the  USAREC  regression  model.  This  is  discussed  in  Chapter  IV. 

There  were  many  variables  in  the  USAREC  file  that  could 
be  applied  to  individuals  on  the  MEPS  Applicant  File.  For 
instance,  the  USAREC  file  lists  the  unemployment  rate  for 
the  specific  DRC's  for  specific  quarters.  This  is  an  impor¬ 
tant  variable  because  it  reflects  the  unemployment  rate  that 
was  influencing  the  applicant  to  make  contact  with  a  recruit¬ 
er.  The  problem  became  to  merge  desired  variables  in  the 
two  files  into  a  new  file.  This  process  is  discussed  in  de¬ 
tail  in  Chapter  V.  The  major  stumbling  block  in  this  step 
was  the  fact  that  the  variables  in  the  USAREC  file  relate 
to  DRC's,  whereas  the  variables  in  the  ARI  file  relate  to 
individuals,  and  the  individual  MEPS  records  are  not  catego¬ 
rized  by  DRC.  The  task  became  to  equate  an  individual  MEPS 
record  to  a  DRC.  Once  this  task  was  accomplished,  meaning¬ 
ful  analysis  was  performed  on  the  merged  file  using  discrim¬ 
inant  and  factor  analysis  routines  from  SPSS.  These  results 
are  discussed  in  Chapter  VI. 

Another  objective  of  this  study  was  to  perform  time 
series  analysis  on  accession  data.  The  USAREC  file  contained 
only  22  quarterly  data  points  which  is  not  an  adequate 
amount  to  perform  time  series  analysis.  A  third  data  set  of 
monthly  accession  figures  from  October  1975  to  July  1982  was 
obtained  from  the  Defense  Manpower  Data  Center  ( DMDC ) ,  US 
Army  MILPERCEN  Liaison  Officer,  Monterey,  California.  The 
81  data  points  were  utilized  in  the  analysis  performed  in 
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Chapter  VII.  When  tru, ielling  accessions  using  leading  indi¬ 
cator  techniques,  the  following  variables  were  used  as  input 
variable: 


Variable  Journal  Source 

National  Unemployment  Rates  Federal  Reserve  Bulletin 

Producer  Price  Index  Monthly  Labor  Journal 

Prime  Rate  Monthly  Labor  Journal 

There  was  a  slight  problem  to  be  resolved  when  using 
these  variables.  From  month  to  month  the  rate  can  change 
due  to  revised  calculations.  Therefore,  the  last  time  a 
rate  or  index  appears  in  a  journal  is  the  rate  used  in  the 
analysis.  For  example,  in  the  September  issue,  for  the 
first  time,  is  listed  the  unemployment  rate  for  August. 
However,  the  August  rate  will  appear  in  several  future  edi¬ 
tions  and  may  be  changed.  In  June  of  the  following  year  the 
August  rate  will  no  longer  be  shown.  Thus,  this  analysis 
used  the  rate  from  the  May  issue. 

A  sample  of  the  USAREC  file  is  shown  in  Appendix  A-l. 
Complete  listings  of  accessions,  unemployment  rates,  the 
producer  price  index,  and  the  prime  rate  are  shown  in 
Appendices  A-2  through  A-5. 

Analysis  Techniques 

All  analyses  were  performed  on  a  Cyber  CDC  6600  using 
SPSS  for  the  regression,  factor  and  discriminant  analysis. 
Biomedical  Computer  Programs  P-Series  ( BMDP )  (Ref  24)  was 
used  along  with  FORTRAN  time  series  computer  programs  (re¬ 
ferred  to  as  TSP )  that  utilize  various  International 


Mathematical  and  Statistical  Library  (IMSL)  subroutines 
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Regression 

It  is  felt  that  regression  analysis  is  a  common  method 
of  analysis,  therefore,  it  is  not  expounded  upon  in  this 
paper. 

File  Merge 

The  USAREC  file  and  the  MEPS  applicant  file  were  merged 
in  order  to  perform  multivariate  data  analysis.  Thi  process 
is  discussed  in  detail  in  Chapter  V.  The  merged  file  is  re¬ 
ferred  to  as  MU  for  MEPS-USAREC. 

Factor  Analysis 

Factor  analysis  was  performed  on  the  MU  file  using  the 
SPSS  subroutine  Factor  to  discover  how  the  various  variables 
in  the  analysis  "hang  together”  to  form  what  are  called  fac¬ 
tors  or  principal  components.  The  SPSS  method  of  factoring 
chosen  was  "Principal  Component  Analysis",  i.e.,  TYPE=PA1. 
Factor  analysis  can  be  used  to  reduce  the  dimensionality  of 
the  data  sample  by  condensing  out  principal  components  from 
the  patterns  and  relationships  found  among  the  observed 
variables  in  the  sampled  data.  This  type  of  analysis  is 
similar  to  using  "cluster  analysis"  on  a  correlation  matrix 
to  pick  out  clusters  of  variables  that  have  high  correlation 
coefficients  among  each  other,  and  calling  the  resulting 
cluster  a  "factor."  The  MU  file  was  examined  using  factor 
analysis  in  order  to  add  to  the  overall  MU  data  profile  and 
to  enhance  the  cross-tabulation  and  discriminant  analysis. 


l3 


Thus,  this  study  supports  the  thesis  objective  of  giving 
ARI  and  USAREC  a  detailed  overview  of  the  important  factors 
relative  to  recruitment.  Chapter  VI  contains  the  results 
of  the  factor  analysis  performed  on  the  MU  file. 


Discriminant  Analysis 

Discriminant  analysis  was  performed  on  the  MU  file 
using  the  SPSS  subroutine  Discriminant  to  examine  how  well 
the  different  MU  variables  could  predict  that  an  individual 
would  join  the  Army  over  another  branch  of  the  service,  and 
to  discover  the  differences  among  these  predictive  variables 
by  region.  The  fundamental  aim  of  discriminant  analysis  is 
to  formulate  a  rule  that  can  be  used  to  correctly  classify 
individual  cases  into  two  or  more  groups  based  on  attribute 
variables  associated  with  the  given  cases.  In  SPSS  the  rule 
that  is  formulated  consists  of  assigning  both  standardized 
and  unstandardized  weights  to  the  attribute  variables  that 
can  then  be  used  to  predict  into  which  group  a  randomly 
sampled  case  may  fall.  Under  SPSS  there  are  two  ways  to 
classify  groups,  given  raw  input  data:  either  by  using 
Fisher's  Linear  Discriminant  Coefficients,  or  by  using  the 
Canonical  Discriminant  Function  Coefficients.  If  Fisher's 
Coefficients  are  used,  then  the  highest  of  the  computed 
functional  values  will  show  what  group  the  input  case  is  to 
be  classified.  If  the  Canonical  Coefficients  are  used,  then 
for  the  case  of  two  groups,  one  can  classify  the  input  case 
based  on  the  sign  of  the  resulting  functional  value  since 
the  overall  population  mean  (population  centroid)  will  have 
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a  value  of  zero.  In  this  study  the  Canonical  Coefficients 
rather  than  Fisher's  Coefficients  are  tabled  in  Appendices 
C-9  through  C-12.  When  the  unstandardize  Canonical  Coeffi¬ 
cients  are  used,  only  one  equation  and  set  of  coefficients 
is  required  to  compute  which  group  a  random  case  belongs  to. 
The  result  of  such  a  computation  on  raw  input  data  will  be 
either  positive  showing  it  belongs  to  one  group  or  negative 
showing  it  belongs  to  the  other  group.  A  zero  result  would 
point  to  neither  group.  The  standardized  weights  show  the 
intensity  with  which  a  particular  variable  is  able  to  help 
predict  into  what  group  a  case  belongs.  The  unstandardized 
weights  are  coefficients  which  are  very  much  like  the  coef¬ 
ficients  in  regression  analysis  that  are  used  to  compute  a 
dependent  variable;  but  in  discriminant  analysis  the  comput¬ 
ed  variable  is  compared  to  the  group  means  to  show  where  it 
falls  relative  to  them,  thus  predicting  to  which  group  a 
case  "belongs."  ' 

To  develop  and  later  test  the  discriminating  equation, 
the  MU  file  is  randomly  divided  into  two  halves.  The  first 
half  of  the  file  is  used  to  model  the  discriminating  equa¬ 
tion,  then  the  first  half  is  run  against  the  formulated 
equation  to  see  what  percentage  of  the  records  are  correct¬ 
ly  classified  by  the  equation.  The  second  half  of  the  file 
is  likewise  run  against  the  same  formulated  equation  to 
determine  how  well  the  equation  correctly  classifies  a  ran¬ 
dom  data  sample  other  than  the  random  sample  that  created 
the  equation. 
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The  SPSS  procedure  includes  a  listing  of  both  the 
standardized  coefficients  as  well  as  the  values  of  the  group 
centroids,  i.e. ,  the  values  of  the  discriminant  function 
evaluated  at  the  group  means.  There  are  also  five  impor¬ 
tant  statistical  values  generated  by  the  discriminant  sub¬ 
routine: 

1.  Canonical  Correlation, 

2.  Wilks  Lambda, 

3.  Chi-Squared, 

4.  Degrees  of  Freedom,  and 

5.  the  corresponding  Significance. 

The  last  value  is  a  percentage  used  in  hypothesis  testing 
such  that  if  a  test  is  conducted  at  alpha  equal  .05  and  the 
Q:  resulting  Significance  is  greater  than  .05  then  the  null 

hypothesis  will  be  accepted.  Specifically,  in  the  case  of 
discriminant  analysis,  the  null  hypothesis  is  that  the  group 
centroids  are  equal,  i.e.,  there  is  no  significant  differ¬ 
ence  between  the  groups.  However,  if  the  Significance  is 
less  than  .05  then  the  null  hypothesis  will  be  rejected. 
Specifically,  it  is  concluded  that  the  discriminant  func¬ 
tions  provide  a  significant  differentiation  among  the  groups. 
The  Significance  value  in  this  study  is  a  resulting  measure 
of  the  following  test  of  hypothesis: 

H0  !  U1  -  U2 
Hi  :  ui  *  u2 

where  =  the  group  centroid  for  non-Army  records 
and  U2  =  the  group  centroid  for  Army  records. 
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To  test  this  hypothesis  a  chi-squared  statistic  is 
computed  by  SPSS  and  compared  to  the  tabled  value  of  chi- 
squared  at  a  given  alpha  and  with  p*(g-l)  degrees  of  freedom, 
where  p  is  the  number  of  parameters  and  g  is  the  number  of 
groups.  If  the  computed  chi -squared  is  greater  than  the 
table  value  of  chi -squared,  the  null  hypothesis  will  be  re¬ 
jected.  If  the  null  hypothesis  is  rejected,  then  it  can  be 
concluded  that  the  discriminant  functions  do  indeed  provide 
a  significant  differentiation  among  the  groups. 

The  discriminant  coefficients,  centroid  values,  and 
forementioned  statistics  are  the  primary  tools  used  in  this 
portion  of  the  study  to  examine  the  MU  file  and  to  try  to 
discover  what  variables,  if  any,  had  an  impact  on  differen¬ 
tiating  the  Army  from  the  non- Army  prospective  recruits. 

The  Mahalonobis  distance  stepwise  method  was  employed 
in  all  discriminant  analyses  performed  in  this  study.  This 
method  allowed  only  the  most  significant  variables  to  enter 
in  each  individual  analysis,  thus  permitting  very  refined 
differentiation  by  DRC.  Since  only  the  most  important  vari¬ 
ables  within  each  DRC  were  selected,  a  clearer  definition 
of  the  demographic  differences  amoncNaRC 1 s  can  also  be 


developed.  The  selection  criteria  for  the 


:y  of  a  given 


variable  in  the  stepwise  procedure  depended  on  having  a 
partial  F  ratio  of  greater  than  or  equal  to  1.0,  while  the 
tolerance  level  was  set  at  .001  (Ref  20:448-454).  Chapter  VI 
contains  the  results  of  the  discriminant  analyses  performed 
on  the  MU  file. 


The  goal  of  the  time  series  analysis  was  to  develop  a 
model  which  predicts  accessions  based  upon  previous  values 
of  accessions.  A  FORTRAN  program  called  TSP  was  written 
based  upon  Box-Jenkins  (Ref  1)  time  series  analysis  tech¬ 
niques.  This  program  uses  IMSL  subroutines  to  calculate 
autocovariances,  autocorrelations  and  partial  autocorrela¬ 
tions.  These  values  are  graphed  and  examined  for  signifi¬ 
cant  values.  Additionally,  Gray,  Kelley  and  Mclntire's  R  and 
S  array  techniques  for  model  identification  (Ref  9)  have  been 
incorporated  into  this  program  listed  in  Appendix  A-6.  Once 
the  model  has  tentatively  been  identified,  BMDP  time  series 
routines  are  used  to  further  analyze  the  model.  Once  a 
final  model  has  been  chosen,  adequacy  checks  are  conducted, 
using  FORTRAN  programs  and  BMDP.  These  checks  included  the 
Portmanteau  lack  of  fit  test  and  the  cumulative  periodogram 
plot  (Ref  1:290-298)  and  analysis  of  the  residual  series  for 
significant  autocorrelations.  This  time  series  approach  is 
applied  in  Chapter  VII. 

Leading  Indicator  Analysis.  In  modelling  accessions 
using  a  leading  indicator,  the  desired  results  are  two-fold: 
first,  determination  of  an  input  variable  that  leads  acces¬ 
sions,  in  order  to  allow  use  of  current  values  of  the  input 
variable  to  predict  future  values  of  accessions,  and  second, 
tightening  of  the  standard  error  band  about  the  predicted 
values  of  accessions.  Once  again  BMDP  leading  indicator  tech¬ 
niques  were  used.  The  steps  in  this  modelling  process  are 
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as  follows: 

1*  Model  the  input  variable  as  an  Autoregressive  Inte¬ 
grated  Moving  Average  Process  ( ARIMA )  to  obtain  the  pre¬ 
whitened  input  series,  say  R^tt). 

2.  Apply  the  ARIMA  model  to  accessions  and  denote  the 
filtered  output  series  as  Ry(t). 

3.  Cross  correlate  R  (t)  and  R  (t)  attempting  to  iden- 

x  y 

tify  the  functional  form  of  the  transfer  function  component. 

4.  Identify  an  ARIMA  model  for  Ry(t)  and  cobine  the 
ARIMA  component  with  the  transfer  function  component  to  ob¬ 
tain  a  tentative  transfer  function  model. 

5.  Apply  autocorrelation  and  cross  correlation  checks 
(Ref  1:392-298)  to  test  for  model  adequacy. 

6.  Forecast  future  values  of  accessions. 

All  of  these  steps  were  performed  using  BMDP.  This 
technique  is  applied  in  Chapter  VII,  using  accessions  as  the 
output  variable  and  unemployment  rates,  producer  price  index, 
and  the  prime  rate  as  input  variables. 

Statistical  Significance 

In  any  study,  results  could  appear  that  are  attribut¬ 
able  solely  to  chance  or  sampling  variance.  In  order  to 
minimize  this  possibility,  tests  described  in  this  report 
are  conducted  at  an  alpha  level  of  .05  in  Chapter  IV  and 
Chapter  V.  Therefore  when  something  is  referred  to  as  being 
significant  or  statistically  significant,  this  equates  to  a 
less  than  five  percent  likelihood  that  the  result  is  due  to 
chance.  Unless  otherwise  clearly  indicated,  significance 
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relates  to  an  alpha  equals  five  percent  level.  In  the  time 
series  analysis  (Chapter  VII)  when  discussing  autocorrela¬ 
tion  plots,  significance  relates  to  an  approximate  two  stan¬ 
dard  error  band. 

Assumptions 

The  following  assumptions  are  made  in  this  report: 

1.  The  relationships  that  exist  between  the  variables 
examined  over  the  time  periods  examined  are  not  changing. 

2.  The  results  from  the  BMDP  time  series  analysis  are 
valid. 

3.  The  results  from  TSP  are  valid. 

The  second  and  third  assumptions  are  of  particular 
concern.  In  order  to  validate  BMDP  and  TSP,  the  results  of 
Box -Jenkins  Gas  furnace  example  (Ref  1:381-409)  were  com¬ 
pared  to  results  of  the  same  example  from  BMDP  and  TSP  anal¬ 
ysis.  This  test  validated  the  acceptability  of  both  BMDP 
and  TSP. 

Choice  of  Dependent  Variable 

There  are  numerous  categories  of  accessions  but  the 
time  series  and  regression  segment  of  this  report  studies 
only  one:  Male  High  School  Graduate  Mental  Category  I-IIIA 
accessions  (CAT13AM ) •  This  is  done  for  two  reasons:  first, 
for  continuity  in  the  analytic  phases  of  the  study,  and. 


second,  because  the  other  categories  of  accessions  such  as 
female,  mental  category  IV,  and  prior  service  accessions 
are  demand  constrained  whereas  CAT13AM  accessions  are  supply 


constrained.  Recently  the  Army  has  been  able  to  focus  al¬ 
most  all  of  its  recruiting  effort  on  CAT13AM  accessions  and 
has  eliminated  the  accession  of  non-high  school  graduate 
and  mental  category  IV  personnel.  One  of  the  underlying 
objectives  of  this  study  is  for  it  to  be  useful  for  Army 
personnel  management  policy  makers  and  planners.  Therefore 
the  dependent  variable  studied  in  this  report  is  CAT13AN. 


REGRESSION  ANALYSIS 


One  of  the  main  objectives  of  this  report  is  to  take 
the  current  model  that  the  United  States  Army  Recruiting  Com 
mand  (USAREC)  uses  to  predict  Male  HSG  Mental  Category  I- 
IIIA  accessions  and  improve  on  its  capability  to  account  for 
the  variance  in  accessions.  The  regression  model  is  called 
the  Enlistment  Projection  Model  (EPM)  and  it  is  the  first 
step  in  a  process  that  USAREC  planners  and  managers  use  to 
assign  quarterly  recruiting  missions  to  the  various  District 
Recruiting  Commands  (DRC).  The  EPM  and  data  base  were  pro¬ 
vided  for  use  in  this  thesis  by  USAREC  Plans  and  Analysis 
Section.  The  focal  point  of  this  study  is  essentially  Male 
HSG  Mental  Category  I-IIIA  accessions  so  that  we  only  con¬ 
cern  ourselves  with  this  variable  as  a  dependent  variable. 
Table  1  is  a  summary  of  the  variables  used  in  the  model  and 
a  brief  explanation  of  their  meaning.  In  Appendix  A-l  is  a 
sample  of  the  USAREC  data  base  and  its  associated  input  for¬ 


mat. 


The  regression  design  statement  used  by  USAREC  in  their 


SPSS  Program  was  modified  so  as  to  specify  an  inclusion  and 
exclusion  level  for  variables  to  enter  and  leave  the  model 


and  is  as  follows: 

REGRESSION  METHOD=STEPWISE/ 

VARIABLES® AREA, REACT ,HSSNR ,QMA ,RCTRS ,NPSMACC, 
HSDMACC .UNEMP , RCTREX , INCOME , BLACKS , PROPEN , QlVAR , 
Q2VAR , Q3 VAR , DODRCR , AIDES , DODN , DODH , TOTCON , HSGCON , 
ADVCOST ,DRCADV,CAT13AM .PSTOTACC .NPSFACC/ 
REGRESSI0N=CAT1 3AM (*,2.0, *,1.99)  WITH  AREA.HSSNR, 
QMA , UNEMP , RCTREX , INCOME , BMA , PROPEN , QlVAR , Q2 VAR , 

Q3 VAR , DODRCR , AIDES , ADVCOST , DRCADV , ACTR5 , REACT/ 


The  model  was  run  for  the  five  recruiting  regions  with 


following  results: 


Reaion 

Adjusted  r 

NE 

.75220 

SE 

.73807 

SW 

.68286 

MW 

.62862 

WEST 

.75133 

Final  models  for  each  region  are  in  Appendix  B-l.  It  is 
important  to  note  that  we  are  not  attempting  to  analyze  the 
resulting  equations  that  model  accessions  in  the  various 
regions.  However,  analysis  of  the  results  and  the  signifi¬ 
cance  of  various  variables  in  the  models  did  aid  in  the 
following  chapters  of  this  report. 

Studying  the  results  from  the  USAREC  model  and  consider 
ing  the  meaning  of  some  of  the  variables,  we  decided  to  de¬ 
velop  a  new  variable  called  density  that  would  be  a  reflec¬ 
tion  of  some  of  the  variables  already  in  the  model.  Our 
perception  was  that  several  of  the  variables  alone  had  hard- 
to-interpret  meanings  whereas  in  some  functional  relation¬ 
ship  with  other  variables  they  would  have  more  meaning.  For 
instance,  in  several  models  DRC  area  in  square  miles  (AREA) 
is  significant  as  is  on-production  recruiters  ( RCTRS )  and 
qualified  military  available  (QMA).  It  seems  to  make  sense 


that  accessions  might  be  tied  to  some  sort  of  relationship 

RCTRS 

between  these  variables  such  as  Density  =  or  Density 

RCTRS 

AREA-"  *  Density  was  computed  as  indicated  and  added 

to  the  data  and  the  regressions  were  run  again.  Addition¬ 


ally,  consideration  was  given  to  the  possibility  that 


VARIABLE 

NAME 

TIME  QTR 

RGN , DRC 

AREA 

QMA 

REACT 

HSSNR 

RCTRS 

NPSMACC 

HSDGMACC 

UNEMP 

RCTREX 

RCTRSPCT 

DODNPS 

DODHSDG 

INCOME 

BMA 

TOTCON 

HSGCON 

PROPEN 
AIDES 
C ATI 3 AM 
PSTOTACC 
NPSFACC 

DRCADV 


DATABASE 

COLUMNS 

1-4 

5-6 

7-11 

12-16 

17-21 

22-26 

27-30 

31-36 

37-41 

42-44 

45-47 

48-50 

51-56 

57-62 

63-67 

68-72 

73-76 

77-80 

81-83 

84-87 

88-91 

92-95 

96-99 

100-103 


VARIABLE  DESCRIPTION 

FISCAL  YEAR  &  QUARTER 

DISTRICT  RECRUITING  COMMAND  CODE 

DRC  AREA  IN  SQUARE  MILES 

QUALIFIED  MILITARY  AVAILABLE 

REACT  LEADS  (MAGAZINE  ADS) 

MALE  HIGH  SCHOOL  SENIORS 

ON-PRODUCTION  ARMY  RECRUITERS 

ARMY  NON-PRIOR  SERVICE  MALE  ACCESSIONS 

ARMY  HIGH  SCHOOL  DIPLOMA  GRADUATE 
MALE  ACCESSIONS 

DRC  OVERALL  UNEMPLOYMENT 

PERCENT  OF  DRC  RECRUITERS  WITH  ONE 
OR  MORE  YEARS  EXPERIENCE 

ARMY  RECRUITERS  AS  PERCENT  OF  DOD 
RECRUITERS  BY  DRC 

DOD  NON-PRIOR  SERVICE  MALE  ACCESSIONS 

DOD  HIGH  SCHOOL  DIPLOMA  GRADUATE 
MALE  ACCESSIONS 

MEDIAN  DISPOSABLE  FAMILY  INCOME 

BLACK  MILITARY  AVAILABLE 

ARMY  TOTAL  CONTRACTS 

ARMY  HIGH  SCHOOL  DIPLOMA  GRADUATE 
MALE  CONTRACTS 

ACTIVE  ARMY  ENLISTMENT  PROPENSITY 

HOMETOWN  RECRUITER  AIDES  (AVERAGE) 

ARMY  MC  1-3A  MALE  ACCESSIONS 

ARMY  PRIOR  SERVICE  TOTAL  ACCESSIONS 

ARMY  NON-PRIOR  SERVICE  FEMALE 
ACCESSIONS 

DRC  LOCAL  ADVERTISING  EXPENDITURES 


Table  I.  Variables  in  USAREC  Model 


accessions  are  not  driven  by  a  linear  function  but  rather 
that  it  might  be  some  sort  of  a  multiplicitive  or  Cobb-Douglas 
type  production  function  of  the  form: 


where  Y  represents  accessions,  the  xn's  the  various  indepen¬ 
dent  variables,  and  the  an's  their  associated  parameters. 

The  equation  is  solved  by  ordinary  least  squares  regression 
after  utilizing  a  logarithmic  transformation  as  follows: 

LnY  =  a  +  a.lnX.  +  a-lnX-  +  ...  +  a  lnX 
01122  n  n 

The  coefficients  generated  by  SPSS  for  the  an's  represent 
each  of  the  variables  elasticity  with  relationship  to  acc¬ 
essions.  All  of  the  models  were  run  again  taking  a  log¬ 
arithmic  transformation  of  all  variables.  In  addition,  a 
model  was  estimated  with  both  densities  in  the  equation. 

The  resulting  adjusted  r-square  for  all  models  is  in  Table  II. 
The  SPSS  output  for  only  the  final  iteration  of  optimal 
models  discussed  below  is  in  Appendix  B-2. 

From  Table  H  it  is  clearly  seen  that  in  all  cases  the 
USAREC  Model  has  been  improved.  By  saying  " improved”  it  is 
meant  that  more  of  the  variance  in  the  dependent  variable 
is  being  accounted  for  by  the  independent  variables.  The 
model  that  does  the  best  overall  job  is  Model  4,  the  log 
transformation  of  the  USAREC  model  including  the  variable 

Density  =  (RCTRS*QMA)/AREA.  For  the  northeast,  sourthwest 

2 

and  west  regions  the  adjusted  r  value  is  the  highest  value 

2 

and  for  the  midwest  region  the  adjusted  r  value  is  within 


•0082  of  the  largest  value.  That  is  not  to  imply  that  the 

key  is  the  density  variable.  The  summary  of  the  final  model 

for  the  northeast  and  midwest  regions  shows  that  density  did 

not  enter  the  equation  modelling  accessions  (Model  1  and  4) 

and  thus  shows  that  the  logarithmic  transformation  is  respon- 

2 

sible  for  the  higher  adjusted  r  value.  For  the  southeast 

2 

region  model  3  optimizes  the  adjusted  r  value.  Model  3  is 
the  USAREC  model  plus  a  Denisty  =  (RCTRS*QMA)/( AREA) .  In 
this  model,  Denisty  proved  to  be  a  very  strong  variable, 
entering  the  model  at  step  6  with  a  F  value  of  44.84. 

For  the  sake  of  simplification,  models  3  and  4  will  be 
combined  into  one  SPSS  regression  program  with  two  regression 
design  statements:  one  with  a  logarithmic  transformation  of 
the  USAREC  Model  plus  Density  =  (RCTRS*QMA)/AREA  which  is  to 
be  used  to  predict  accessions  for  all  regions  except  the 
southeast  region  and  the  second  with  the  USAREC  model  plus 
Denisty  =  (RCTRS*QMA)/AREA  which  is  to  be  used  for  the  south¬ 
east  region.  A  complete  listing  of  this  program  is  at 
Appendix  B-3.  This  model  will  be  provided  to  USAREC  for  use 
as  a  new  Enlistment  Projection  Model. 

Having  evaluated  all  these  models,  there  are  some  inter¬ 
esting  relationships  between  accessions  and  the  various  in¬ 
dependent  variables.  It  might  be  expected  that  the  more 
dollars  spent  on  advertising  the  more  accessions  the  Army 
would  have,  yet  it  appears  that,  in  the  midwest,  dollars 
spent  on  advertising  has  a  negative  effect  on  accessions. 

Table  3  is  a  listing  of  variables  that  have  unexpected  negative 


effects  on  accession  rates.  Additionally,  in  the  southeast 
and  southwest  regions,  the  number  of  react  leads  has  a  neg¬ 
ative  impact  on  accessions.  A  react  lead  is  counted  when 
an  individual  responds  to  any  form  of  magazine  advertising. 
It  is  these  types  of  relationships  that  prompted  the  next 
portion  of  this  study.  We  attempt  to  determine  what  vari¬ 
ables  are  factors  affecting  accessions  and  to  classify  pos¬ 
sible  recruits  as  being  likely  to  join  the  Army  or  not. 

Results  and  Conclusions 

The  USAREC  model  was  enhanced  by  computing  a  new  vari¬ 
able  Density  *  (RCTRA*QMA)/AREA  for  the  Southeast  region  and 
for  all  other  regions  by  taking  a  logarithmic  transformation 
of  the  previous  model.  However,  there  were  several  interest 
ing  and  unexplainable  relationships  that  appear  in  the 
models . 

Recommendation  for  Future  Research 

There  are  many  other  variables  that  might  further  en¬ 
hance  the  USAREC  Model  but  due  to  time  constraints  we  were 
not  able  to  pursue  them.  Some  examples  are  a  wage  dif¬ 
ferential  variable  and  a  percent  annual  pay  raise  variable. 
From  this  analysis  questions  arise  as  to  what  factors  relate 
to  Army  accessions  and  how  do  they  compare  with  the  other 
services?  What  type  of  a  person  is  likely  to  enter  the 
Army?  These  questions  will  be  addressed  in  Chapter  VI  of 
this  report. 
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Table  III.  Variables  With  Regression  Coefficients 
With  Unexpected  Negative  Sign 
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Data  Base  Preparation 
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Topic  Selection  and  Data  Gathering 

One  of  the  many  ongoing  activities  of  the  Army  Research 
Institute  (ARI)  in  Alexandria,  Virginia,  is  its  investigation 
of  the  Military  Entrance  Processing  Stations  (MEPS)  data 
base  to  aid  in  answering  a  continuing  chorus  of  questions 
at  the  highest  levels  of  personnel  management  in  the  Army. 
The  Personnel  Research  and  Analysis  Section  at  ARI  (office: 
ARI— PERI— RP )  suggested  that  their  long  run  investigation 
could  be  aided  through  the  employment  of  multivariate  anal¬ 
ysis  techniques  on  a  sample  of  the  MEPS  file.  It  was  recom¬ 
mended  that  this  study  could  most  effectively  be  carried  out 
by  first  building  a  "profile"  of  the  general  structure  and 
content  of  the  data  base.  Several  areas  of  current  interest 
to  ARI  were  suggested  as  objectives  for  the  investigation. 
Two  of  these  objectives  formed  the  basis  of  this  thesis 
study: 

1.  Determine  what  factors  influence  fully  qualified 
individuals  to  join  the  Army  over  other  branches  of  the 
Military  Service ,  and 

2.  Improve  the  Army's  current  predictive  model  of  re¬ 
cruitment  projections. 

To  aid  in  the  pursuit  of  these  objectives,  it  was  re¬ 
commended  that  more  data  and  information  be  obtained  from 
the  operations  researchers  of  the  US  Army  Recruiting  Command 
(USAREC)  at  Fort  Sheridan,  Illinois.  USAREC  (office:  USAREC' 
PAE-RE)  provided  a  complete  copy  of  their  recruiting  model 
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input  file,  while  ARI-PERI-RP  was  able  to  provide  a  useful 
random  sample  of  the  MEPS  data  base. 

The  data  gathering  phase  also  included  an  examination 
and  modeling  of  a  typical  MEPS  operation;  and  the  MEPS  site 
chosen  for  this  research  was  Cincinnati.  There,  information 
was  obtained  on  the  characteristics  of  the  MEPS  operation, 
personnel  flow,  and  the  data  being  entered  into  the  system. 

To  verify  the  format  of  the  elements  being  entered  into  the 
system,  the  MEPS  Headquarters  (MEPCOM)  Data  Automation  Cen- 
ter  was  contacted.  Throughout  the  study,  both  the  ARI-PERI- 
RP  and  the  USAREC  operations  research  personnel  aided  the 
analysis  with  information  and  encouragement. 

The  data  sample  received  from  ARI  was  drawn  from  the 
MEPS  Applicants  File  which  contains  690  characters  of  per¬ 
sonnel  data  on  each  individual  who  has  applied  for  testing 
and  entrance  in  any  branch  of  the  United  States  Aimed  Forces 
since  1976.  The  time  slice  of  data  that  was  made  available 
for  this  study  centered  around  personnel  entering  the  Mili¬ 
tary  Services  in  the  1980  to  1981  time  frame.  This  partic¬ 
ular  sample  consisted  of  a  file  containing  48520  records 
extracted  at  random  from  the  MEPS  file,  selecting  only  re¬ 
cords  having  a  "9"  as  the  last  digit  of  the  Social  Security 
Number. 

Selecting  Variables  that  are  Kev  to  Recruitment  Investigation 

The  first  task  involved  identifying  and  extracting  the 
"target"  elements  within  the  sample  data,  that  is,  to  select 
and  save  on  a  separate  file  those  variables  in  each  record 


that  would  be  required  to  create  the  data  "profile" ,  and 
would  later  be  used  in  the  multivariate  analysis.  The  file 
was  meticulously  screened  to  save  all  variables  relevant  to 
the  study ,  while  at  the  same  time  reducing  the  size  of  the 
records  by  eliminating  variables  that  were  either  redundant 
or  did  not  offer  clear  discriminating  categories  for  future 
investigation.  The  sample  was  examined  using  the  SPSS  sub¬ 
routines  Frequencies  and  Crosstabs  and  using  several  short 
FORTRAN  programs  in  order  to  develop  a  profile  of  the  data 
base.  The  aims  of  this  data  examination  were: 

1.  To  eliminate  those  variables  having  too  many  blank 
data  fields, 

2.  To  eliminate,  if  necessary,  records  with  invalid 
entries  critical  to  the  analysis, 

3.  To  get  a  "feel"  for  the  percentages  within  the 
categories  that  were  available  for  testing, 

4.  To  examine  the  relationships  among  variables, 

5.  To  examine  the  computed  raw  Chi -Squared  scores  among 
the  various  categories,  and 

6.  To  discover  any  obvious  patterns  existent  in  the 
cross-tabulation  tables. 

There  were  a  number  of  variables  in  the  initial  study 
of  the  MEPS  file  that  were  eliminated  from  the  final  model: 

1 •  Number -of -Dependents , 

2.  Religious-Denomination, 

3.  Years-of-Education, 

4.  Mental-Category, 
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5 •  Youth-Progr am-and-Youth-Program-Conducted-by , 

6.  Program-Enlisted-for-Option, 


7.  Designated-Option , 

8.  Enlistment -Option, 

9.  Enlistment-Option-Guaranteed, 

10.  Enlistment-Bonus -Level , 

11.  Training/Enlistment-Military-Occupational-Specialty, 

12.  Term-of -Enlistment, 

13.  Pay-Grade,  and 

14.  Citizenship-Code. 

These  apparently  useful  discriminating  variables  were 
deleted  for  the  following  reasons: 

1.  Large  numbers  of  blank  data  fields 

Religious-Denomination,  Years-of -Education,  and  Citizen¬ 
ship-Code  had  a  majority  of  uncoded  records,  rendering  most 
multivariate  analyses  useless  during  the  modeling  of  factors 
and  discriminants. 

2.  Incommensurable  codes  among  both  intra  and  inter¬ 
service  categories  of  variables  examined 

All  of  the  "Option"  code  variables  (6  through  9),  and 
Enlistment-Bonus-Level  had  a  significant  number  of  codes 
that  could  not  easily  be  measured  or  be  generally  related  to 
other  variables.  Within  Program-Enlisted-for-Option  there 
were  33  different  codes;  the  Army  codes  alone  contained  sev¬ 
eral  levels  of  Veterans  Educational  Assistance  Program  (VEAP) 
kicker  enlistment  options,  a  number  of  different  Tuition 
Assistance  choices,  several  assignment  electives,  and  many 


combinations  of  the  foregoing  alternatives.  Designated- 
Option  has  35  different  codes,  while  Enlistment-Option  has 
at  least  518  different  codes  among  which  include  many  per¬ 
mutations  of  such  things  as:  Advanced  Enlistment  Grade, 

Unit  or  Location  Guarantee,  Training  or  Skill  Guarantee,  and 
Accelerated  Promotion.  Additionally,  these  variables  have  a 
less  common  measure  or  basis  for  comparison  when  trying  to 
relate  their  numerous  permutations  between  the  different 
Services.  It  also  became  unwieldy  to  attempt  to  break  these 
given  variables  down  into  their  separate  categories  and  model 
them  as  indicator  variables.  Another  difficult  variable  to 
categorize  and  compare  within  and  among  Services  is  the 
Enlistment/Training-Military-Occupational-Specialty  code : 
the  Army  alone  has  about  300  specialty  designations. 

All  of  the  aforementioned  variables  have  a  great  deal 
of  intuitive  appeal  for  comparing  what  draws  individuals  to 
enlist  in  one  branch  of  the  Service  as'  opposed  to  another. 
Indeed,  the  kinds  of  educational  benefits,  training,  assign¬ 
ments,  advancement  potential,  and  job  opportunities  should 
be  the  very  discriminating  variables  that  would  differentiate 
individuals  from  joining  one  branch  of  the  Service  over  an¬ 
other.  However,  the  problem  boils  down  to  an  almost  impossi¬ 
ble  task  of  sorting  out,  classifying,  and  subjectively  weight¬ 
ing  the  numerous  categories  within  the  variables  involved. 
Thus,  these  particular  variables  were  dropped  from  further 
consideration. 
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3.  Frequent  occurrence  of  erroneous  data  fields 

Mental-Category  when  compared  against  the  Armed  Forces 
Qualification  Test  (AFQT)  scores  had  a  substantial  amount  of 
miscoded  data  entries.  The  Mental  Category  codes  did  not 
agree  with  the  limits  set  by  AFQT  scores  for  10.6  percent  of 
the  records.  ARI-PERI-RP  personnel  confirmed  a  problem  with 
this  data,  and  suggested  that  AFQT  was  a  more  accurate  indi¬ 
cator  of  mental  distinction.  The  "Option"  codes  also  con¬ 
tained  miscodes  of  "S"  for  "5",  zero  for  character  "O" ,  and 
"2"  for  "Z"  in  many  cases,  along  with  a  number  of  unrec  gniz- 
able  category  codes.  Also,  the  Number -of -Dependents  variable 
contained  an  erroneous  range  of  dependents  in  a  range  from 
10  to  99  in  about  two  percent  of  the  file. 

4.  Inability  to  substantially  enhance  the  predictability 
of  the  final  model 

Number-of -Dependents  may  have  had  a  tolerable  level  of 
erroneous  data,  but  it  was  decided  to  drop  it  due  to  the 
fact  that  95  percent  of  the  records  had  zero  dependents, 
attributable  to  the  high  percentage  of  young,  single  individ¬ 
uals  on  the  file.  Thus,  Number-of -Dependents  provided  very 
little  variance  to  use  for  factor  analysis  and  predictive 
purposes. 

Years-of -Education  (discounting  blank  fields)  ranged 
from  1  to  24,  with  the  distribution  of  data  approaching  the 
shape  of  a  Normal  curve,  having  a  mean  of  11.5  and  a  standard 
deviation  of  1.3.  However,  years  of  education  is  not  as  ex¬ 
act  a  measure  of  educational  achievement  as  is  the  variable 
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Level-of -Education,  since  the  "worth"  of  a  given  number  of 
years  of  education  is  only  valued  in  terms  of  the  type  of 
diploma  one  has,  not  by  how  many  years  it  took  one  to  aquire 
that  diploma.  In  addition,  the  Armed  Forces  place  much  em¬ 
phasis  on  recruiting  high  school  graduates  as  opposed  to  re¬ 
cruiting  those  who  have  twelve  years  of  education  but  no  high 
school  diploma.  And  enlistees  with  a  college  degree  have  a 
much  better  opportunity  for  enlisting  into  uniquely  qualified 
Occupational  Specialties  as  well  as  officer  training  than  do 
individuals  with  16  or  more  years  of  education  but  no  college 
diploma.  Those  who  have  completed  a  given  "level  of  educa¬ 
tion"  are  by  definition  higher  achievers  than  those  who  have 
received  as  many  years  of  education,  but  who  have  not  achieved 
the  same  level  of  education.  Also,  the  Army  prefers  to  re¬ 
cruit  high  school  graduates  over  non-high  school  graduates. 

For  these  reasons.  Year s-of -Education,  being  less  than  an 
interval  level  variable,  was  deleted  from  the  study,  and 
Level-of-Education  was  retained  in  the  study  as  a  much  more 
precise  ordinal  level  variable. 

Youth-Program-and-Youth-Program-Conducted-by  was  an 
insignificant  variable  because  of  the  very  small  number  of 
occurrences  of  individuals  involved  in  such  programs.  Term- 
of-Enlistment  was  more  like  a  dependent  variable  than  a  pre¬ 
dictor  variable;  the  term  depends  primarly  on  which  compo¬ 
nent  of  the  services  an  individual  enlists  into  (National 
Guard,  Reserves,  or  Regular),  and  on  what  options  are  select¬ 
ed  which  may  require  special  service  obligations.  Finally, 
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Pay-Grade  was  deleted  from  the  model  since  it  too  depends  on 
options  selected  and  whether  an  individual  was  prior  service; 
also,  most  records  tended  to  be  coded  "El"  Trainee. 

The  majority  of  variables  on  the  original  MEPS  file  and 
those  used  in  the  final  analyses  are  nominal  level  variables. 
Several  ordinal  level  arrangements  were  applied  without  suc¬ 
cess  such  as  ordering  Service  according  to  an  arbitrary  pro¬ 
gression  of  "desirability" ,  or  arranging  Race  according  to 
Census  Bureau  statistics  of  median  family  income.  Also,  no 
new  information  was  gained  from  assigning  interval  values  to 
ordinal  level  data  such  as  treating  Level-of -Education  cate¬ 
gories  as  integers.  Since  none  of  these  variable  alterations 
were  fruitful,  it  was  decided  to  recode  each  of  the  given 
variable  categories  as  "dummy"  (indicator)  variables. 

Index  Variables:  Time  and  Location 

Early  in  the  investigation  a  need  was  recognized  to  per- 

\ 

form  regional  analysis  on  the  MEPS  file  as  well  as  to  manip¬ 
ulate  time  variables  in  the  data.  Part  of  the  goal  of  this 
study  was  to  determine  whether  there  existed  a  variance  in 
military  recruiting  depending  on  time  and  geographical  loca¬ 
tion.  Both  the  time  and  the  location  variables  presented 
difficulty  in  early  attempts  to  correctly  isolate  and  cate¬ 
gorize  their  values  and  the  variables  related  to  them.  There 
are  a  number  of  dates  on  the  MEPS  file  associated  with  every 
record  to  include: 


1.  Date-of -Action , 

2.  Date-of-Birth, 


3.  Date-of -Entry , 

4.  Date-of -Entry-2, 

5 .  Entry-or-Discharge-Date , 

6.  Date-of -Grade, 

7.  Projected-Active-Duty-Date, 

8.  Delayed-Entry-Program-Date-of -Entry, 

9 .  Advanced-Individual-Training-Graduation-Date , 

10.  Date-Processing-Took-Place,  and 

11.  Cycle -Number -Julian -Date. 

The  last  two  dates  (items  10  and  11)  each  have  up  to 
ten  occurrences  depending  on  the  number  of  transactions  pro¬ 
cessed  against  any  given  individual,  such  as  mental  examina¬ 
tion  on  one  day,  physical  examination  on  another  day,  entry 
into  Delayed  Entry  Program  (DEP)  another  day,  and  possible 
enlistment  into  the  Regular  force  at  a  later  date.  Based  on 
the  Date-of -Entry  variable  in  the  original  MEPS  data  sample, 
37.9  percent  of  the  records  belong  to  1980,  while  54.8  per¬ 
cent  belong  to  1981,  and  7.3  percent  are  blank.  Based  on  the 
Date-of-Action  variable,  60.7  percent  of  the  records  belong 
to  1980,  while  39.3  percent  belong  to  1981.  Using  the  earli¬ 
est  of  the  Date-Processing-Took-Place  entries  in  each  record, 
4.2  percent  of  the  records  belong  to  1979,  while  78.9  percent 
belong  to  1980,  with  16.9  percent  in  1981. 

Choosing  the  correct  date  was  important  to  the  profile 
building  of  the  initial  investigation,  but  it  became  even 
more  crucial  for  the  final  analyses  when  the  MEPS  and  USAREC 
files  were  unified  in  a  study  of  the  key  variables  from  both 


files.  It  became  imperative  in  that  final  study  for  the  vari¬ 
ables  to  be  combined  into  a  meaningful  sequence.  ARI-PERI-RP 
faced  the  same  problem  of  arriving  at  a  date  variable  that 
would  properly  represent  prospective  accessions  at  a  definite 
point  in  time. 

ARI-PERI-RP  representatives  visited  a  Military  Entrance 
Processing  Station  (MEPS)  to  examine  the  flow  of  personnel 
through  the  system  and  to  understand  the  dates  associated 
with  the  various  transactions  on  the  file.  These  representa¬ 
tives  determined  that  using  Date-of -Action,  Date-of -Entry , 
or  a  combination  of  the  two  was  not  capturing  the  time  clos¬ 
est  to  the  prospective  enlistee's  decision  to  join  the  Ser¬ 
vice.  Rather  they  suggested,  following  their  thorough  re¬ 
search,  that  the  best  date  to  use  for  this  purpose  was  the 
earliest  entry  from  Date-Processing-Took-Place.  This  pre¬ 
ferred  date  was  used  in  all  subsequent  analyses  in  this  study, 

\ 

to  include  the  crucial  MEPS  and  USAREC  file  merging  phase. 

The  location  variable  was  not  as  elusive  as  the  time 
variable,  but  choosing  workable  regional  boundaries  proved 
to  be  virtually  impossible  given  the  original  regional  vari¬ 
ables  on  the  MEPS  file.  The  MEPS  data  base  has  the  following 
location  variables  within  each  individual  record: 

1.  Sector-ID  (Eastern,  Central,  and  Western), 

2.  MEPS-ID , 

3.  Home-of -Record,  state  and  county, 

4.  Home-of -Record,  ZIP  Code, 

5.  Present-Address,  state  and  county  , 
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6.  Present-Address,  ZIP  Code, 

7.  Transfer-to-Code  (unit  of  assignment), 

8.  MEPS-Code  (up  to  ten  transactions), 

9.  Recruiting-Station-Identification, 

10.  Advanced-Individual-Training-Location,  and 

11.  Test-Site. 

One  location  variable  used  by  ARI  to  generate  regionality 
within  their  data  is  the  MEPS-ID.  There  are- 71  MEPS-ID's, 
each  representing  one  of  the  processing  centers  where  the 
potential  enlistees  are  examined  and  inducted  into  the  vari¬ 
ous  branches  of  the  Military  Service.  The  problem  with  using 
this  location  variable  is  that  it  cannot  be  directly  related 
to  the  Recruiting  Regions  and  Districts  that  are  the  locus 
of  Army  recruiting  efforts  as  well  as  the  locus  of  the  acces¬ 
sion  results  of  that  recruiting.  Compounding  this  problem 
is  the  fact  that  a  group  of  individuals  may  all  be  recruited 
from  the  same  ZIP  Code  and  county  but  -be  processed  separately 
at  several  different  MEPS.  Thus,  there  is  no  homogeneity 
for  developing  characteristics  of  population  based  on  MEPS-ID. 
Recruiting-Station-ID  cannot  be  used  since  it  is  not  recorded 
on  many  of  the  records  in  the  file.  In  addition,  for  those 
records  in  which  Recruiting-Station-ID  does  appear,  it  is 
coded  differently  for  each  branch  of  the  Service,  and  even 
the  Army  codes  do  not  directly  represent  the  five  USAREC 
Regions  or  their  57  District  Recruiting  Commands  ( DRC ) . 

Initial  regional  analysis  of  the  file  had  to  be  based 
on  MEPS-ID,  but  breaking  the  country  down  into  these  71 


processing  locations  was  fruitless,  especially  given  the  nom 
inal  nature  of  roost  of  the  remaining  MEPS  variables,  and 
given  the  fact  that  the  MEPS-ID's  could  not  be  directly  re¬ 
lated  to  recruiting  results  of  any  given  USAREC  Region,  In 
addition,  any  kind  of  multivariate  analysis  obtained  from 
this  type  of  regional  breakdown  could  not  easily  be  related 
to  analysis  of  the  USAREC  data  base. 

Determining  How  to  Link  the  MEPS  and  USAREC  Files 

The  USAREC  data  base  is  organized  by  DRC  within  year 
and  Quarter,  thus  giving  each  year  a  total  of  228  records 
(57  DRC's  times  4  Quarters).  Each  of  these  records  contains 
summary  statistics  for  the  individual  DRC's  concerning  demo¬ 
graphics,  economics,  recruiter  data,  advertising,  and  acces¬ 
sions.  The  only  time  variable  on  the  USAREC  file  is  Fiscal- 
Year-and-Quarter.  The  only  location  variable  is  DRC;  yet 
there  is  no  regional  variable  on  the  MEPS  file  that  is  di¬ 
rectly  comparable  or  easily  translatable  into  DRC.  ARI  was 
not  aware  of  a  conversion  technique,  but  gladly  welcomed  our 
support  in  an  effort  to  accomplish  such  a  translation  since 
several  of  their  studies  could  benefit  from  having  the  abil¬ 
ity  to  analyze  the  MEPS  file  by  DRC.  This  was  incentive 
enough  in  itself  to  launch  an  effort  to  translate  DRC's  for 
use  in  merging  the  two  files  and  to  eventually  analyze  the 
combined  data  base. 


The  first  stage  in  DRC  translation  was  to  examine  all 
prospective  location  variables,  testing  for  adaptability  in 
conversion  and  for  completeness  of  record  entries.  That  is, 


any  adopted  conversion  variable  must  be  able  to  be  precisely 
translated  into  a  DRC  code,  and  information  must  be  consis¬ 
tently  recorded  in  the  conversion  variable's  data  fields, 
meaning  that  there  must  be  few,  if  any,  blank  entries.  One 
thing  to  keep  in  mind  in  this  research  is  that  many  MEPS 
data  fields  are  intentionally  left  blank  in  order  to  econo¬ 
mize  coding  effort,  depending  on  the  status  of  the  individual 
at  various  stages  in  the  testing-to-enlistment  process.  ZIP 
Codes  were  studied  for  a  time  until  it  was  realized  that  ZIP 
Codes  do  not  correspond  to  county  boundaries,  sometimes  having 
a  ZIP  Code  overlap  different  counties.  Also,  both  USAREC 
and  the  Census  Bureau  organize  their  statistics  along  county 
boundaries,  not  by  ZIP  Code.  USAREC-PAE-RE  provided  a  map 
of  these  county  boundaries  by  DRC  and  Region,  which  could  be 
used  to  recode  the  MEPS  state  and  county  codes  into  DRCs. 

The  applicable  location  variables  for  state/county  on  the 
MEPS  file  are:  Home-of -Record,  state  and  county,  and  Present- 
Address,  state  and  county.  Examination  of  the  data  within 
these  fields  revealed  that  Home-of -Record  had  much  more  con¬ 
sistency  in  information  being  present  than  did  Present-Address. 
However,  it  was  felt  that  in  those  cases  where  Present-Address 
had  valid  information  different  from  Home-of -Record,  Present- 
Address  should  then  be  used.  The  logic  in  this  is  that 
Present-Address  is  a  better  indicator  of  the  location  from 
which  an  individual  is  recruited.  With  the  proper  state  and 
county  determined,  the  next  step  was  to  relate  each  county 
to  its  proper  DRC. 
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DRC's  do  not  overlap  county  borders,  however,  DRC's  do 
overlap  state  borders,  that  is,  any  given  DRC  may  contain 
counties  from  more  than  one  state,  and  the  respective  states 
may  contain  parts  of  several  DRC's,  This  means  that  a  county 
to  DRC  conversion  must  be  carried  out  on  a  county  by  county 
b^isis.  Each  county,  all  United  States  possessions,  and  cer- 
tain  foreign  countries  were  recorded  by  their  proper  DRC  and 
then  coded  into  a  FORTRAN  program.  This  program  was  used  to 
convert  the  county/states  to  DRC's  while  also  converting  and 
purifying  several  other  variables  on  the  MEPS  data  base,  such 
as  Level-of -Education ,  Race,  Marital-Status,  Service,  and 
Sex,  as  well  as  creating  or  computing  others  such  as  Age, 
Hispanic,  and  a  composite  code  for  Status.  The  output  gen¬ 
erated  by  this  program  allows  ARI  for  the  first  time  to  anal¬ 
yze  the  MEPS  files  using  DRC  codes.  This  is  a  capability 
beneficial  to  both  ARI  and  USAREC  since  it  facilitiates  the 
further  exchange  of  information  between  the  MEPS  and  USAREC 
data  bases. 


Proaram  for  Merqinq  the  Data  Bases 

The  groundwork  having  been  laid  for  analyzing  the  MEPS 
file  based  on  DRC  and  for  relating  the  two  data  bases,  the 
next  task  was  to  create  a  program  that  would  do  the  actual 
merging  of  the  MEPS  and  USAREC  files.  Since  the  USAREC  file 
is  composed  of  summary  statistics  organized  by  Quarter,  while 
the  MEPS  file  is  organized  as  individual  personnel  records, 
there  are  two  approaches  that  can  be  used  to  merge  the  files: 
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1.  Condense  the  MBPS  data  into  summary  format  by  DRC 
by  Quarter ,  or 

2.  Attach  applicable  USAREC  data  fields  to  each  MEPS  re¬ 
cord  by  individual  by  DRC  by  month  from  the  applicable  Quarter. 
This  given  approach  was  taken  since  the  primary  concern  was 

to  analyze  recruiting  variables  as  they  affected  individuals. 

The  USAREC  file  contains  chiefly  ratio  level  data  com¬ 
posed  of  demographic ,  economic,  and  accession  statistics. 

Many  of  the  variables  pertain  only  to  the  Army  and  thus  can¬ 
not  be  directly  applied  against  the  Navy,  Air  Force,  or  Marine 
records  in  the  MEPS  file.  Among  the  variables  not  applied  or 
merged  with  the  MEPS  file  are:  Percent-of-DRC-Army-Recruiters- 
with-One-or-More-Years-Experience ,  Army-Non-Prior-Service- 
Female-Accessions ,  and  DRC-Army-Local-Advertising-Expenditures. 
However,  these  variables  and  others  in  the  USAREC  file  were 
considered  and  used  in  the  improvement  of  the  USAREC  Regres¬ 
sion  Model.  Seven  other  variables  werfe  entered  into  the  mul¬ 
tivariate  factor  and  discriminant  analyses  since  these  seven 
variables  pertained  directly  to  all  of  the  Armed  Services: 

1.  Qualified-Military-Available, 

2.  Male-High-School-Seniors, 

3.  DRC-Overall-Unemployment , 

4.  DOD-Non-Prior-Service-Male-Accessions , 

5 .  DOD-High-School-Diploma-Gr aduate-Male-Accessions , 

6.  Median-Disposable-Family-Income,  and 

7.  Black-Military-Available. 


There  is  an  eighth  Army  variable  on  the  USAREC  file  that 
has  a  direct  inverse  relationship  with  the  other  branches  of 
Service: 

8.  Army-Recruiters-as-Percent-of-DOD-Recruiters-by-DRC. 

The  information  from  these  eight  variables  was  merged 
with  the  data  from  the  MEPS  file  based  on  the  MEPS  conversions 
of  state/county  to  DRC  and  the  conversion  of  Month  to  Quarter. 
The  program  that  was  written  to  do  this  merging  function  can 
merge  a  period  of  up  to  three  years  of  data  at  a  time,  and 
can  easily  be  modified  to  merge  more  years  on  a  given  run 
or  modified  and  rerun  several  times  to  merge  data  for  a  . 
group  of  years.  This  program  is  the  final  in  a  series  of 
three  programs  that  can  be  used  by  anyone  having  access  to 
both  the  MEPS  file  and  the  USAREC  file.  These  programs  are 
written  in  FORTRAN  so  that  the  input  end  output;  formats  can 
be  easily  modified  to  incorporate  more  or  fewer  of  a^y  of 
the  variables  located  on  either  file.  Thus,  expansion  or 
contraction  of  the  numbers  and  kinds  of  variables  to  be  anal¬ 
yzed  is  permitted  with  little  recoding  effort. 

An  option  also  exists  to  use  only  the  MEPS  file  by  run¬ 
ning  only  the  first  two  programs  in  the  series#  Thus,  analysis 
of  the  MEPS  data  base  can  be  done  independently  of  the  USAREC 
data  and  USAREC  variables,  but  with  the  benefit  of  DRC  recod¬ 
ing.  This  translated  data  gives  ARI  the  capability  of  per¬ 
forming  several  of  the  DRC  analyses  that  they  could  not  do 
in  the  past.  Output  from  both  the  second  and  third  programs 
is  highly  compatible  with  SPSS  input  format  and  the  data 


manipulation  requirements  since  it  is  almost  all  numeric  in 
format,  unlike  the  unrecoded  MEPS  file. 

The  output  of  the  third  program  is  a  merged  MEPS  USAREC 
(MU)  file  containing  the  following  variables: 

1.  Year, 

2.  Month, 

3.  DRC  (first  character  is  a  Region  code:  1,3, 4, 5, or  6), 

4.  Status,  a  composite  variable  based  on  the  MEPS  Status 
Code  for  an  individual  record  and  on  the  MEPS-Entry  Status. 

See  Appendix  a7  for  an  explanation  of  these  variables.  The 
following  are  the  MU  Status  codes  and  their  parent  codes  from 
the  MEPS  file: 

1  =  Qualified,  Not  Joined  (MEPS  Status-Code  D,  MEPS 

Entry-Status  blank ) , 

2  =  Enlisted  into  Delayed  Entry  Program  ( DEP )  (MEPS 

Status-Code  A,  MEPS  Entry-Status  3), 

3  =  Enlisted  without  DEP  (MEPS  Status-Code  A,  MEPS 

Entry-Status  0) , 

4  *  Shipped  Reservist  from  DEP  (MEPS  Status-Code  B 

or  C,  MEPS  Entry-Status  2), 

5  =  Shipped  Enlistment  from  DEP  (MEPS  Status  Code 

B  or  C,  MEPS  Entry-Status  1), 

6  =  Shipped  without  DEP  (MEPS  Status-Code  B  or  C, 

MEPS  Entry  Status  0), 

0  *  all  other  combinations  of  MEPS  Status -Code  and 

MEPS  Entry-Status. 

5.  Service  (Army,  Navy,  Air  Force,  or  Marines), 


6.  Component  (National  Guard,  Reserve,  or  Regular), 

7.  Prior-Service  Indicator, 

8.  Mental-Category  (for  reference  only), 

9.  Armed-Forces  Qualification-Test  (AFQT)  score, 

10.  Level-of-Education  (Less  Than  HS  Diploma,  HS  Equiv¬ 
alency,  HS  Senior,  1  year  College,  Associates  or 
Nurse,  Baccalaureate,  Masters,  Post  Masters,  or 
Doctorate/Professional) , 

11.  Age  (computed  using  MEPS  Date-of -Birth  and  earliest 
Date-Processing-Took-Place ) , 

12.  Sex, 

13.  Race, 

14.  Hispanic  (recoded  from  MEPS  Ethnic-Group), 

15.  Marital-Status  (Married,  Never  Married,  or  Other), 

16.  Qualified-Military-Available, 

17.  Male-High-School-Seniors, 

18.  DRC-Overall-Unemployraent , 

19.  Army-Recruiters-As-Percent-of-DOD-Recruiters-by-DRC 

20.  DOD-Non-Prior-Service-Male-Accessions , 

21.  DOD-High-School-Diploma-Graduate-Male-Accessions , 

22.  Median-Disposable-Family-Income,  and 

23.  Black-Military-Available. 

These  23  variables  became  the  heart  of  all  further  re¬ 
search  in  the  profile  and  multivariate  analyses.  The  MU  file 
merges  the  data  base  resources  of  ARI  and  USAREC,  and  can  be 
used  as  a  model  for  exchange  of  data  resources  between  these 
organizations  through  the  common  link  of  DRC  codes.  Listings 


Will 


Several  SPSS  Cross-tabulation  runs  were  carried  out  on 
the  various  MU  variables,  breaking  the  data  out  mainly  by 
Service,  by  Region,  and  by  Level-of -Education.  The  ARI-PERI 
RP  operations  researchers  encouraged  this  type  of  investiga¬ 
tion  since  it  would  help  in  focusing  on  the  key  variables 
and  the  categories  within  variables  to  be  evaluated.  The 
researchers  at  ARI  periodically  examine  the  MEPS  data  base 
in  this  way  in  response  to  the  many  studies  undertaken  by 
their  organization. 

Based  on  the  MEPS  earliest  Date-Processing-Took-Place 
entry  that  was  used  as  the  Month  variable  in  the  new  MU  file 
only  4.2  percent  of  the  data  fell  within  1979,  with  70.5 
percent  of  this  being  skewed  to  the  last  four  months.  The 
1980  data  was  also  skewed  with  the  first  four  months  con¬ 


taining  only  10.0  percent  of  the  1980  data.  The  1981  data 
included  only  January  and  February  records.  Some  initial 
analysis  was  conducted  on  the  1980  records  alone,  while  all 
later  investigation  down  to  the  DRC  level  was  conducted  on 
the  latest  one  year  period  from  March  1980  to  February  1981. 
For  all  practical  purposes  this  means  that  only  four  quar¬ 
ters  of  USAREC  data  were  involved  in  the  study:  from  the 
second  quarter  of  1980  to  the  first  quarter  of  1981.  The 
month  of  March  in  the  first  quarter  of  1980  contains  only 
2.1  percent  of  the  last  12  months  of  data. 
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The  Year  variable  was  combined  into  the  Month  variable 
by  recoding  January  1981  as  month  13  and  February  1981  as 
month  14.  The  MU  Status  variable  was  not  considered  in  the 
final  factor  and  discriminant  analyses  since  it  is  more  of  a 
dependent  variable  than  a  predictor  of  whether  an  individual 
will  join  the  Army  over  another  branch  of  service.  Under 
factor  analysis  the  MU  Service  variable  was  recorded  into 
four  indicator  variables,  each  representing  a  branch  of  the 
Service.  Under  discriminant  analysis  the  MU  Service  variable 
was  recoded  into  groups  as  Army  and  Non-Army.  The  MU  Compo¬ 
nent  variable  was  used  to  break  the  study  down  into  three 
separate  analyses,  so  that  the  factors  affecting  recruiting 
of  National  Guard,  Reserves,  and  Regular  enlistees  could  be 
investigated  separately. 

At  the  DRC  level,  the  MU  Prior-Service  variable  was  used 
to  examine  only  those  without  prior  service  since  there  are 
certain  Service  unique  restrictions  imposed  on  prior  service 
enlistees.  The  Air  Force,  for  instance,  will  not  accept 
prior  service  individuals  unless  they  enlist  into  critical 
skills.  This  fact  tends  to  "push"  unemployed  prior  service 
individuals  toward  enlistment  into  the  Army,  as  is  borne  out 
in  the  initial  analysis  by  region. 

Additionally,  the  Prior-Service  variable  was  used  to 
delete  prior  service  records  from  the  final  model  because 
all  branches  of  the  Service  are  interested  primarily  in  those 
factors  that  affect  the  recruiting  of  non-prior  service  male 
accessions.  For  this  same  reason  the  variable  Sex  was  used 


to  delete  all  female  records  from  the  final  DRC  level  anal¬ 
yses;  the  Services  do  not  have  to  gear  advertising  or  re¬ 
sources  toward  the  recruitment  of  females  because  the  quotas 
for  female  enlistments  are  easily  being  filled  without  spe¬ 
cial  recuitment  efforts  in  that  area. 

Records  with  blank  Mental  Category  fields  were  elimi¬ 
nated  from  consideration  in  the  analyses.  The  MU  Level-of- 
Education  variable  was  recoded  into  indicator  variables  as 
indicated  in  the  variable  list  below.  Race  was  also  recoded 
to  allow  for  indicator  variables  in  both  White  (75.9  percent 
of  MU  file)  and  Black  (19.8  percent  of  MU  file)  records. 
Finally,  Married  was  established  as  a  separate  indicator 
variable  based  on  the  Marital-Status  variable. 

The  data  on  the  MU  file  was  not  altered,  rather  certain 
variables  were  recoded  using  SPSS  "IF"  statements  to  create 
"dummy"  (indicator)  variables,  while  certain  other  variables 
ceased  to  be  considered  in  the  final  analyses  at  DRC  level. 
The  following  is  a  list  of  all  the  recoded  variables  used  as 
input  to  the  multivariate  analyses: 

1 .  Month , 

2.  DRC, 

3.  Army, 

4.  Navy, 

5.  Air-Force, 

6.  Marines, 

7.  AFQT  score, 

8 .  No-High-School-Diploma , 


9.  High-School-Senior, 

10.  High-School-Diploma,  hut  less  than  college  degree, 

11 .  College-Degree , 

12.  Age, 

13.  White, 

14.  Black, 

15.  Hispanic, 

16.  Married, 

17  -  24.  the  eight  USAREC  variables. 

Once  the  key  variables  had  been  selected  and  the  data 
base  constructed,  the  next  phase  of  the  study  could  be  per¬ 
formed,  namely  the  multivariate  classification  analyses: 
Factor  Analysis  and  Discriminant  Analysis  using  SPSS  sub¬ 
routines. 


VI.  Multivariate  Analyses 


Factor  Analysis 

The  factor  relationships  among  the  variables  were  stud¬ 
ied  at  the  national,  regional,  and  DRC  levels  to  focus  on 
and  uncover  any  tangible  demographic  differences  among  pro¬ 
spective  recruits.  Armed  with  this  data  profile  and  the 
related  statistics ,  the  planners  at  USAREC  might  better  be 
able  to  channel  resources  to  meet  the  needs  of  a  particular 
region  and  DRC  based  on  known  critical  factors  influencing 
a  particular  segment  of  the  nation.  Understanding  where 
there  do  and  do  not  exist  differences  in  a  particular  vari¬ 
able's  relationship  with  other  variables  and  what  those  re¬ 
lationships  are  can  often  give  the  decision  maker  a  key  to 
the  realignment  of  scarce  resources. 

Initial  analysis  and  the  cross-tabulation  runs  showed 
that  there  existed  a  substantial  difference  in  the  data  by 
Component  (National  Guard,  Reserves,  and  Regular).  There 
were  no  Navy  or  Marine  National  Guard  and  very  few  Air  Force 
Reserve  records  on  the  MBPS  file.  Also,  the  different  com¬ 
ponents  were  logically  best  analyzed  in  separate  categories 
since  the  Regular  forces  have  a  different  mission  than  the 
Guard  or  Reserves.  When  the  forces  were  broken  down  by 
component,  there  were  not  enough  record  samples  in  the  Guard 
to  do  separate  analysis  of  this  component.  In  the  Reserve 
there  were  only  enough  records  to  carry  out  analysis  down 
to  the  regional  level.  There  were,  however,  enough  Regular 
forces  records  so  that  both  factor  and  discriminant  analysis 


could  be  performed  down  to  the  DRC  level  on  the  Regular 
component . 


Regional  level  analysis  was  conducted  both  with  and 
without  prior  service  and  female  records,  while  all  analysis 
at  the  DRC  level  was  performed  only  on  male  non-prior  ser¬ 
vice  records. 


How  to  Read  the  Factor  Analysis  Tables 

The  factor  analysis  results  are  given  in  tabular  form 
by  Region  and  by  DRC.  Explanations  of  variables  and  codes 
for  multivariate  analysis  are  given  in  Appendix  C-l.  The 
complete  set  of  Factor  Tables  is  given  in  Appendices  C-2 
through  C-4: 

C-2  =  Factor  Analysis  Tables  by  DRC, 

C-3  =  Factor  Analysis  Table  by  Region  for  Regular 
Forces ,  and 

C-4  =  Factor  Analysis  Table  by  Region  for  Reserve  Forces. 

Appendix  C-2  should  be  referred  to  while  reading  the 
following  explanation  of  the  Factor  Tables.  The  various 
principal  components  (factors)  are  listed  across  the  top  of 
the  table.  In  the  upper  half  of  the  table  are  listed  the 
variables  (by  row,  e.g..  Army,  Navy)  which  make  up  the  vari¬ 
ous  factors.  In  the  bottom  of  the  table  are  the  different 
DRC's  for  a  given  Region  (e.g.,  1A,  IB):  and  included  on  the 
left  hand  side  is  a  column  of  numbers  representing  the  per¬ 
centage  of  information  that  the  factors  have  retained  from 
the  original  variables.  Interpretation  of  the  factors  is 


accomplished  by  knowing  from  what  variables  the  given  factor 
is  derived  and  recognizing  how  much  weight  a  given  variable 
contributes  to  the  composition  of  the  given  factor. 

These  variable  weights  are  based  on  the  values  of  the 
factor  loadings,  which  are  given  as  matrix  output  from  SPSS. 
Signs  for  the  factor  loadings  are  important  in  explaining 
positive  or  negative  changes  in  the  value  of  the  factor. 

The  factor  value  will  tend  to  increase  when  variables  with 
which  it  has  high  positive  correlations  increase,  and  de¬ 
crease  when  negatively  correlated  variables  increase  in 
value.  Thus  the  positive  and  negative  correlations  of  the 
given  factor  can  be  thought  of  as  positive  and  negative 
"poles"  of  the  factor,  and  the  correlated  variables  are 
associated  with  the  positive  or  negative  "pole"  depending 
on  their  correlation  sign.  Since  values  in  the  factor  load¬ 
ings  matrix  are  correlations ,  they  must  be  squared  to  deter¬ 
mine  the  percentage  of  contribution  they  provide  to  the 
explanation  of  any  given  factor.  The  tabular  output  results 
take  this  into  consideration  by  assigning  "measures"  within 
the  cells  of  the  Variable  by  Factor  part  of  the  table  as 
follows: 

"++"  =  50  to  99  percent  of  the  variable  is  like  the 

positive  "pole"  of  the  factor* 

" — ••  -  50  to  99  percent  of  the  variable  is  like  the 

negative  "pole"  of  the  factor, 

"+"  =  25  to  50  percent  of  the  variable  is  like  the 

positive  "pole"  of  the  factor, 


5 


=  25  to  50  percent  of  the  variable  is  like  the 

negative  “pole"  of  the  factor. 

The  table  is  intended  to  be  read  using  the  following 
algorithm: 

1.  Select  a  DRC  of  interest. 

2.  Read  across  the  selected  DRC  row  until  an  "X"  is 
encountered.  This  signifies  that  a  given  factor  in  the 
column  above  is  associated  with  the  selected  DRC. 

3.  Read  up  the  "X"  column  until  a  " — ",  ,  or 

is  encountered.  This  signifies  that  a  given  variable 

is  associated  with  the  given  factor. 

4.  Read  across  the  row  to  the  left  to  find  the  asso¬ 
ciated  variable  and  note  the  variable  name. 

5.  Go  back  along  the  row  to  the  given  factor  column. 

6.  Repeat  steps  3  through  5  until  the  top  of  the  table 
is  reached. 

7.  Record  a  preliminary  interpretation  for  the  combi¬ 
nation  of  the  noted  variables  and  call  this  interpretation 
a  factor. 

8.  Go  back  down  the  column  to  the  original  DRC  row. 

9.  Repeat  steps  2  through  8  until  the  right  hand  side 
of  the  table  is  reached. 

10.  Analyze  the  recorded  factors  among  themselves  for 
a  within  DRC  comparison. 

11.  Repeat  steps  1  through  10  until  all  individual 


DRC's  of  interest  are  analyzed. 


12.  Analyze  all  DRC's  of  interest  between  DRC's. 
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With  these  Factor  Tables  one  can  compare  DRC's  with 
each  other  or  by  comparing  one  table  with  another  table, 
compare  different  Regions  with  each  other.  Thus,  one  can 
quickly  refer  to  the  factors  facing  recruiters  within  any 
given  DRC  and  then  make  comparisons  by  DRC  or  Region  to  see 
what  factors  other  recruiters  face.  At  the  managerial  level, 
one  may  be  able  to  use  the  tables  to  aid  in  the  allocation 
of  scarce  resources  depending  on  the  profile  characteristics 
of  a  given  Region  or  DRC. 

Factor  Analysis  Interpretation 

The  first  DRC  in  Appendix  C-2  is  DRC  1A  (Albany)  with 
associated  factors  in  columns  1,  10,  15,  17,  20,  26,  32,  and 
39.  For  convenience,  column  1  will  be  referred  to  as  "factor 
1" ,  column  2  will  be  referred  to  as  "factor  2",  and  so  on. 
Factor  1  is  composed  of  a  strong  increase  in  accessions  ac¬ 
companied  by  a  strong  increase  in  unemployment  rates,  with 
a  slight  positive  correlation  associated  with  an  increase  in 
the  number  of  Army  recruiters.  Thus,  one  can  interpret  fac¬ 
tor  1  as  a  barometer  for  accessions  which  will  increase  or 
decrease  in  the  same  direction  as  the  increase  or  decrease 
in  the  percentage  of  unemployment,  and  slightly  "influenced" 
by  movement  in  the  same  direction  of  a  change  in  the  percent¬ 
age  of  Army  recruiters  in  the  DRC.  It  should  be  noted  at 
this  point  that  only  DRC's  1A  (Albany),  1G  (New  Haven),  and 
IN  (Syracuse)  have  this  particular  relationship  within 
Region  1.  However,  DRC  ID  (Concord)  in  factor  2  has  a 
stronger  correlation  between  the  increase  in  percentage  of 


Army  recruiters  and  the  increase  in  accessions.  Continuing 
in  this  vein,  factor  3  shows  that  the  other  DRC's  in  Region 
1  have  no  noticeable  link  between  an  increase  in  percentage 
of  Army  recruiters  and  the  increase  in  accessions.  However, 
factor  3  still  strongly  correlates  increases  in  unemployment 
with  an  increase  in  accessions  in  all  of  the  DRC's  in  Re¬ 
gion  1.  As  one  reads  across  to  see  what  factors  are  present 
in  a  given  DRC ,  one  should  freely  perform  a  simultaneous 
analysis  of  related  factors  in  other  DRC's  in  order  to 
quickly  compare  how  the  DRC's  differ  among  one  another. 

This  gives  the  user  a  better  grasp  of  the  distribution  of 
factor  influences  within  the  entire  Region;  and  it  is  an 
aid  to  understanding  the  structure  of  a  given  factor  in  a 
particular  DRC  by  relating  it  to  similar  factors  in  other 
DRC's. 

Factor  10  in  DRC  1A  (Albany)  shows  that  as  time  in¬ 
creases  there  is  an. accompanying  increase  in  the  median 
family  income  as  well  as  an  increase  in  the  number  of  male 
high  school  seniors  in  the  given  DRC.  One  should  keep  in 
mind  the  period  over  which  this  analysis  was  performed,  that 
is,  March  1980  to  February  1981.  Only  DRC  1A  (Albany)  and 
DRC  1G  (New  Haven)  have  the  factor  10  relationship,  while 
in  DRC  IE  (Harrisburg)  factor  9  shows  a  strong  negative 
correlation  between  number  of  male  high  school  seniors 
versus  the  month  and  income  variables.  That  is,  in  DRC  IE 
(Harrisburg)  as  months  and  income  increase  the  number  of 
seniors  decreases.  Another  closely  associated  factor  is 


factor  8  in  DRC  IB  (Baltimore)  and  in  DRC  IP  (Port  Monmouth), 
which  is  really  an  enhancement  of  factor  10  in  that  the  per¬ 
centage  of  Army  recruiters  strongly  moves  in  the  same  direc¬ 
tion  with  months,  income  and  the  number  of  male  high  school 
seniors.  In  this  same  manner,  factor  7  in  DRC  1L  (Pitts¬ 
burgh)  is  an  extension  of  factor  9,  showing  that  an  increase 
in  percentage  of  Army  recruiters  is  correlated  with  a  de¬ 
crease  in  the  number  of  seniors.  Other  factors  reflecting 
the  month  and  income  variable  influences  are  factors  4,  5, 

6,  and  11. 

Factor  15  in  DRC  1A  (Albany)  is  the  race  factor  which 
is  a  relationship  that  is  present  in  31  of  the  57  DRC’s. 
Factor  16  is  another  race  factor  but  it  contains  the  slight 
influence  of  AFQT.  Factor  17,  college  degree  is  a  separate 
ingredient  or  principal  component  in  DRC  1A  (Albany).  Fac¬ 
tor  17  is  separate  in  most  of  the  DRC's,  however  it  is  some¬ 
times  correlated  with  other  variables  such  as  with  age  in 
DRC  1G  (New  Haven)  and  in  DRC  1L  (Pittsburgh).  Factor  20 
in  DRC  1A  (Albany)  shows  a  dichotomy  between  high  school 
senior  recruits  and  high  school  diploma  recruits  with  a 
strong  linkage  to  the  age  of  the  recruit.  Factor  26,  the 
Army  variable  correlated  with  the  No-High-School-Diploma 
variable,  is  also  evident  in  DRC  1A  (Albany).  This  factor 
implies  that  non-high  school  diploma  recruits  tend  to  join 
the  Army  in  DRC  1A  (Albany).  Other  factors  in  this  DRC 
include  the  Marines  as  a  factor  by  themselves.  Also,  the 
Navy  and  Air  Force  are  one  category,  although  the  signs 
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between  these  Service  variables  are  opposite.  The  Air  Force 
versus  Navy  factor  by  itself  simply  implies  that  recruits 
from  these  two  services  have  virtually  opposite  association 
among  the  other  variables  in  the  sample  from  DRC  1A  (Albany). 
That  is,  one  can  distinguish  the  recruits  between  these  two 
services  based  on  the  other  variables  in  the  analysis. 

As  a  cautionary  note,  the  preceding  interpretation  of 
factors  must  be  combined  with  all  the  other  facts  available 
concerning  the  DRC  for  the  given  period  of  analysis.  Common 
sense  and  good  management  are  essential  in  any  use  of  these 
findings.  Also,  exact  interpretation  of  the  factors  is  not 
always  possible.  All  factors  in  a  given  DRC  are  orthogonal, 
that  is,  each  factor  is  independent  of  all  the  other  factors. 
Possibly  the  safest  way  to  handle  comparisons  of  large  num¬ 
bers  of  factors  is  to  use  the  technique  of  hyphenation  be¬ 
tween  the  variables  for  a  given  factor  to  form  a  chain  of 
variables  naming  the  given  factor.  This  is  possibly  a  more 
precise  method  of  carrying  meaning,  but  it  lacks  interpre¬ 
tation  and  ease  of  handling.  Also  to  be  noted  is  the  per¬ 
centage  given  in  the  left  hand  column  in  the  bottom  half  of 
the  table.  As  an  example,  the  percentage  for  DRC  1A  is  only 
78,  which  means  that  22  percent  of  the  information  contained 
in  the  original  variables  is  absent  from  the  factors  shown 
in  the  table  for  that  DRC.  These  cautionary  notes  must  be 
kept  in  mind  when  performing  analysis  on  demographic  data 
that  is  hard  to  measure  and  prone  to  subjective  interpreta¬ 
tion. 
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The  analysis  of  factors  points  out  the  principal  com¬ 
ponents  that  recruiters  face  in  a  given  DRC.  At  the  same 
time,  relationships  are  examined  between  DRC's  and  between 
similar  factors  as  an  aid  to  the  understanding  of  the  prin¬ 
cipal  components  which  affect  the  given  DRC.  This  same 
procedure  can  be  applied  to  regional  level  analysis  as  well 
Regional  level  managers  may  be  able  to  use  the  factor  anal¬ 
ysis  to  redirect  resources  through  an  understanding  of  the 
differences  between  DRC's  based  on  the  different  factors 
affecting  the  individual  DRC's.  It  is  important  that  these 
factors  be  interpreted  within  and  between  DRC's  in  order  to 
realize  the  full  fruit  of  the  factoring  technique.  That  is 
why  the  tables  presented  in  this  study  are  constructed  for 
ease  of  this  type  of  “within  and  between"  comparative  anal¬ 
ysis  of  factors  by  DRC  and  Region. 

Conclusions  about  the  factor  analysis  are  easily  drawn 
by  extending  beyond  the  interpretation  of  the  principal  com 
ponents  or  factors.  For  instance,  the  Race  factor  which 
has  the  slight  correlation  with  AFQT  in  certain  DRC's  can 
be  examined  nationwide.  If  one  plots  the  geographical  lo¬ 
cation  of  the  presence  of  the  Race-AFQT  factor  on  a  map, 
one  will  see  a  prevalence  of  this  factor  in  the  southern 
states,  and  in  certain  heavily  populated  areas.  This  could 
possibly  mean  that  AFQT  test  questions  are  geared  toward 
one  cultural  segment  of  the  country.  Another  factor,  the 
Army  correlation  with  No-High-School-Diploma,  factor  26  in 
Region  1,  shows  that  non-high  school  graduate  recruits  are 


principally  associated  with  the  Army.  This  may  be  an  un¬ 
palatable  interpretation  for  the  Army  to  accept.  However, 
the  idea  is  to  make  the  recruiting  manager  aware  of  the 
presence  of  a  factor  with  this  type  of  configuration  and 
thus  aid  him  in  aligning  or  reallocating  scarce  resources 
such  as  advertising  in  the  given  DRC.  For  instance,  it  may 
tell  the  recruiting  manager  that  this  particular  DRC  needs 
more  advertising  showing  high  technology  Army  weapons.  Army 
computer  technicians  at  work,  and  the  numerous  educational 
benefits  available  to  Army  recruits. 

The  overall  aim  of  this  analysis  is  to  provide  the 
manager  with  an  additional  tool  to  focus  on  the  factors 
affecting  recruiting  by  individual  DRC  and  Region.  With 
this  added  information  the  manager  should  be  able  to  make 
better  decisions  on  the  allocation  of  scarce  recruiting  re¬ 
sources;  and  the  higher  level  recruitment  planners  should 
be  better  able  to  determine  whether  to  continue  or  alter  re¬ 
cruitment  policies. 


Discriminant  Analysis 

The  primary  aim  of  discriminant  analysis  in  this  study 
was  to  discover  the  identity  and  magnitude  of  the  principal 
variables  that  could  be  used  to  distinguish  Army  accessions 
from  non-Army  accessions.  This  analysis  was  carried  out  at 
the  national  level,  the  regional  level,  and  the  DRC  level. 
Only  Regular  Forces  were  examined  at  the  DRC  level.  USAREC 
and  ARI  planners  should  be  able  to  use  the  discriminant 
analysis  procedures  developed  in  this  study  to  better 


understand  the  differences  between  Army  and  non-Army  acces¬ 
sions.  The  comparative  analysis  techniques  presented  herein 
can  aid  operations  research  analysts  at  ARI  and  USAREC  to 
differentiate  variables  that  are  critical  to  Army  recruit¬ 
ment  in  the  various  DRC's  and  Recruiting  Regions. 

Factor  analysis  and  cross-tabulation  runs  had  pointed 
to  a  difference  among  Recruiting  Regions  and  among  Components 
(National  Guard,  Reserves,  and  Regular)  for  all  analyses  in¬ 
volving  differentiation  by  the  Service  variables  (Army  and 
non-Army).  It  was  therefore  decided  to  begin  this  portion 
of  the  study  by  doing  discriminant  analysis  on  the  MU  file 
Regions  and  Components.  That  is,  the  variable  "Region"  was 
specified  with  the  SPSS  keyword  "GROUPS="  as  the  variable 
to  be  evaluated  under  discriminant  analysis.  Likewise,  the 
variable  "Component"  was  evaluated  under  an  independent 
analysis.  The  purpose  of  this  examination  was  to  test 
whether  these  two  variables  exhibited  "group  differences  as 
they  had  done  in  the  factor  and  cross-tabulation  analysis. 
Finally,  based  on  the  results  of  the  analysis  by  Region  and 
Component,  "Service"  was  evaluated  as  the  "GROUPS="  variable 
to  form  the  core  of  the  discriminant  analysis. 

The  results  of  testing  the  Regions  as  groups  and  test¬ 
ing  the  Components  as  groups  showed  that  there  was  evidence 
of  differences  both  among  Regions  and  among  Components. 

The  discriminant  analysis  test  of  Regions  revealed  that 
72.6  percent  of  the  cases  could  be  correctly  classified  by 
Region  based  on  the  discriminant  equation  using  all  the  MU 
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file  variables  except  Status  and  Component.  The  first  func¬ 
tion  had  an  eigenvalue  which  accounted  for  75.10  percent  of 
the  variance.  The  canonical  correlation  for  this  first 
function  was  very  strong  at  .8569  and  the  Wilks  Lambda  value 
was  extremely  significant  at  .1219.  The  Component  variable 
had  less  dramatic  differences  with  58.4  percent  of  the  sam¬ 
pled  cases  correctly  predicted,  a  first  function  percent  of 
variance  at  88.77,  a  canonical  correlation  of  .302,  and  a 
Wilks  Lambda  of  .8973.  The  results  of  these  two  tests  con¬ 
firm  the  suspicion  that  the  MU  variables  are  influenced  by 
regionality  and  by  component  of  service.  Thus,  as  in  the 
case  of  Factor  Analysis,  it  was  decided  to  break  the  dis¬ 
criminant  analyses  down  into  separate  tests  by  Region  and 
by  DRC  for  all  tests  concerning  the  determination  of  dif¬ 
ferences  between  the  Army  and  the  non-Army  groups. 

Sir  Maurice  Kendall  in  Multivariate  Analysis  (Ref  15: 
177)  says  that  "...  the  amalgamation  of  classes  may  give 
rise  to  misleading  conclusions" ,  as  is  demonstrated  in 
"Simpson's  paradox."  Simpson's  paradox  can  be  exhibited 
through  contingency  tables  where  two  sub-tables  are  created 
from  a  parent  table,  each  table  having  the  same  row  and 
column  categories.  The  two  sub-tables  may  be  tested  and 
yield  one  conclusion  about  the  association  of  the  tabled 
categories:  however  when  the  parent  table  is  similarly 
tested  it  yields  a  completely  opposite  conclusion  about  the 
association  between  the  categories.  For  example,  the  ta¬ 
bled  categories  may  be  Service  versus  Education  Level,  the 


two  sub-tables  may  be  Region  1  and  Region  3,  while  the 
parent  population  would  be  the  combination  of  the  two  Re¬ 
gions.  Simpson's  paradox  in  this  case  would  occur  if  the 
contingency  tables  for  Region  1  and  for  Region  3  showed 
there  was  a  difference  between  Service  and  Education  Level, 
yet  the  parent  contingency  table  composed  of  the  two  Regions 
added  together  showed  nt>  difference  existed  between  Service 
and  Education  Level.  Sir  Kendall  recommends  that  classifi¬ 
cation  type  of  analyses  be  carried  out  at  the  lowest  level 
feasible,  that  is,  that  the  analysis  be  broken  down  into 
sub-categories  to  the  lowest  level  that  still  permits  enough 
cases  per  cell  and  is  not  prohibitively  costly  in  terms  of 
computer  resources. 

Following  Sir  Kendall fs  advice,  this  study  broke  the 
analysis  down  to  the  DRC  level  in  the  case  of  the  Regular 
Forces.  However,  there  were  not  enough  differentiable  re¬ 
cords  in  the  Reserve  Forces  to  do  analysis  down  to  the  DRC 
level;  nor  were  there  enough  cases  in  the  National  Guard 
Component  to  accomplish  discrimination  at  the  Regional  level 
There  are  no  Navy  or  Marine  National  Guard  records,  and 
only  a  small  percentage  of  Air  National  Guard  on  the  file. 

It  should  also  be  noted  that  some  Regions  had  no  Air  Force 
Reserve  records  on  the  MU  file;  however,  there  were  enough 
Army  versus  non-Army  Reserve  cases  to  accommodate  analysis 
by  Region.  The  Prior  Service  and  Sex  Variables  were  used 
to  further  classify  the  final  analyses;  Regional  level  scru¬ 
tiny  of  the  MU  file  involved  discriminant  analysis  both  with 


and  without  the  prior  service  and  female  cases.  DRC  level 
study  was  done  without  prior  service  and  female  records 
since  the  main  Army  recruitment  focus  is  on  non-prior  ser¬ 
vice  male  accessions#  thus  following  Sir  Kendall's  recom¬ 
mendation  concerning  categories. 

The  USAREC  11  Intensity"  Variables 

The  USAREC  variables  on  the  MU  file  are  treated  as 
"barometers"  of  demographic  and  economic  conditions  which 
affect  each  individual  whose  records  are  on  the  MU  file. 
These  barometers  indicate  by  DRC  and  by  quarter  the  "in¬ 
tensity"  that  certain  demographic  and  economic  variables 
have  on  individuals.  Six  of  the  USAREC  variables  are  fairly 
straight-forward  predictor  variables  which  affect  the  in¬ 
dividual?  the  other  two  variables  are  usually  used  as  de¬ 
pendent  variables  in  recruitment  modelling.  The  latter  two 
variables  on  the  MU  file  are  DOD-Non-Prior-Service-Male- 
Accessions ,  and  DOD-High-School-Diploma-Graduate-Male- 
Accessions.  To  accommodate  both  the  role  of  predictor  and 
non-predictor  for  these  accession  variables,  two  sets  of 
discriminant  analyses  were  generated  down  to  the  DRC  level, 
so  that  the  user  may  choose  the  equations  and  relations 
among  variables  that  are  most  suitable  to  his  interests. 

The  following  is  a  short  summary  of  how  the  eight  "in¬ 
tensity"  variables  affect  the  individual: 

1.  Qualified  Military  Available  and  Male  High  School 
Seniors:  These  two  variables  give  an  indication  of  the 
supply  side  of  the  manpower  market  relating  to  how  much 


the  individual  must  compete  to  get  into  the  military.  The 
Army  in  1982  did  not  have  to  accept  non-high  school  diploma 
recruits.  "The  Marines  are  looking  for  only  a  few  good  men. 

2.  Black  Military  Available:  Same  as  part  1  above, 
plus  there  is  a  demographic  intensity  in  terms  of  an  im¬ 
portant  racial  group  in  today's  Army.  The  density  of  the 
black  population  in  a  given  DRC  may  well  have  an  impact  on 
whether  an  individual  (black  or  white)  will  join  the  Army 
from  a  given  DRC. 

3.  DRC  Overall  Unemployment  Rate  and  Median  Disposable 
Family  Income:  These  variables  exhibit  the  economics  of  a 
given  DRC  that  will  "drive"  individuals  toward  or  away  from 
military  service. 

4.  Army  Recruiters  As  Percent  of  DOD  Recruiters  by 
DRC:  The  "intensity"  of  Army  recruiters  as  a  ratio  of  the 
overall  number  of  DOD  recruiters  in  a  given  DRC  supplies  an 
indication  of  the  inter-service  "competition"  for  recruit¬ 
ment  in  that  DRC  and  the  individual ' s  likelihood  of  contact 
with  Army  type  recruiters  as  opposed  to  other  service  re¬ 
cruiters. 

5.  Finally,  the  two  accession  variables:  These  are 
technically  output  of  the  recruitment  process,  not  true 
predictive  variables.  However,  the  number  of  accessions 
closely  parallels  accession  quotas  which  are  based  on  pro¬ 
jections  of  a  DRC's  ability  to  provide  recruits  in  a  quarter 
Therefore,  these  accession  variables  could  be  treated  as 
accession  quotas,  giving  a  measure  of  the  "intensity"  of 


the  demand  side  of  the  recruiting  economics.  As  in  all 
supply  and  demand  systems ,  the  intensity  of  demand  in  the 
market  place  has  an  impact  on  how  many  individual  suppliers 
will  be  willing  to  come  in. 

The  Discriminant  Analysis  Tables  are  organized  in  rows 
of  standardized  coefficients  by  DRC.  This  facilitates  com¬ 
parison  within  any  given  DRC  among  the  standardized  coeffi¬ 
cients  to  show  which  variables  are  of  greatest  relative  im¬ 
portance  in  predictive  ability.  Comparison  between  DRC's 
can  be  made  by  noting  which  variables  are  important  in  one 
DRC  as  opposed  to  the  variables  that  are  important  in  an¬ 
other  DRC.  This  is  somewhat  similar  to  the  comparative 
techniques  that  were  employed  to  relate  factor  differences 
in  the  factor  analysis.  However,  variable  relations  and 
intensities  under  factor  analysis  will  not  correspond  ex¬ 
actly  with  the  importance  of  given  variables  in  discriminant 
analysis.  This  is  due  to  the  fact  that  factor  analysis 
attempts  to  determine  how  each  of  the  factors  relate  to  all 
of  the  others,  while  discriminant  analysis  tries  to  find  in 
selected  variables  the  discriminating  features  for  one 
"group"  type  variable.  Therefore,  for  the  factor  analysis 
in  this  study,  the  aim  was  to  find  relationships  among  the 
important  variables  that  may  relate  to  recruiting,  while 
discriminant  analysis  was  concerned  with  finding  the  vari¬ 
ables  and  their  relative  intensities  that  differentiated 


Army  recruits  from  non-Army  recruits. 


Appendix  C-5  through  Appendix  C-16  contain  the  numerical 
results  of  the  discriminant  analysis  performed  in  this  study. 
These  Appendices  contain  the  following  classes  of  information 
C-5  =  Discriminating  Power  of  Discriminant  Functions 

by  DRC,  No  Accession  Variables, 

C-6  =  Discriminating  Power  of  Discriminant  Functions 

by  DRC,  With  DOD  Accession  Variables, 

C-7  =  Discriminating  Power  of  Discriminant  Functions 

by  Region,  No  Accession  Variables, 

C-8  =  Discriminating  Power  of  Discriminant  Functions 

by  Region,  With  DOD  Accession  Variables, 

C-9  =  Standardized  Canonical  Discriminant  Function 

Coefficients  by  DRC,  No  Accession  Variables, 
C-10  =  Standardized  Canonical  Discriminant  Function 

Coefficients  by  DRC,  With  DOD  Accession 
Variables, 

C— 11  =  Standardized  Canonical  Discriminant  Function 

Coefficients  by  Region,  No  Accession 
Variables , 

C-12  =  Standardized  Canonical  Discriminant  Function 

Coefficients  by  Region,  With  DOD  Accession 
Variables , 

C-13  =  Percentages  of  Cases  Correctly  Classified  into 

Army  and  Non-Army  by  DRC,  No  Accession 


Variables 


C-14  =  Percentages  of  Cases  Correctly  Classified 

into  Army  and  Non-Army  by  DRC,  With  DOD 
Accession  Variables, 

C-15  =  Percentages  of  Cases  Correctly  Classified 

into  Army  and  Non-Army  by  Region,  No 
Accession  Variables,  and 

C-16  =  Percentages  of  Cases  Correctly  Classified 

into  Army  and  Non-Army  by  Region,  With  DOD 
Accession  Variables. 

Appendices  C-5  through  C-8  are  each  organized  into  five 
columns  containing:  Canonical  Correlation,  Wilks'  Lambda, 
Chi-Squared,  Degrees  of  Freedom,  and  Significance.  Appen¬ 
dices  C-6,  C-8,  C-10,  C-12,  C-14,  and  C-16  contain  only 
those  DRC 1 s  or  Regions  whose  discriminating  functions  were 
modified  with  the  addition  of  the  DOD  Accession  variables, 
implying  that  DRC's  or  Regions  which  are  not  listed  in  C-6, 
C-8,  C-10,  C-12,  C-14,  and  C-16  have  no  increase  in  predict¬ 
ability  when  the  Accession  variables  are  included.  The 
national  and  regional  level  results  in  Appendices  C-7 ,  C-8, 
C-ll ,  C-12,  C-15,  and  C-16  contain  output  that  falls  into 
two  categories: 

1.  Analysis  including  prior  service,  non-prior  ser¬ 
vice,  male,  and  female  fully  qualified  applicants,  and 

2.  Analysis  including  only  non-prior  service  male 
fully  qualified  applicants. 

The  first  of  these  categories  will  have  a  "PS,F"  (prior 

service  and  female)  indicator  in  the  column  marked 


"Component  &  Region."  The  components  are  abbreviated  "RSV" 
for  Reserve  and  "REG"  for  Regular.  The  regions  are  the 
five  standard  regional  numbers:  1,  3,  4,  5,  and  6.  Thus,  a 
typical  "Component  &  Region"  mnemonic  would  be  "RSV3  PS,F" 
meaning  Reserve,  Region  3,  with  prior  service  and  female 
applicant  cases  included.  In  all  cases  when  prior  service 
and  female  cases  are  added  to  the  analysis,  both  the  "Prior 
Service"  and  the  "Sex"  indicator  variables  come  into  the 
models  as  significant  predictors  in  the  discriminating  func¬ 
tions  for  Army  versus  non-Army.  As  was  pointed  out  earlier, 
detailed  analysis  at  the  DRC  level  focuses  only  on  the  most 
sought  after  "target"  group,  non-prior  service  males. 

The  canonical  correlation  values  in  Appendices  C-5 
through  C-8  give  one  measure  of  discriminating  power  by  show 
ing  the  relative  association  of  the  discriminant  function 
with  the  dummy  variables  Army  and  non-Army.  The  second 
measure  of  discriminating  power  is  given  in  the  value  asso¬ 
ciated  with  the  Wilks'  Lambda,  its  conversion  to  a  Chi- 
Squared  approximation,  and  the  resulting  significance  value. 
The  Wilks'  Lambda  is  an  inverse  measure  of  the  discriminat¬ 
ing  power  of  the  original  variables,  so  that  the  lower  val¬ 
ues  of  Wilks'  Lambda  imply  greater  separation  of  the  Army 
versus  non-Army  groups.  The  results  contain  a  Chi-Squared 
measurement  based  on  a  transformation  of  Wilks'  Lambda 
which  is  asymptotically  Chi-Squared  distributed.  Finally, 
the  results  are  converted  to  a  significance  value  for  hypoth 
esis  testing:  the  lower  the  value  in  the  Significance  column 


the  greater  is  the  probability  that  the  Army  versus  non-Army 
groups  are  distinguishable.  Appendix  C-5  can  be  examined  at 
alpha  equal  .01,  meaning  that  values  in  the  Significance 
column  above  .01  imply  the  rejection  of  the  hypothesis  that 
the  Army  versus  non-Army  groups  are  indistinguishable  with 
99  percent  confidence. 

Appendices  C-9  through  C-12  give  only  an  abbreviated 
form  of  the  standardized  coefficient  rounded  to  the  hun¬ 
dredths  decimal  place.  The  purpose  of  these  four  appendices 
is  to  give  the  reader  a  tool  for  comparative  discriminant 
analysis  among  DRC's  and  Regions,  not  to  give  final  exact 
equations.  Since  this  study  only  involved  a  sample  from  the 
MEPS  data  base,  it  would  be  wise  to  compute  exact  coeffi¬ 
cients  from  the  whole  data  base  rather  than  rely  on  this 
sample.  The  procedures  and  programs  as  outlined  in  this 
thesis  are  designed  to  accommodate  this  type  of  exact  anal¬ 
ysis.  Therefore,  the  results  of  the  MEPS  sample  data  are 
confined  to  comparative  analysis  within  and  among  the  DRC's 
and  Regions,  and  highlight  those  variables  most  critical  to 
these  Regions  and  DRC's. 

The  standardized  coefficients  in  Appendices  C-9  through 
C-12  show  the  relative  weights  that  the  different  variables 
have  in  predicting  Army  versus  non-Army  accessions.  All 
positive  coefficients  imply  that  as  the  value  of  the  respec¬ 
tive  variable  goes  up  so  does  the  likelihood  that  the  given 
individual  or  "target"  case  will  join  the  Army,  All  nega¬ 
tive  coefficients  say  that  the  respective  variables  point 
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toward  the  individual  having  a  tendency  to  join  another 
branch  of  service.  By  comparing  standardized  coefficient 
magnitudes  of  variables  within  a  DRC ,  one  can  come  to  under¬ 
stand  which  variables  have  a  greater  impact  in  the  DRC,  and 
by  how  much  this  impact  will  influence  an  individual  to  join 
the  Army.  To  compare  DRC's  with  each  other,  one  simply  has 
to  look  down  the  columns  of  Appendices  C-9  through  C-12  and 
see  which  DRC's  do  or  do  not  have  standardized  coefficients 
in  the  different  variable  columns.  However,  relative 
weights  of  coefficients  can  not  be  assigned  when  comparing 
between  DRC's  since  the  values  of  the  centroids  are  not  the 
same  in  each  DRC. 

Appendices  C-13  through  C-16  show  the  percentages  of 
individuals  that  were  correctly  predicted  to  go  into  the 
Army  or  not  based  on  using  Fisher ' 3  Linear  Discriminant 
Function  Coefficients  and  plugging  in  the  variable  values 
of  each  individual  on  a  case  by  case  basis.  This  final 
measure  of  the  effectiveness  of  the  classification  process 
is  divided  into  two  parts  on  each  page  of  Appendix  C-13 
through  Appendix  C-16.  The  first  part  shows  the  results  of 
using  the  classification  function  to  predict  the  groups  of 
the  analyzed  cases  which  created  the  function  in  the  first 
place.  The  second  part  shows  the  results  of  calculations 
performed  on  a  random  sample  of  the  data  different  from  the 
data  that  created  the  function. 


Results  of  the  Discriminant  Analyses 

Most  of  the  canonical  correlations  have  values  between 
.300  to  .599,  showing  that  the  derived  discriminating  func¬ 
tions  are  moderately  correlated  with  the  Army  versus  non- 
Army  groups.  DRC's  having  lower  than  "moderate"  correlation 
of  less  than  .300  are:  4D  (Denver),  4J( Oklahoma  City),  5B 

(Cincinnati),  and  6E  (San  Francisco).  When  the  DOD  Acces- 

» 

sion  Variables  are  added,  only  4D  (Denver)  increases  in 
predictability  enough  to  move  into  the  "moderate"  range. 

At  the  national  and  regional  levels,  three  analyses  have 
canonical  correlations  above  .599,  namely:  RSV3  PS,F  (Re¬ 
serve,  Region  3,  with  prior  service  and  female):  RSV3  (Re¬ 
serve,  Region  3,  no  prior  service  nor  female):  and  RSV6  PS,F 
(Reserve,  Region  6,  with  prior  service  and  female).  The 
addition  of  the  Accession  variables  at  the  national  and  re¬ 
gional  level  had  no  notable  effect  on  the  results. 

The  second  measurement  of  discriminating  power,  sum¬ 
marized  in  the  Significance  value,  shows  that  the  following 
DRC's  have  a  questionable  amount  of  discrimination  between 
the  Army  versus  non-Army  groups:  3F  (Louisville),  4D  (Den¬ 
ver),  4J  (Oklahoma  City),  5E  (Des  Moines),  6E  (San  Francisco), 
and  6H  (Phoenix).  When  the  DOD  Accession  variables  are  in¬ 
cluded,  only  4D  (Denver)  improves  in  predictability.  In 
addition,  even  at  an  alpha  value  of  .05  (95  percent  confi¬ 
dence),  there  would  still  be  questionable  discriminating 
ability  in  40  (Oklahoma  City),  6E  (San  Francisco),  and  6H 
(Phoenix).  Researchers  using  the  discriminant  function 
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coefficients  generated  for  these  particular  DRC's  should 
keep  in  mind  that  there  is  a  significantly  greater  chance 
for  error  in  the  predictive  ability  for  these  DRC's.  At  the 
national  and  regional  levels,  the  Significance  values  show  a 
confidence  level  of  better  than  99.99  percent. 

Appendices  C-9  through  C-12  show  standardized  coeffi¬ 
cient  relationships  within  and  between  the  various  DRC's  and 
Recruiting  Regions.  The  most  obvious  pattern  throughout 
these  appendices  is  the  strong  negative  relationship  that 
the  Armed  Forces  Qualification  Test  (AFQT)  scores  have  with 
people  joining  the  Army  versus  the  other  services.  Only  a 
few  DRC's  do  not  exhibit  this  relationship:  1J  (Niagara 
Falls),  3C  (Charlotte),  4A  (Albuquerque),  4C  (Dallas),  4J 
(Oklahoma  City),  5C  (Cleveland),  5E  (Jacksonville),  61  (Port¬ 
land),  and  6J  (Sacramento),  The  next  strongest  association 
of  a  large  number  of  DRC's  is  to  the  variable  "Age",  which 
shows  relatively  large  coefficient  weights  in  most  cases. 

The  Age  variable  demonstrates  that  the  older  an  individual 
is,  the  greater  is  the  likelihood  that  the  individual  will 
choose  the  Army  in  most  DRC's,  with  the  reverse  taking  place 
in  1A  (Albany),  4H  (Little  Rock),  and  5E  (Des  Moines).  At 
the  national  level  both  the  "Prior-Service"  and  the  "Sex" 
variables  are  linked  to  the  Army,  meaning  that  being  male 
and  having  had  prior  service  tends  to  help  predict  that  the 
individual  will  join  the  Army. 

An  important  variable  that  must  be  read  with  caution  is 
ArmyR  (Army  Recruiters  as  Percent  of  DGD  by  DRC).  This 
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critical  variable  can  show  by  DRC  how  much  effect  the  percent¬ 
age  of  recruiters  has  upon  an  individual  choosing  the  Army. 
However,  great  care  must  be  taken  in  analyzing  a  limited  num¬ 
ber  of  data  samples:  only  four  quarters,  that  is  four  separate 
data  points,  per  DRC  were  merged  from  the  USAREC  file  since 
the  MEPS  file  data  contained  a  limited  time  period.  There 
were  nine  DRC's  that  actually  showed  a  negative  trend  as  the 
percentage  of  Army  recruiters  increased:  that  is,  as  rela¬ 
tively  more  Army  recruiters  were  present  in  certain  DRC's, 
the  applicants  tended  to  join  the  other  services.  These 
nine  DRC's  were:  IE  (Harrisburg),  4C  (Dallas),  4D  (Denver), 

4E  (Houston),  5E  (Des  Moines),  5J  (Milwaukee),  6L  (Seattle), 
and,  when  the  Accession  variables  are  added  to  the  predictive 
model,  II  (Long  Island),  and  4F  (Jackson).  No  definite  con¬ 
clusions  should  be  drawn  from  these  results  until  further 
analysis  has  been  conducted  over  a  longer  time  period.  As 
has  been  repeatedly  stated  in  this  study,  the  driving  purpose 
of  this  thesis  effort  is  to  establish  a  set  of  programs  and 
procedures  for  the  in-depth  analysis  of  these  sort  of 
phenomena. 

The  last  set  of  appendices  on  discriminant  analysis, 

C-13  through  C-16,  give  the  results  of  the  percentage  of 
cases  which  were  correctly  classified  by  the  derived  dis¬ 
criminant  functions.  The  obvious  trend  in  these  results  is 
that  the  discriminant  variables  are  slightly  more  often  cor¬ 
rect  in  predicting  that  an  individual  will  join  the  other 
branches  of  service,  than  when  these  same  variables  are 
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used  to  predict  that  an  individual  will  join  the  Army.  That 
is,  correct  classification  of  Army  applicants  is  less  than 
the  correct,  classification  of  non-Army  applicants,  in  most 
DRC's.  However,  six  DRC's  show  consistently  better  predic¬ 
tive  ability  for  Army  applicants  over  non-Army  applicants: 

3A  (Atlanta),  3D  (Columbia/Fort  Jackson),  4A  (Albuquerque), 

5D  (Columbus),  5L  (Omaha),  and  6F  (Honolulu).  Thus,  in  these 
six  DRC's  the  percentage  of  correctly  predicted  Army  appli¬ 
cants  is  higher  than  the  percentage  of  correctly  predicted 
non-Army  applicants.  At  the  national  and  regional  levels 
these  relative  percentages  are  more  evenly  distributed,  with 
no  obvious  trend. 

In  summary,  the  results  of  Appendices  C-5  through  C-16 
lay  down  a  format  for  research  into  the  critical  variables 
affecting  Army  versus  non-Army  recruitment.  The  analysis 
that  surrounds  these  appendices  provides  a  set  of  procedures 
to  be  used  in  examining  the  key  relationships  among  the  vari¬ 
ables  and  their  weighted  effect  upon  the  dependent  group 
variables  Army  and  non-Army.  The  analysis  carried  out  down 
to  the  DRC  level  gives  the  researcher  a  clear  picture  of  the 
effects  of  the  specific  set  of  variables  that  affect  recruit¬ 
ment  within  a  given  DRC.  At  the  same  time,  the  analyst  is 
able  to  compare  coefficient  weights  to  determine  relative 
strengths  among  the  critical  variables.  At  the  global  level, 
the  structure  of  the  appendices  allows  the  researcher  to 
analyze  the  differences  between  DRC's,  to  note  which  vari¬ 
ables  deserve  further  investigation  over  different  periods 


of  time,  and  which  DRC's  exhibit  noteworthy  countermovements 
to  overall  trends  in  the  data  base.  The  techniques  and  for¬ 
mats  developed  in  this  study  can  easily  be  expanded  through 
the  incorporation  of  more  variables  from  Census  Bureau  sta¬ 
tistics  or  other  demographic  and  economic  information  avail¬ 
able  to  USAREC  and  ARI.  These  results  can  then  be  further 
refined  using  longer  time  periods,  and  more  data  samples  for 
the  National  Guard  and  Reserve.  The  ultimate  aim  would  be  to 
create  very  precise  predictive  models  at  the  DRC  level  that 
could  be  used  in  ARI 1 s  continuing  research  concerning  re¬ 
cruitment  indicators,  and  in  USAREC' s  efforts  to  optimally 
allocate  scarce  resources. 


VII.  TIME  SERIES  ANALYSIS 


Introduction 

This  chapter  presents  the  time  series  analysis,  utiliz¬ 
ing  the  techniques  of  Box  ar 1  Jenkins,  as  performed  on  the 
accessions  data  shown  in  Appendix  A-2.  The  modelling  was 
done  in  two  phases.  First,  Male  High  School  Graduate  Mental 
Category  I-IIIA  accessions  { ACC13AM )  was  modelled  as  a  simple 
time  series  and,  second,  unemployment  rates,  the  producer 
price  index  and  the  prime  rate  were  used  as  input  variables 
with  ACC13AM  as  the  output  variable  in  three  leading  indi¬ 
cator  models. 

General  Information 

It  was  initially  hoped  to  use  the  data  directly  from 
the  USAREC  file  to  model  ACC13AM.  However  there  were  only 
22  data  points  (quarterly  from  mid  1976)  and,  in  order  to 
perform  meaningful  analysis  using  Box  and  Jenkins  techniques, 
at  least  50  data  points  are  needed  (Ref  1:18).  Thus  began 
the  "Great  Data  Chase." 

In  earlier  readings  on  the  subject,  to  validate  the 
hypothesis  that  the  Army  advertising  dollar  was  well  spent, 
N.W.  Ayer  had  used  Male  High  School  Graduate  Mental  Category 
I-IIIA  monthly  accession  data  from  April  1976  to  March  1980. 
Accession  data  was  requested  from  the  United  States  Army 
Recruiting  Command  (USAREC)  from  January  1980  to  the  present. 
The  overlapping  data  was  requested  in  order  to  insure  a 
"good"  set  of  data  when  paired  with  the  N.W.  Ayer  data. 
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Unfortunately ,  the  new  data  from  USAREC  did  not  match  the 
N.W.  Ayer  data.  This  anomoly  was  explained  by  a  USAREC 
official  as  a  problem  of  renorming  the  accession  rates. 

In  the  late  seventies  and  early  eighties  tests  that 
determined  a  potential  recruit's  mental  category  were  com¬ 
promised  and  personnel  were  categorized  in  mental  categories 
higher  than  they  should  have  been.  Recently,  in  an  attempt 
to  rectify  the  bookkeeping,  records  have  beeh  checked  and 
the  accession  rates  adjusted  to  reflect  a  more  accurate  re¬ 
presentation  of  the  quality  of  recruits  entering  the  services 
These  renormed  rates  were  not  available  when  N.W.  Ayer  con¬ 
ducted  their  study.  USAREC  personnel  suggested  contacting 
the  "keepers  of  the  data,”  the  Defense  Data  Manpower  Center 
(DMDC)  in  Monterey,  California.  From  the  Army  Liason  Officer 
in  DMDC,  monthly  renormed  Male  High  School  Graduate  Mental 
Category  I -I II A  accession  data  (from  October  1975  to  July 
1982)  was  received. 

Once  the  data  was  received,  the  81  points  were  analyzed 
in  two  ways.  First,  the  points  were  analyzed  as  a  simple 
time  series.  This  was  done  in  order  to  compare  predictions 
with  the  second  phase  of  the  time  series  analysis.  The 
second  phase  utilized  Box  and  Jenkins  leading  indicator 
techniques  with  the  unemployment  rate,  producer  price  index 
and  prime  rate  as  input  or  independent  variables  and  the 
accession  data  ( ACC13AM)  as  the  output  or  dependent  variable. 

Modelling  (Phase  I)  -  ACC13AM  as  a  Time  Series 

Appropriate  computer  output  for  this  phase  of  modelling 


is  in  Appendix  D-l  through  D-6. 

Accession  rates  were  modelled  as  each  value  being  a 
function  of  previous  values,  in  order  to  have  a  basis  for 
comparison  with  the  leading  indicator  models  to  be  discussed 
later.  Two  methods  were  used  to  analyze  the  data;  FORTRAN 
programs  that  were  written  while  studying  Box  and  Jenkins 
techniques  (referred  to  as  TSP)  aided  in  model  identification 
and  residual  analysis,  and  a  standard  time  series  analysis 
package  in  BMDP  was  utilized  to  fine-tune  models  and  fore¬ 
cast  future  values.  The  TSP  programs  employed  various  sub¬ 
routines  from  the  IMSL  library. 

Initially  the  simple  and  partial  autocorrelations  were 
examined  for  significant  lags  or  a  recognizable  pattern  such 
CL'  as  a  damped  exponential  or  a  sine  wave.  The  TSP  pxi^grains 

generated  plots  of  both  the  simple  and  partial  autocorrela¬ 
tions.  These  are  shown  in  Figures  la  and  lb.  The  two  stan¬ 
dard  error  band  was  estimated  as: 

2  SE  Band  =  ♦  — — 

”<\j  n 

where  n  is  the  number  of  observations  in  the  time  series. 

The  simple  autocorrelation  plot  might  be  construed  as 
either  a  linearly  decreasing  function  indicating  nonstation- 
arity  or  a  very  long  sine  wave  where  only  a  small  portion  of 
the  wave  is  plotted.  The  partial  autocorrelation  plot  has 
a  significant  autocorrelation  at  lags  1  and  15.  This  could 
be  an  autoregressive  process  of  order  1,  assuming  that  the 
significance  at  lag  15  can  be  attributed  to  spurious 
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Figure  la 

Autocorrelation  Plot  of  ACC13AM 
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Partial  Autocorrelation  Plot  of  ACC 13AM 
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movement.  Before  any  model  was  chosen  additional  techniques 
were  used  in  model  identification. 

In  order  to  analyze  the  periodogram  for  significance, 
the  average  of  the  intensity  values  is  calculated  to  compare 
all  intensity  values  for  significance.  This  is  done  by  us¬ 
ing  Table  7.1.2  in  Fuller  (Ref  8:284)  to  find  the  appropriate 
multiple  (based  on  number  of  ordinates  and  the  desired  con¬ 
fidence  level).  The  average  intensity  times  the  multiple 
(hereafter  referred  to  as  periodogram  multiple)  is  computed 
and  compared  with  all  intenisty  values.  The  average  was 
approximately  2,421,741  and  the  periodogram  multiple  for  a 
40-value  periodogram  at  the  95  percent  confidence  level  is 
6.295.  Therefore  any  intensity  value  over  15,244,860  is 
significant.  The  only  intensity  value  greater  than  the 
computed  statistic  (15,244,860)  is  at  a  period  of  81  (Ap¬ 
pendix  D— 1 ) .  This  does  not  indicate  a  period  of  81  but 
hints  that  the  full  period  of  the  data  has  not  yet  been  com¬ 
pleted  in  the  81  data  points. 

Another  tool  used  in  model  identification  was  developed 
by  H.L.  Gray,  G.  Kelly,  and  D.  Mclntire  (Ref  9).  This  tech¬ 
nique  utilizes  the  R  and  S  array  method  to  aid  in  ARIMA 
model  identification.  By  recognizing  certain  patterns  in 
these  arrays  the  p,  d,  and  q  (auto  regressive,  number  of 
differences,  and  moving  average)  operators  can  be  easily 
identified  in  the  ARIMA  model. 

Using  this  methodology,  the  R  and  S  arrays  (at  low  and 
frigh  frequency)  were  examined  (Appendix  D-2).  The  low 
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frequency  values  of  the  S  array  in  Column  1  are  all  very 
close  to  -2  thus  also  indicating  a  possible  nonstationarity 
in  the  data  and  the  need  to  non-seasonally  difference  the 
data  once  to  remove  this  nonstationarity.  This  was  accom¬ 
plished  using  the  TSP  programs  as  in  earlier  analysis. 

For  the  differenced  data,  there  were  no  significant 
simple  autocorrelations,  although  there  were  several  border¬ 
line  cases  (Appendix  D-3).  The  only  significant  partial 
autocorrelation  was  at  lag  two,  indicating  that  the  differ¬ 
enced  data  might  be  modelled  by  an  AR(2)  process.  Analysis 
of  the  periodogram  of  the  differenced  data,  using  the  same 
periodogram  multiple  as  before,  indicated  that  any  value 
greater  than  5.7  million  is  significant  at  the  95  percent 
confidence  level.  There  were  no  such  values,  thus  there 
are  no  significant  periods  to  be  modelled  out  of  the  differ¬ 
enced  data  (Appendix  0-3). 

The  R  and  S  arrays  for  the  differenced  data  indicated 
several  possible  models:  ARMA  (2,6),  ARMA  (4,2),  and  AR(6), 
all  of  which  seemed  to  be  an  overfit  of  the  simple  AR(2) 
model.  Appendix  D-4  is  a  listing  of  the  autocorrelation 
plots,  periodogram  values,  and  the  R  and  S  arrays  for  the 
differenced  data.  At  this  point,  the  BMDP  time  series  pack¬ 
age  was  used  to  help  identify  the  model  and  perform  various 
adequacy  checks.  Since  our  initial  thought  was  that  the 
process,  once  differenced,  was  autoregressive  of  order  two, 
the  data  was  differenced  and  modelled  as  AR(2). 


Coefficients  were  estimated  for  the  autoregressive 
parameters  by  the  conditional  least  squares  method  and  by 
backcasting.  The  coefficients  for  both  the  AR(1)  and  AR(2) 
terms  were  statistically  significant.  The  residuals  were 
estimated  using  these  coefficients  and  then  the  residual 
series  was  examined  for  adequacy  of  the  fit  of  the  proposed 
model.  If  the  model  is  adequate  the  residual  series  would 
not  exhibit  any  significant  simple  or  partial  autocorrela¬ 
tions.  The  residual  series  from  the  ARI(2(1)  model  had  no 
significant  simple  autocorrelations  (but  large  at  lag  14) 
and  one  significant  partial  autocorrelation  at  lag  14.  For 
a  listing  of  BMDP  output  for  this  model  see  Appendix  D-5. 

The  foregoing  results  can  be  taken  in  two  ways.  When 
working  at  a  95  percent  level  of  confidence  in  100  lags  one 
might  expect  5  to  be  significant  but  still  have  an  adequate 
model,  or  as  in  this  case  in  20  lags  one  might  expect  1 
significant  autocorrelation.  Taking  this  point  of  view  one 
would  accept  this  model  as  being  adequate  and  the  model  for 
differenced  data  would  be: 

Zfc  =  -.2763Zt-1  -  .3269Zt_2  +  afc 

where  a^  represents  a  white  noise  residual  and  Z^  are  the 
various  values  of  the  differenced  data  at  time  t.  The  mean 
of  this  process  was  determined  to  not  be  statistically  sig¬ 
nificant  so  it  is  not  in  the  equation.  Since  Zt  are  differ* 
enced  data,  the  equation  for  undifferenced  data  is 

(1-B)  (1  +  . 2763B  +  .3269B2)  X.  =  a„ 


Xt  =  .7237Xt-1  -  ,0506Xt_2  +  .3269Xt_3  +  afc  (7-1) 

However,  if  the  significance  of  the  partial  autocorrelation 
at  lag  14  is  deemed  to  be  unacceptable,  an  autoregressive 
term  of  order  14  would  be  added  to  the  model  and  the  residual 
series  checked  for  any  remaining  significant  simple  or  par¬ 
tial  autocorrelations.  This  was  done  using  BMDP.  See  Ap¬ 
pendix  D-6  for  a  listing  of  output  for  this  model.  Once 
again  all  coefficients  were  statistically  significant  and 
none  of  the  partial  or  simple  autocorrelations  were  signifi¬ 
cant  for  the  residual  series.  Thus  a  multiplicative  seasonal 
ARIMA  (2,1,0)*(1,0,0)^4  is  also  tentatively  considered  as  an 
adequate  model  with  the  following  equation: 

Zt  *  -.lSSSZj.^  -  ,2429Zt_2  +  ,2670Zt_14  +  afc 

where  Zt  once  again  represents  differenced  data  and  at  the 
white  noise  residual  series.  Thus  the  equation  for  undif¬ 
ferenced  data  is: 

( 1— B )  (1  +  . 1835B  +  .2429B2  -  .2670B14)  Xfc  -  at 
or 

Xfc  *  .8165Xt-1  -  .0594Xt_2  +  ,2429Xt_3  +  .2670Xt_14 
-  ,2670Xt_15  +  at  (7-2) 

Before  deciding  on  a  final  model  for  accessions  there 
are  two  other  checks  to  apply  to  the  model.  The  Portmanteau 
lack  of  fit  test  (Ref  1:290-293)  looks  at  the  autocorrela¬ 
tions  of  the  residuals  taken  as  a  whole  rather  than  consider¬ 
ing  then  individually. 

K  7  A 

=  n2  r l  (§) 

k»l 


Computing  Q 


where  n  is  the  number  of  data  points  (minus  the  number  of 
times  differenced),  K  is  the  number  of  autocorrelations  con¬ 
sidered  and  r£  (a)  is  the  estimated  value  of  the  residual 

2 

series  autocorrelation  at  lag  k.  Q  has  an  approximate  X 
distribution  with  (K-p-q)  degrees  of  freedom,  with  p  and  q 
representing  the  number  of  components  in  the  model  average 
components  of  the  model.  If  the  model  is  inappropriate, 
the  average  value  of  Q  will  be  inflated  or  larger  than  the 
X2  statistic  with  (K-p-q)  degrees  of  freedom  at  the  desired 
confidence  level. 

For  the  two  models  considered,  the  Portmanteau  results 
at  a  95  percent  confidence  level  were: 

ARI  (2,1)  ARIMA  (2,1,0)*(1,0,0)14 

Q  -  24.298  Q  *  25.723 

X(38)  *  53-36  X(37)  ®  52«16 

therefore  both  models  given  by  equations  (7-1)  and  (7-2) 
are  adequate. 

A  final  model  check  is  a  graph  of  the  cumulative  period- 
ogram  (Ref  1:294-297).  This  is  done  to  insure  that  the  model 
has  accounted  for  all  periodic  components.  The  plots  of  the 
cumulative  periodograms  for  both  models  are  at  Figure  2  and 
3  along  with  the  autocorrelations  of  the  residual  series 
used  in  the  calculation  of  the  Portmanteau  Q  value.  Both 
plots  were  within  the  acceptable  Kolomogorov-Smirnov  limits. 
Additionally,  both  plots  were  almost  identical  in  nature. 

This  is  expected  because  the  ARIMA  (2,1,0)*(1,0,0)^4  is  an 
over fit  of  the  ARIMA  (2,1,0)  model. 


Since  both  models  were  adequate  fits  for  the  data,  the 
BMDP  Time  Series  package  was  used  to  predict  future  results 
for  both  models  as  well  as  for  predicting  values  at  various 
points  in  time  throughout  the  data  (Appendix  E-9).  In  anal¬ 
yzing  these  forecasts,  one  might  conclude  the  ARIMA  (2,1,0)* 
(1,0,0)^4  model  is  better  because  it  has  a  smaller  standard 
error  band.  However,  in  much  the  same  manner  one  would  ex- 
pect  a  higher  r  in  regression  analysis  when  an  additional 
variable  is  entered,  a  tighter  standard  error  band  is  ex¬ 
pected  when  a  time  series  model  is  overfit.  In  most  in¬ 
stances,  the  first  one  or  two  predictions  are  quite  accurate 
but,  as  one  proceeds  further  away  from  the  start  point,  the 
predictions  as  compared  to  actual  values  showed  marked  dif¬ 
ferences.  This  helps  to  point  out  the  reason  for  using  the 
leading  indicator  technique.  When  one  has  a  model  that  uses 
predicted  values  to  predict  values  further  into  the  future, 
the  standard  error  band  will  get  enormously  wide.  What  is 
hoped  for  in  the  leading  indicator  portion  is  to  find  an  in¬ 
dependent  variable  that  leads  accession  rates.  This  allows 
for  use  of  actual  current  values  of  the  independent  vari¬ 
able  to  predict  future  values  of  the  dependent  variable. 


Leading  Indicator  Modelling  (Phase  II) 

Unemployment  Rates  as  the  Input  Variable.  Many  studies 
have  been  done  attempting  to  relate  fluctuations  in  the 
economy,  particularly  the  unemployment  rate,  to  accessions. 
Recently  Dale  and  Gilroy  (Ref  5)  U3ed  a  time  series  analysis 
technique  to  show  a  correlation  between  unemployment  rate 
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Figure  2 

Cumulative  Periodogram  for  ARI(2,1)  ACC13AM  Model 
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Figure  3 

Cumulative  Periodogram  for  ARIMA( 2,1,0 )*(1 ,0,0)^ 


and  accessions.  They  used  regression  analysis  with  various 
lags  on  unemployment  rate  as  independent  variables.  Their 
dependent  variable  value,  Male  HSG  Mental  Category  I-IIIA 
Accessions,  are  the  same  DMDC  data  used  in  this  analysis. 

This  phase  of  the  analysis  will  extend  the  Dale  and  Gilroy 
study  by  using  Box-Jenkins  leading  indicator  techniques. 

Hopefully  the  forecasts  on  accessions  rates  can  be 
considerably  improved  by  associating  these  forecasts  with 
another  variable.  This  will  be  the  case  if  changes  in  ac¬ 
cession  rates  are  anticipated  by  changes  in  unemployment 
rates.  The  first  step  in  the  leading  indicator  technique 
is  to  model  the  input  variable,  unemployment  rates.  TSP  was 
used  in  model  identification.  The  simple  autocorrelation 
plot  of  the  raw  data  appeared  (Appendix  E-l )  to  be  either 
the  middle  of  a  long  sine  wave  or  a  linearly  decreasing 
function.  This  pattern  is  recognized  as  indicating  nonsta- 
tionarity  in  the  data  which  necessitates  differencing.  The 
periodogram  implies  the  frequency  of  the  data  might  be  longer 
than  the  81  data  points  utilized.  The  S-Array  at  low  fre¬ 
quency  tended  to  the  value  two  in  column  one,  thus  indicat¬ 
ing  that  one  root  of  the  polynomial  function  that  explains 
the  data  lies  on  the  unit  circle,  which  also  suggests  a  need 
to  difference  the  data  before  proceeding  with  the  modelling. 

After  differencing,  the  unemployment  rate  data  was 
analyzed  for  significant  autocorrelations.  There  were  no 
simple  or  partial  autocorrelations  significant  thus  indicat¬ 
ing  that  all  that  is  needed  is  to  difference  unemployment 


rate  data  and  the  pattern  of  the  remaining  residuals  is  white 
noise*  However,  closer  examination  shows  that  there  are 
several  large,  though  not  significant,  simple  and  partied, 
autocorrelations  (Appendix  E-2).  Knowing  that  the  2  stan¬ 
dard  error  band  is  an  approximation  (Ref  1:32-34),  it  was  de¬ 
cided  to  further  model  the  data.  The  average  intensity  from 
the  periodogram  was  *129*  This  multiplied  by  the  periodo- 
gram  multiple  of  6.295,  for  40  periodogram  values  at  o C  = 

.05  gives  a  value  of  approximately  .81.  There  are  no  inten¬ 
sities  greater  than  this  value  therefore  seasonal  differenc¬ 
ing  is  not  necessary.  The  S-arrays  of  the  differenced  data 
initially  identified  an  AR(1),  AR(2)  and  ARMA  (4,2)  as 
possible  models.  Appendix  E-2  is  a  listing  of  the  appro¬ 
priate  TSP  computer  output.  BMDP  was  used  to  estimate  the 
coefficient  of  the  autoregressive  parameter,  calculate  re¬ 
siduals  and  check  the  residual  series  for  significant  auto¬ 
correlations.  As  expected  there  were  no  significant  simple 
or  partial  autocorrelations  and  the  estimated  AR(l)  coef¬ 
ficient  of  .1332  was  not  statistically  significant  (T  values 
all  less  than  2.0).  Next  the  ARMA  (4,2)  model  was  investi¬ 
gated.  Remembering  that  any  of  these  models  are  an  overfit 
to  an  ARIMA  (0,1,0),  one  does  not  expect  to  find  any  signi¬ 
ficant  autocorrelations  and  there  were  none.  Once  again 
none  of  the  estimated  coefficients  were  statistically  signi¬ 
ficant.  These  mode  s  were  herefore  rejected.  Next  the 
AR(2)  was  investigate:;.  There  were  no  significant  auto¬ 
correlations.  The  Portmanteau  Q  value  was  well  within  the 


acceptable  region  and  the  plot  of  the  cumulative  periodo- 
gram  values  shows  no  model  deficiencies.  However,  the 
estimated  coefficients  were  not  statistically  significant 
and  this  model  was  also  rejected. 

Since  differencing  the  data  yielded  no  significant 
autocorrelations,  it  was  attempted  to  model  the  unemployment 
rate  data  without  differencing.  First,  an  AR(1)  model  was 
applied  to  the  data.  The  AR(1)  coefficient  was  significant 
and  there  were  no  remaining  significant  autocorrelations  in 
the  residual  series.  The  Portmanteau  Q  statistic  was  cal¬ 
culated  (22.8)  and  compared  to  a  chi-squared  table  value 
with  39  degrees  of  freedom  and  this  model  appeared  to  be 
adequate.  Second,  an  AR( 2 )  model  was  applied  to  the  data. 
Both  AR  coefficients  were  significant  and  there  were  no 

significant  autocorrelations  in  the  residual  sequence.  The 

2 

Portmanteau  Q  was  calculated  (21.5)  and  compared  to  a  X 
table  value  with  38  degrees  of  freedom  and  this  model  ap¬ 
peared  to  be  adequate.  The  cumulative  periodogram  values 
for  these  two  models  were  plotted  with  no  values  outside 
the  K-S  limits,  and  a  visual  comparison  with  the  AR(1)  plot 
clearly  indicates  a  better  fit  on  the  AR(2)  model. 

A  comparison  of  the  adequate  models  is  shown  in  Table 
3  below.  See  also  Appendix  E-3  for  pertinent  output  in 
model  identification  for  unemployment  rates. 


Coefficients 
(T-values ) 


AR(1)  .9940 
(1273) 


AR( 1 )  1.1322 

(1071.85) 


AR(2)  -.1363 
(29.18) 


ACF/PACP 


No  significant  NSL 
Lags  (NSL) 


Q  22.8  21.5 

Significance  (*)  .025  .015 


Cumulative-Periodogram  Within  K-S 
Plot  Limits 


Very  tight  about 
the  center  line 


Always  above 
center  line 


The  ARIMA  (2,0,0)  model  (Model  B)  was  selected  as  the 
best  model  for  unemployment  rates  because  it  has  a  smaller 
Q  value,  its  cumulative  periodogram  plot  is  visually  better 
than  Model  A  and  it  has  a  higher  significance  level.  There 
fore  the  resulting  AR(2)  model  for  unemployment  rates  iss 
(1  -  1.1322B  +  .1363B2)  =  afc 

or 

Xt  =  1.1322Xfc_1  -  .1363Xt_2  +  afc  (7-3) 

where  Xt  is  unemployment  rate  at  time  t  and  afc  is  the  as¬ 
sociated  white  noise  residual. 

The  next  phase  in  the  leading  indicator  modelling  pro¬ 


cess  is  to  filter  the  output  series  (accessions)  by  the 
ARIMA  model  for  tho  input  (unemployment  rates)  series. 


In  other  words ,  the  same  exact  formula  as  (7-3)  above  is 
implied  to  accessions: 

Yt  -  1.1322Yt_1  -  ,1363Yt_2  +  bfc  (7-4) 

The  residuals  (bfc)  in  this  model  are  not  likely  to  be  white 
noise  as  they  (afc)  cure  in  the  input  model.  The  residuals 
from  the  filtered  output  are  then  checked  for  cross  correla¬ 
tions  with  the  white  noise  residuals  from  the  input  model. 
The  cross  correlation  plot  aids  in  the  identification  of 
the  functional  form  of  the  transfer  function  components. 
Additionally  the  simple  and  partial  autocorrelation  func¬ 
tion  of  the  filtered  output  residuals  are  examined  for  any 
additional  significant  ARIMA  components  to  be  included  in 
the  transfer  function. 

These  steps  were  applied  to  accession  rates  using  BMDP. 
The  resulting  cross  correlation  and  autocorrelation  plots 
are  at  Figures  4,  5a  and  5b.  The  cross  correlation  plot 
showed  no  significant  positive  lags.  However,  if  the  2 
S.E.  band  were  tightened  slightly  there  might  be  significant 
cross  correlations  at  lags  1,  3,  7  and  8.  Thus,  from  the 
sample  cross  correlations  between  the  prewhitened  input 
(unemployment  rate)  series  and  the  filtered  accession  rates 
the  preliminary  indications  were  a  transfer  function  of  the 
following  form: 

Yt  ®  <u1b1  +  U2fi3  +  U3B?  +  U4b8)  Xt  +  Et 
where  Y^  and  are  the  respective  output  and  input  series 
with  their  means  removed.  is  a  white  noise  residual. 
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Figure  4 

Cross  Correlation  Plot  of  Filtered  Unemployment  Rate  Model 
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Figure  5a  r 

Simple  Autocorrelation  Plot  of  Filtered 
Unemployment  Rate  Model 
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Figure  5b 

Partial  Autocorrelation  Plot  of  Filtered 
Unemployment  Rate  Model 


A  full  explanation  of  the  derivation  of  a  transfer  func¬ 
tion  model  is  found  in  Reference  1.  Additionally,  the 
simple  autocorrelation  function  is  significant  at  lags  1 
and  14  suggesting  a  possible  MA(1)  component  in  the  model 
while  the  partial  autocorrelation  function  is  significant 
at  lags  1,  2  and  14  suggesting  a  possible  AR(2)  component 
in  the  model.  Attempts  to  model  out  these  significant  lags 
by  using  the  moving  average  of  order  1  failed. 

Fine-tuning  the  transfer  function  model  by  adding  an 
AR(2)  model  along  with  the  U  coefficients  at  lags  1,  3,  7 
and  8  had  a  remaining  significant  partial  autocorrelation 
at  lag  6.  Thus,  an  autoregressive  term  of  order  6  was 
added  to  the  model  and  the  model  was  reestimated,  residuals 
calculated  and  autocorrelation  plots  of  the  residual  series 
examined.  The  coefficient  for  the  second  order  autoregres¬ 
sive  term  and  the  U  coefficient  at  lag  3  were  not  signifi¬ 
cant  so  that  they  were  dropped  from  the  model.  The  remain¬ 
ing  components,  AR  orders  (1,6)  and  UP  orders  (1,7,8)  were 
estimated  and  the  resulting  BMDP  output  is  at  Appendix  E-4. 

A  summary  of  the  model  follows: 

Coefficient 

Variable  Order  Estimate  T-Ratio 


ACC13AM 


ACC13AM 


.ei 


Both  the  simple  and  partial  autocorrelation  plots  are 


within  the  two  standard  error  bands.  These  results  lead 
one  to  the  conclusion  that  this  is  a  tentative  transfer  func 
tion  model.  This  is  almost  exactly  what  we  were  looking  for 
a  variable  that  leads  accession  rates.  This  model  says  the 
important  unemployment  rates  to  predict  current  accession 
rates  are  one,  seven  and  eight  months  back  from  the  time 
you  are  trying  to  predict.  It  was  hoped  that  by  dropping 
the  first  lag  and  only  estimating  the  U  order  coefficient  at 
lags  7  and  8  the  model  would  still  be  adequate.  This  was 
done,  however  the  7th  U  order  coefficient  was  no  longer 
statistically  significant,  thus,  it  was  dropped.  The  trans¬ 
fer  function  model  was  reevaluated  only  with  AR  (1,6)  and 
UP(  8 )  components.  In  the  resulting  model  the  lag  at  8  was 
determined  to  no  longer  be  statistically  significant.  There¬ 
fore  we  returned  to  the  AR  (1,6),  UP  (1,7,8)  as  a  tentative 
transfer  function  model.  The  resulting  transfer  function 
is: 


Yt  +  *u1b1  +  U2B?  +  U3bS )  xt  + 


I"**!3.  2B 


6  at 


or 


Yfc  ■  702.9Xt_1  -  441.6Xt-?  +  378.9Xt_Q  + 


i-. 5288b1-. 4605b6 


Diagnostic  Checking 

Before  accepting  the  above  transfer  function  as  an 
adequate  representation  of  the  system  modelled,  autocorrela¬ 
tion  and  cross  correlation  checks  should  be  applied  (Ref  1: 
393-397 ) .  Both  the  Q  and  S  criterion  are  derived  in  the 


same  manner  as  the  Portmanteau  Q  statistic  discussed  earlier 
The  first  24  lags  of  the  residual  autocorrelation  are  used 
in  the  computation  of  the  Q  criterion: 


24 


Q  =  mS  <a(K) 
K=1  aa 


where  m  is  the  actual  number  of  residuals  available  for  the 
computation: 

Q  =  71  (.2752)  S  12.44 

2 

Comparison  of  Q  =  12.44  with  the  X  tabled  value  for  K-p-q 
degrees  of  freedom  (where  p  and  q  are  from  the  noise  model) 
shows  no  grounds  for  questioning  the  model's  adequacy  at  an 
o*  level  of  .005. 

Using  the  first  36  values  of  the  cross  correlations  be¬ 
tween  the  output  residuals  and  the  prewhitened  input  resi¬ 
duals  yields  the  S  criterion: 


35 


S  -  m  S  r*  (K)  *  71(. 157994)  S  11.218 
K*=0  aa 


Comparison  of  S  with  the  tabled  value  with  (K  +  1)  - 
(r  +  s  +  1)  degrees  of  freedom  at  the  .005  level  again  shows 
no  grounds  for  questioning  the  model's  adequacy. 


Forecasting  Under  Transfer  Function  Models 

In  order  to  forecast  accession  rates,  it  is  necessary 
to  provide  forecasts  of  the  unemployment  rates.  This  is 
done  using  BMDP.  Once  the  unemployment  rates  are  forecast¬ 
ed,  the  transfer  function  is  used  to  predict  accession  rates 
A  BMDP  listing  of  this  output  is  at  Appendix  E-5. 


A  comparison  of  these  results  with  the  results  of 
modelling  accession  rates  alone  (Phase  I)  shows  that  the 
forecasts  for  short  lead  times  under  the  leading  indicator 
model  have  a  much  narrower  error  band  thus  increasing  the 
accuracy  of  the  forecast.  These  results  are  shown  in 
Table  V. 

Table  V 

Comparison  of  LI  Model  and  Simple 
Time  Series  Forecasts 


Period 

Forecasts 

Forecasts 

Ahead 

Accessions 

Leading 

X 

Alona  ARK 2.1) 

SE 

Indicator 

SE— LI 

1 

4900 

599 

4777 

519 

2 

4745 

740 

5045 

626 

3 

4852 

792 

5969 

688 

4 

4873 

878 

4802 

735 

5 

4837 

966 

4494 

775 

6 

4849 

1032 

4731 

811 

Model  Enhancing 

If  the  reader  refers  to  Appendix  E-4  the  simple  and 
partial  autocorrelations  at  lag  1  are  both  large  but  not 
significant.  An  effort  was  made  to  reduce  the  size  of  this 
autocorrelation  hopefully  reducing  even  further  the  stan¬ 
dard  error  band  on  future  predictions.  Since  there  were 
already  first  and  second  order  autoregressive  components  in 
the  transfer  function,  a  moving  average  component  of  order 
one  (MA(  1 ) )  was  added  to  the  model.  Coefficients  were 


estimated,  residuals  calculated  and  the  residual  series 
checked  for  significant  simple  or  partial  autocorrelations. 
Appendix  E-6  is  a  listing  of  BMDP  output  for  this  model. 


The  model  is  summarized  as  follows: 


Coefficient 


Variable 

Order 

Estimate 

T-Ratio 

ACC 

MA 

1 

.3169 

1.91 

ACC 

AR 

1 

.6832 

6.64 

ACC 

AR 

6 

.3114 

3.07 

Unemp 

UP 

1 

701.92 

5.26 

Unemp 

UP 

7 

-493.81 

2.04 

Unemp 

UP 

8 

426.99 

1.93 

All  of  the 

coefficients  were 

significant 

at  the  .05 

.  and  there 

were  no 

significant  simple  or 

partial  auto 

correlations. 

The  resulting  transfer  function  model  is: 

Yfc  =  701.92Xt_1  -  493.81Xt_?  +  426.99Xfc-8  + 

( 1**»  3169B )  __  a 

.  1-. 6832B-. 3114B6  t 

Since  this  is  an  overfit  of  a  previously  calculated 
model  there  is  no  doubt  about  its  adequacy.  A  comparison 
between  the  forecasts  of  the  two  models  is  shown  in 
Table  VI. 

Table  VI 

Comparison  of  Unemployment  Rates 
Leading  Indicator  Models 


Period 

Ahead 

X 

Model  1 
AR(1,6) 
UP(1,7 ,8) 
Forecast 

SE 

1 

4777 

519 

2 

5045 

626 

3 

4969 

688 

4 

4802 

735 

5 

4494 

775 

6 

4731 

811 

Model  2 
AR( 1 , 6 ) 

MA(  1 ) 

UP( 1 ,7 , 8 ) 

Forecast  SE _ Actual 


4612 

523.6 

4912 

596.5 

4822 

640.1 

4696 

673.4 

4419 

702.2 

4708 

728.1 

103 


It  should  be  noted  that  in  all  cases  Model  2's  pre¬ 
dictions  are  less  than  Model  l's  and  the  desired  effect  of 
tightening  the  SE  band  happened  at  all  predictions  except 
the  first. 

However,  when  considering  a  2  SE  band  about  Model  l's 
first  prediction  it  can  be  said  with  approximately  95  per¬ 
cent  confidence  that  the  true  value  of  the  accessions  for 
the  next  month  will  lie  between  3739  and  SSl^. 

The  conclusion  reached  from  this  point  is  that,  re¬ 
gardless  of  the  fact  that  both  transfer  function  models  are 
adequate  representations  of  the  data  and  they  have  tight¬ 
ened  the  2  SE  band  when  compared  to  the  model  from  Phase  I, 
the  models  have  too  wide  a  standard  error  band,  even  on  the 
first  prediction,  to  be  models  that  Army  planners  can  use 
to  predict  future  accession  rates. 

Producer  Price  Index  as  Input  Variable.  Since  unemp¬ 
loyment  rates  failed  to  produce  a  usable  leading  indicator 
model,  a  different  type  of  variable  was  chosen  to  model  with 
accessions.  The  producer  price  index  (PPI)  much  like  the 
consumer  price  index  is  a  measure  of  inflation.  PPI  was 
chosen  as  the  independent  variable  in  the  leading  indicator 
technique.  Since  the  same  modelling  process  was  used  as 
with  unemployment  rates,  only  the  final  results  will  be 
discussed. in  this  section.  The  producer  price  index  was 
successfully  modelled  as  follows: 

( 1-B ) ( 1-. 2985B3  -  .4965B12 )  Xt  »  (1+.8789B  + 

.3169B5)  at 


or 

(1-.2985B3  +  . 2985B4  -  .4965B12  +  .4965B13)  x  _ 
(1  +  . 8789B  +  . 3169B5 )  t 

This  equation  was  used  to  filter  accessions  and  there 
were  three  transfer  function  models  that  were  found  to  suc¬ 


cessfully  model  the  data.  They  are: 


1.  Yt  *  (73.3B5-70.5B6)  Xfc  + 


(1  +  .9072B) 


( 1-. 5949B2-. 2420B6 )  t 


or 


Yt  -  73.3Xfc_5  -  70.5Xfc_6  + 


(1  +  .9072B ) 

( 1-. 5949B2- . 2410B6 )  t 


2.  Yt  =  Xt  + 


( 1-.7531B-. 1119B°  ) 


3.  Y.  =  (65.5  -  65. 2B)  X.  +  - - - r- 

Z  *  (1-.7202B  -.1522B*) 

or 

at 

Y.  *  65. 5X.  -  65. IX.  ,  +  - £ - =- 

*  *  C  1  (1-.7202B  -  ,1522Ij ) 


All  three  of  these  models  have  some  coefficients  that 

\ 

are  not  statistically  significant  at  the  .05  probability 
level.  This  causes  the  standard  error  bands  on  the  fore¬ 
casts  to  be  even  wider  than  those  in  the  unemployment  model 
Additionally,  the  only  model  that  allows  you  to  look  ahead 
with  current  values  of  the  PPI  is  the  first  and  the  lag  co¬ 
efficients  on  and  Xfc_6  are  not  significant.  The  other 

models  use  values  of  Xt  to  predict  Yfc.  This  defeats  the 
purpose  of  modelling  with  a  leading  indicator. 

Predictions  were  made  using  the  three  transfer  func¬ 
tion  models  and  the  results  are  in  Table  VII. 
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Table  VII 


P 


Comparison  of  PPI  Leading  Indicator  Models 


Period 


Ahead 

i 

Transfer 

Forecast 

Fn.  1 
SE 

Transfer 

Forecast 

Fn.  2 
SE 

Transfer 

Forecast 

Fn.  3 
SE 

1 

4936 

614 

4673 

619 

4874 

629 

2 

4739 

819 

4462 

762 

4716 

771 

3 

4870 

948 

4248 

867 

4732 

833 

4 

4724 

1037 

4127 

941 

4600 

862 

5 

4798 

1100 

3947 

995 

4522 

877 

6 

4691 

1149 

3821 

1037 

4640 

884 

The  best  two  standard  error  band  for  the  first  predic- 
tion  is  using  transfer  function  1  with  a  band  of  +  1228. 

For  a  list  of  BMDP  output  see  Appendix  E-7.  Once  again 
this  is  not  a  model  that  would  be  readily  adopted  by  Army 


I 

i 

« 

I 


planners. 

Prime  Rate  as  Input  Variable.  The  Prime  Rate  was  also 
chosen  as  the  input  variable  in  the  leading  indicator  model 
with  accession  rates  remaining  the  output  variable.  As 
with  the  PPI,  only  the  final  results  will  be  discussed. 

The  prime  rate  yielded  two  models: 

(1-B)(1-.3053B14)  Xt  =  (1+.6640B  -  .3625B3  -  .5574B4)  afc 
or 

(1-B-.3053B14  +  .3053B15)  Xfc  =  (1  +  .6640B  -  .3625B3  - 
•5574B4)  at 
or 

(1-B-.3053B14  +  .3053B15) _  X  =  a 

(1  +  . 6640B  -  . 3625B3  -  .5574B4)  t 


l[  and  the  second  model: 

{ 1-B ) ( 1- . 67 60B  +  .3685B2  +  .3756B4  -  .3729B5  +  .2573B10) 

f 
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(1-  1.6760B  -  1.0445B2  -  .3685B3  +  .3756B4  -  .7485B5  + 
•8421B6  -  .4692B7  +  .2573B10  -  .2573B11)  Xfc  =  afc 


Xfc  -  1.6760Xt_1  -  1.0445Xt_2  -  .3685Xt_3  +  .3756Xt_4 
,7485Xt-5  +  .8421Xt_6  -  .4692Xt_7  +  ,2573Xt_1()  - 

•2573Xt-ll  "  at 

where  in  both  models  Xt  is  the  input  variable «  prime  rate, 
and  afc  is  a  white  noise  residual. 

These  models  were  used  to  filter  accession  rates  and 
they  resulted  in  the  corresponding  transfer  functions: 


(1-  . 7502B  -  . 1048B- ) 


v  _  v  a.  (1  +  .  3455B 
Yt  t  +  (1  -  .7177B ) 


Neither  of  these  transfer  function  models  are  usable 
because  they  use  values  of  the  input  to  predict  the  same 
time  frame's  accession.  BMDP  computer  listings  for  prime 
rate  Leading  Indicator  Modelling  are  at  Appendix  E-8.  A 
summary  of  the  two  models'  predictions  for  accessions  rates 
is  at  Table  VIII. 


.1: 


Comparison  of  Prime  Rate  Leading  Indicator  Models 


Period 


Ahead 

Model 

Forecast 

1 

SE 

Model  2 
Forecast 

SE 

1 

4674 

613 

4720 

606 

2 

4356 

766 

4322 

744 

3 

4167 

840 

4176 

805 

4 

3970 

879 

3751 

834 

5 

3828 

900 

3582 

848 

6 

3788 

912 

3669 

856 

Once  again  the  two  standard  error  bands  are  too  wide 
on  the  first  predictions  for  either  of  these  models  to  be 
considered  acceptable. 


Results  and  Conclusions 

The  objective  of  modelling  using  a  leading  indicator 
is  to  find  a 'variable  that  leads  accession  rates  and  to 
provide  a  better  predictive  model  than  by  simply  modelling 
accession  rates  as  a  time  series.  Leading  indicator  models 
were  developed  where  accessions  lagged  behind  unemployment 
rates  at  1,  7  and  8  months  and  behind  Producer  Price  Index 
at  5  and  6  months.  The  unemployment  model  provided  a  tight¬ 
er  two  standard  error  band  than  by  modelling  accessions  as 
a  time  series.  However,  in  all  cases  the  standard  error 
bands  were  too  wide  to  provide  a  realistically  usable  model 
for  Army  planners.  The  conclusion  is  that  unemployment 
rates,  the  producer  price  index  and  the  prime  rate  are  not 
reasonable  input  variables  for  a  leading  indicator  model. 
There  are  many  reasons  for  this;  the  basic  assumption  un¬ 
derlying  all  time  series  modelling  is  that  the  basic 


a 


stochastic  process  inherent  in  the  data  is  not  changing. 
Currently,  unemployment  rates  sure  skyrocketing  and  the  prime 
rate  is  falling;  there  may  be  forces  at  work  that  lead  these 
variables  which  our  models  fail  to  account  for.  However, 
we  are  attempting  to  model  nationwide  accession  rates  with 
nationwide  unemployment  rates.  This  may  be  the  wrong  way 
to  model  accessions.  People  are  affected  in  various  parts 
of  the  country  by  any  number  of  influences  to  join  the  Army. 
The  unemployment  rate  is  one  such  influence.  In  the  South 
the  unemployment  rate  has  been  running  much  lower  than  in 
the  industrial  North  East  and  Mid  West.  One  would  expect 
unemployment  rates  to  have  a  different  effect  on  accessions 
based  on  the  local  unemployment  rate  rather  than  on  the 
national  rate. 
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We  are  led  to  the  obvious  conclusion  that  this  subject 
area  needs  to  be  disaggregated  and  studied  on  a  regional 
basis.  Additionally,  the  input  variables  should  be  consi¬ 
dered  as  output  variables  and  leading  indicators  for  these 
variables  (especially  unemployment  rates)  should  be  sought. 
There  are  a  multitude  of  economy  related  variables  to  be 
found  so  that  continuing  this  effort  with  other  input  vari¬ 
ables  should  not  be  a  difficult  task.  Finally,  there  seems 
to  be  problems  in  the  predictions  for  the  input  variables 
( they  all  seem  to  be  going  in  the  wrong  direction  from  the 
actual  predicted  values).  It  is  recommended  that  this  sub¬ 
ject  be  studied  by  applying  the  frequency  domain  approach. 
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VIII.  Summar 


Research  Objectives 

The  objectives  of  this  research  effort  were  as  follows: 

1.  Enhance  the  USAREC  Enlistment  Prediction  Model. 

2.  Study  accessions  in  an  attempt  to  understand  what 
factors  influence  them. 

3.  Build  a  discriminant  function  to  predict  whether 
potential  accessions  will  join  the  Army  or  other  services. 

4.  Conduct  time  series  analysis  on  accessions  using 
Box  and  Jenkins  techniques. 

5.  Develop  leading  indicator  models  to  predict  acces¬ 
sions  based  on  the  following  input  variables:  unemployment 
rates,  producer  price  index,  and  the  prime  rate. 

6.  Provide  to  ARI  and  USAREC  useable  computer  programs 
to  accomplish  the  above. 

Accomplishments 

The  following  are  the  major  accomplishments  of  this 
thesis: 

1.  Enhanced  Enlistment  Prediction  Model. 

The  predictive  ability  of  the  regression  model  used 
by  USAREC  in  projecting  accessions  has  been  significantly 
improved  through  the  incorporation  of  new  relationships  among 
the  USAREC  variables.  Running  this  EPM  model  is  the  first 
phase  in  USAREC' s  designation  of  recruiting  missions  for  each 
of  the  57  District  Recruiting  Commands  (DRC).  This  enhance¬ 
ment  allows  USAREC  to  formulate  more  accurate  quarterly 


projections  of  accession  rates  and  will  make  all  steps  after 
it  more  exact  with  the  goal  of  assigning  a  more  precise  and 
feasible  recruiting  mission  to  the  DRC's. 

A  comparison  of  the  adjusted  R  from  the  USAREC 
model  and  the  enhanced  model  is  shown  below: 


USAREC 

Enhanced 

Reaion 

Model 

Model 

NE 

.75220 

.76185 

SE 

.73807 

.76889 

SW 

.68286 

.71546 

MW 

.62862 

.63073 

WEST 

.75133 

.81231 

2.  District  Recruiting  Command  (PRC)  reorganization  of 
MEPS  file  and  its  merger  with  the  USAREC  file. 

One  usable  byproduct  resulting  from  this  thesis 
study  is  the  development  of  a  FORTRAN  program  that  will  re¬ 
code  state  and  county  codes  on  the  MEPS  file  into  DRC  codes. 
ARI-PERI-RP  has  been  very  interested  in  having  this  capabili¬ 
ty  in  order  that  they  might  be  able  to  perform  analysis  of 
the  MEPS  file  based  on  recruiting  regions  and  districts.  An 
additional  FORTRAN  program  allows  researchers  to  enhance  this 
DRC  analysis  capability  by  merging  the  information  on  the 
MEPS  file  with  the  data  from  the  USAREC  files.  Thus,  the  re¬ 
searchers  at  ARI-PERI-RP  will  be  able  to  organize  their  MEPS 
file  investigation  by  DRC  codes:  and  as  a  corollary,  a  bridge 
is  provided  for  the  exchange  of  information  between  ARI  and 


USAREC 


3.  Multivariate  analysis  techniques  for  examining  the 
variables  in  the  merged  MEPS/USAREC  (MU)  file. 

Based  on  the  profile  analysis  done  on  the  MEPS  file 
and  the  merged  MU  file,  several  characteristic  variables  were 
selected  and  analyzed  by  Recruiting  Region  and  DRC  using  both 
factor  and  discriminant  analysis.  The  variables  that  emerged 
from  the  profile  analysis  and  the  multivariate  procedures  and 
output  formats  used  in  this  study  can  be  employed  in  the 
continuing  research  efforts  of  both  AR1  and  USAREC. 

4.  Time  Series  Analysis  of  economic  variables  asso¬ 
ciated  with  recruiting. 

The  results  of  this  study's  time  series  analysis 
provides  USAREC  with  new  tools  for  examining  accession  rates 
based  on  certain  economic  indicator  variables.  Box-Jenkins 
time  series  analysis  techniques  were  applied  to  the  USAREC 
data.  In  addition,  unique  procedures  and  programs  were 
specifically  developed  in  this  thesis  effort  for  use  on 
several  economic  variables  available  to  USAREC.  Armed  with 
these  programs  and  procedures,  USAREC  will  be  able  to  carry 
out  extended  research  into  the  economic  and  demographic 
variables  affecting  recruiting  and  accessions. 

The  objective  of  the  time  series  analysis  was  to 
find  variables  that  lead  accession  rates.  This  was  done  in 
the  unemployment  model;  unemployment  rates  were  found  to 
lead  accessions,  i.e.,  coefficients  were  significant,  by  one. 


seven,  and  eight  months.  Also^the  producer  price  index  leads 
accessions  by  5  and  6  months.  However,  there  mciy  be  a 


problem  with  the  leading  indicator  models  since  the  standard 
error  bands  about  the  predicted  values  of  accessions  were 
large.  In  this  case  it  will  be  necessary  to  see  how  well 
the  predictions  fare  over  the  next  six  months. 

Recommendations  for  Further  Analysis 

1.  Extend  time  period. 

The  results  of  this  thesis  effort  open  doors  to 
several  areas  of  investigation  in  the  area  of  Army  recruiting 
The  chief  recommendation  is  that  the  many  procedures  which 
are  outlined  in  this  study  can  and  should  be  applied  against 
existing  data  bases  and  variables  over  an  extended  time 
period  rather  than  the  narrow  bands  of  time  which  were  cov¬ 
ered  in  this  thesis.  From  an  extended  interval  of  time, 
more  accurate  accession  indicators  and  predictive  coeffi¬ 
cients  can  be  developed.  Since  the  MEPS  sample  was  limited 
to  one  year  of  data,  many  USAREC  variables  did  not  have  much 
variance  in  only  four  quarters.  Thus,  shifts  in  such  vari¬ 
ables  as  unemployment,  percentage  of  Army  recruiters  by  DRC, 
and  number  of  blacks  available  cannot  yield  solid  predictive 
models  when  based  on  only  four  quarters  in  the  MU  file  anal¬ 
ysis. 

2.  Add  variables  to  data  base  and  extend  studies  to 
Regions  and  PRC's. 

Now  that  a  method  is  available  for  ARI  to  examine 


the  MEPS  data  using  DRC  codes,  other  demographic  and  economic 
variables  can  be  obtained  from  the  Bureau  of  the  Census  and 
merged  with  the  MEPS  data  by  DRC.  Further  expansion  of  the 


variables  and  analysis  in  this  area  can  now  focus  on  specific 
Recruiting  Regions  and  DRC 1 s  rather  than  be  limited  to  the 
relatively  meaningless  MEPS  ID  codes. 

3.  Perform  multivariate  analyses  on  other  service 
components. 

The  factor  and  discriminant  analysis  procedures 
developed  for  the  study  of  the  MU  file  can  be  employed  in 
further  research  of  the  variables  affecting  the  National 
Guard  and  Reserve  components.  The  data  sample  of  the  MEPS 
file,  although  very  large  at  48,520  records,  did  not  have 
enough  records  to  perform  detailed  analysis  at  the  DRC  level. 
Much  can  be  learned  from  applying  the  multivariate  analysis 
procedures  to  the  other  components  as  was  applied  to  the 
Regular  forces. 

4.  Extend  factor  analysis  with  correlation  analysis. 

The  next  step  in  factor  analysis  of  the  MU  file 

would  be  to  assign  standardized  names  to  the  discovered  fac¬ 
tors  and  then  run  correlation  analysis  such  as  the  SPSS  sub¬ 
routine  PEARSON  CORR  with  the  various  factors  against  a  de¬ 
pendent  variable  such  as  the  Army  variable.  The  results 
should  demonstrate  how  closely  the  new  factors  predict  the 
Army  variable,  that  is,  whether  an  individual  is  "oriented" 
to  join  the  Army  based  on  the  characteristics  of  the  new 
factors. 

5.  Perform  discriminant  analysis  on  entire  data  base 
and  study  variables'  "target"  groups. 

The  next  step  in  discriminant  analysis  of  the  MU 
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file  would  be  to  use  the  entire  data  base  rather  than  just  a 
sample  to  develop  accurate  unstandardized  canonical  correla¬ 
tion  coefficients  by  DRC.  The  coefficients  can  then  be  used 
to  determine  if  certain  "target"  groups  are  more  likely  to 
join  the  Army  or  another  branch  of  service.  This  is  benefi¬ 
cial  in  those  cases  when  a  certain  "target"  group  of  the 
population  is  sought  by  the  Army,  such  as  for  instance  high 
school  graduate,  male  hispanics  from  a  given  DRC.  In  this 
example,  by  knowing  the  current  unemployment  rate,  median 
income,  percentage  of  Army  recruiters,  etc.,  one  can  use  the 
unstandardized  coefficients  to  predict  whether  or  not  the 
"target"  group  will  join  the  Army  over  another  branch.  It 
could  be  that  decreasing  the  percentage  of  Army  recruiters 
0  only  slightly  in  the  equation  will  have  a  strong  impact  on 

the  "target"  group  deciding  to  join  another  branch  instead. 

Several  "target"  populations  with  widely  varying  attributes 

\ 

can  be  run  against  this  model  to  determine  their  likelihoods 
of  joining  the  Army.  Such  techniques  could  be  employed  when 
a  given  Region  or  DRC  is  planning  an  advertising  campaign  in 
a  given  DRC  to  beef  up  accessions,  e.g. ,  it  may  be  found 
that  the  best  advertising  for  the  Army  should  be  based  on 
the  income  variable  rather  than  the  unemployment  variable 
in  a  given  DRC. 

6.  Add  variables  to  discriminant  model. 

Other  variables  can  also  be  added  to  the  discrimi¬ 
nant  model  to  enhance  its  predictive  capabilities,  such  as 
the  addition  of  amount  of  advertising  dollars  spent  by  the 
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various  services.  The  user  of  these  procedures  should  be 
cautious  not  to  approach  the  use  of  these  variables  as  com¬ 
petitive  tools  to  draw  desired  "high  quality"  enlistments 
away  from  the  other  services.  Rather,  the  aim  is  to  attract 
the  same  caliber  of  personnel  into  the  Army  as  is  enjoyed  by 
the  other  branches,  and  to  avoid  expending  scarce  resources 
on  the  wrong  "target"  groups. 

7.  Consider  other  methods  of  preparing  the  data  base 
for  further  analysis. 

Similar  analysis  as  developed  in  this  study  can  be 
performed  on  a  different  mode  of  MEPS/USAREC  file  organiza¬ 
tion  and  file  merging.  One  suggestion  referred  to  in  Chapter 
V,  Data  Base  Preparation,  is  to  develop  summary  totals  within 
the  MEPS  file  by  DRC  and  quarter  based  on  tallies  of  the 
attribute  variables.  For  instance,  ore  may  tally  the  number 
of  individuals  with  AFQT  scores  over  50  by  service  by  DRC 
by  quarter,  or  tally  a  cross  combination  of  several  charac¬ 
teristics.  These  MEPS  file  tallies  could  then  be  merged 
with  the  USAREC  data  base  to  form  a  new  summary  file  to  be 
used  in  further  analysis.  This  merging  can  use  the  same 
FORTRAN  routine  for  state/county  to  DRC  recoding  that  was 
employed  in  the  original  MU  file  development. 

8.  Consider  enhancing  EPM. 

There  are  logical  extensions  to  the  enhancement  of 


the  USAREC  EPM.  In  Dale  and  Gilroy's  study  (Ref  5),  military 
pay  was  found  to  have  a  significant  relationship  with  acces¬ 
sions.  Inclusion  of  wage  differential  variable  or  a  percent 


annual  pay  raise  variable  in  the  model  might  further  enhance 
the  quality  of  its  output. 

9.  Perform  time  series  analysis  at  regional  level 
using  leading  indicators  and  spectral  analysis. 

In  the  time  series  analysis,  we  are  led  to  the 
obvious  conclusion  that  this  subject  area  needs  to  be  dis¬ 
aggregated  and  studied  on  a  regional  basis.  Additionally 
the  input  variables  should  be  considered  as  output  variables, 
and  leading  indicators  for  these  variables  (especially  un¬ 
employment  rates)  should  be  sought  after.  There  is  a  wealth 
of  available  economy-related  variables  to  be  considered,  so 
that  continuing  this  effort  with  other  input  variables  should 
not  be  a  difficult  task.  Finally,  there  seems  to  be  a  pro¬ 
blem  in  the  predictions  for  the  input  variables;  they  all 
seem  to  be  in  error  in  the  same  direction  from  the  actual 

values.  It  is  recommended  that  this  subject  be  studied  by 

\ 

applying  the  time  domain  approach. 

Compendium 

The  overall  importance  of  this  thesis  effort  is  that 
it  provides  a  variety  of  tools  for  the  study  of  the  many 
aspects  of  Army  Recruiting.  In  summary,  the  tools  developed 
and  employed  are: 

1.  State  and  county  recode  to  DRC, 

2.  Merge  of  MEPS  and  USAREC  data  bases, 

3.  Profile  of  MEPS  and  USAREC  variables  important  to 


recruiting, 
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PROGRAM  TSARS t IKPUT , OUTPUT , TAPE5= IHPUT . TftPE6=OUTPUT t TAPED 
COMMON  AC (51) 

COMMON  LH 
COMMON  NAC 
INTEGER  NAC > HP AC 
CHARACTER  IUK(88).IUHAT 
CHARACTER  TITLE*Z# 

REAL  R(4fi4f) iS(4«t4f ) ,Z(3L8) • ACV<5#> iPACV(Sfl) 

REAL  NK(5f) 

CHARACTER  A(-5#:5I) 

CHARACTER  R«1 
REAL  N 
REAL  CC 
CC-f.l 
IFLGM 

PRINT  •••HON  LONG  IS  THE  TINE  SERIES?1 
READ (S»«)  LZ 
09  REMIND  I 
LH-LZ 

DO  1  1=1, LZ 
READUi*)  ZII) 

1  CONTINUE 

88  PRINT* i '  ENTER  TITLE  • 

READ(5il#2)  TITLE 

N*f.f 

CC*#.# 

112  FORMAT (A2I) 

PRINT  •••DO  YOU  I1ANT  TO  DIFFERENCE  OR  TRANSFORM  <T/N1?* 

READ(S.IM)  INHAT 
IM  FORMAT  UAH 

IF l INHAT. ED,’ 1M  THEN 

PRINT  •••NHAT  IS  THE  ORDER  OF  THE  NON-SEASONAL  DIFFERENCE: ID! >=#  » 
READISi*!  101 

PRINT  •••NHAT  IS  THE  ORDER  OF  THE  SEASONAL  DIFFERENCING: ID2>=#  • 
READISi*!  IDZ 

PRINT* •’TRANSFORMATION  COEFF:  IP*#  ILN!,IP>#(ZU)*Z(I)**IP)» 
READISi*)  IP 

PRINT* • ’LENGTH  OF  SEASON:  IS)*#’ 

READISi*) IS 
LZA*LH 

CALL  FTRDIFIIDliIDZiIPiISiLZAiZiSHIFTiLNiIER) 

NRITE  (7i’(AS#!’l>  DIFFERENCED  DATA’ 

NR1TE  (7iMtXiFZf.il)’)  (Z ID i I*iiLH) 

ENDIF 

NAC*LN/Z 

IF  (NAC.CT.4S)  NAC =45 
NPAC*LN/2 

IF  (NPAC.GT.4S)  NPAC*45 
C 

CALL  F(AUTOiZiLHiNACiNPACi7iAMEANiVARiACVtACiPAC9iNK) 
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WRlTEU.tlZ)  TITLE 
IWITC47* i#3)  TITLE 
113  FORMAT!  INI.  1 1.  A2M 

PRINT*. »  DEMI  *  '.AHEAN.'  VARIANCE  -  '.VAR 
PRINT*.'  ' 

URITE<7ilfl)  AHEAN.VAR 

lfl  FORMAT (/>'  AHEAN  *'.F12.6»»  VARIANCE  >  » »FI2.6»/I 

PRINT*. 'DO  YOU  HANT  TO  SEE  THE  FT AUTO  OUTPUT  <I/NI?' 

READtS.INI  IMHAT 
IF  UNHAT.EQ.'Y')  THEN 

PRINT*.'  LAC  A!  i'OCOVARIANCE  AUTOCORRELATION  PARTIAL  AUTOCORRELATI 
ION* 

IF  (NAC.CT.4SI  NAC=45 
11=2.1/ (SORT (REAL(LU) ) I 
DO  2  I’l.NAC 

IF  (ABStAClin .CT.ABS4CC) )  CC=AC(I» 

IF  (ARSIACIin.CT.ABSIim  THEN 

B>i«' 

ELSE 
Is*  * 

ENDIF 

PRINT  7.I.ACVU>.ACU>»B.PACVm 
WRITE  (7.7)1  lACVd) lACil) 

2  CONTINUE  < 

PRINT  *.'  ' 

WRITE  <«.5Ml  'THE  APPRO! IMATEO  2  S.E.  MMD  18  */-Ml 
SM  FORMAT  (1I.A3S.F1I.4I 
PRINT  *.»  ' 

PRINT  «i'«-0EN0TE8  AUTOCORRELATION  OUTSIDE  2  S.E.  BAND' 

PRINT  «.'  CLEAR  SCREEN.  THEN  INPUT  ANT  CHARACTER  FOR  AUTOCORRELATI 
ION  PLOT' 

READ  («.21M>  B 
21M  FORMAT  (All 
PRINT  «.*  ' 

PRINT  «.'  • 

PRINT  (.TITLE 

PRINT  *.'  A 

1UTOCORRCLATION  PLOT' 

PRINT  «.'  • 

IF  (U.CT.ABS(CCI)  CC=«1 

N=ABS(CC)«1.1 

II=ABS(IU*Sf/N 

WRITE  <«»5B1>  'LAC’iABSlCCM(-l.l) . 

1 - - - • . 

I . ’.ABSICCMl.t 

511  FORMAT  (lI.A3.5X»F4.3.Alfl.F6.3) 

DO  2Bff  IsltNAC 
DO  2«1 
A(J) »'*' 

211  CONTINUE 


CCSACU)«5#.I/N 

CC*CC+.5 

ICC*IMT(CC» 

If  UCC.LT.f.)  THEM 
00  6H  J*-5I»ICC-1 
AU)*»  » 

Ml  CONTINUE 

00  MI  JM.5I 
A(JI«»  1 

Ml  CONTINUE 

ELSE 

If  (ICC.EQ.I.)  THEN 
DO  M2 
AUI*'  » 

M2  CONTINUE 

00  M3  JMiS# 

M3  CONTINUE 

ELSE 

If  UCC.CT.I.)  THEN 
00  M4  J=-5f •- 1 
AU)**  * 

M4  CONTINUE 

DO  MS  J=ICC-H>5f 
A<4>*»  ' 

MS  CONTINUE 

ENOlf 

ENOIf 

ENOlf 

A<-Il-.5»*'t* 

AUl*.5I»*l* 

UNITE  U.2fM>  It (A(J) ij=-5i»5l) 

20M  fORMT  (lX«I3tllXtl«lAl> 

2M  CONTINUE 

PRINT  «.»  .  » 

7  FORNATM, 12,61, F12. 6, U,F12.t,lA.  tl.F12.il 

PRINT  «»*  CLEAR  SCREEN.  THEN  INPUT  ANT  CHARACTER  fOR  PARTIAL  AUTOC 
10RRELATION  PLOT* 

READ  («.21M)  B 
PRINT  «,*  • 

PRINT  f,»  *  * 

PRINT  ••TITLE 

PRINT  «>*  PARTIA 

1L  AUTOCORRELATION  PLOT*  ) 

PRINT  «.*  > 

URITE  . . . 

1 . . . .  iN 

1501  FORMAT  IlIiA3iSIiF6.3iAlfl.F6.3I 
00  2M  1*1. NAC 
DO  21?  J>-S0.50 
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A(J!***> 

219  CONTINUE 

CC*PACV<ll*5l.l/N 
CC*CC*.5 
1CCSINT (CCI 
IF  (ICC.LT.I.l  THEN 
00  41 «  J=-5#» ICC- 1 
A(J)*»  ' 

All  CONTINUE 

00  411  J*t.5l 
A(J)>*  • 

411  CONTINUE 

else 

IF  (ICC.EO.f.)  THEN 
00  412  Jt-SI.-l 
AU)*»  • 

412  CONTINUE 

00  413  J*1.SI 
A(J)*'  » 

413  CONTINUE 
ELSE 

IF  (ICC.CT.f.l  THEN 
DO  414  4«-5l»-l 
A(J1*'  * 

414  CONTINUE 

DO  413  J*lCC+l«5f 
AU»«*  * 

413  CONTINUE 

ENQIF 

END  IF 

END  IF 

A4-I1-.5I=»I» 

A(IH.51*'I» 

NRITE  <*,2fH>  I. (A(J) > J='5i •  54) 

2M  CONTINUE 
.  PRINT  «,*  > 

IF(IFLC.EO.l)  CALL  BARTLT 
ENDIF 

C  FOURIER  DECOMPOSITION 

PRINT*. »  DO  fOU  RANT  THE  PERIODQCRAH? 
REAOIS.IMI  1UHAT 
IF  (INHAT. ED.'!* I  CALL  PCRAHIZ.LHl 
PRINT*.*  * 

C  DO  R  AND  S 

C 

DO  4  1*51.2, -1 
ACCIMACCI-U 
4  CONTINUE 
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ACdHl.f 

NXRQNMf 

HR0H*32 

IF  (Ui.LT.32)  HRQN=LN 
DEC- IS 

IF  (Ui.LT.3f)  MCC=LU/2 
KQL-ll 

44  CONTINUE 

PRINT*.*  ENTER  THE  VALUE  OF  IRHO  WHICH  FREO.-ULONI? 

REA0(5.«)  IRHO 

IRS-* 

CALL  RANDSiHXROH.HRONiNEC iKOL. IRHO. AC.R.S. INK. IRS) 

PRINT*.'  DO  YOU  KANT  ANOTHER  4  AND  S  ARRAY?* 

READ(S.IN)  INHAT 

IF  (INHAT.EO.'Y')  CO  TO  44 

PRINT*.'  DO  YOU  UANT  THE  D  STATISTIC?  ' 

READ (5.1*1)  INHAT 
5489  IF  (IHHAT.FVYM  THEN 

PRINT*. •  ENTER  THE  DAI  AR  ORDER  ' 

R£AD<5.«)  IP 

PRINT*.*  ENTER  THE  HAX  HA  ORDER  > 

REA0(5.«)  10 

PRINT*.'  DO  YOU  NANT  TO  SEE  THE  D  STAT?' 

PRINT*.*  f:YES.  l*NOs*>* 

REA0(5>«)  IPNT 

CALL  DSTATINIROH.HROH.NEC.R.S.IP.IO.IPNT) 

END  IF 

PRINT*.',  DO  YOU  NANT  ANOTHER  D  STATISTIC?  ' 

READ (5.  IN)  INHAT 
IF  (INHAT.EO.'Y')  GOTO  5489 
PRINT*.'  00  YOU  NANT  FTRXL?' 

B£AD(5>tM)  IUHAT 

5498  IF  (INHAT.EO.'Y')  CALL  OURHXL(Z.LN) 

PRINT*.'  DO  YOU  NANT  ANOTHER  FTRXL?' 

READI5.1N)  INHAT 
IF  (INHAT.EO.'Y*)  C0T05498 
PRINT* •'  DO  YOU  NANT  ANOTHER  RUN?' 

REA0I5.1N)  INHAT 
IF  (INHAT. HE. ’ Y'l  STOP 

PRINT*.'  WHICH  DATA?  ^ORIGINAL.  1*CURRENT**>' 

R£AD(5.»)  IANS 

IF  (IANS.EO.I)  C0T089 

IF  (IANS.NE.I)  C0T088 

STOP  ’ 

END 


SUBROUTINE  BARTLT 

REAL  R(-lfiM).VR(tf).STDUf) 

INTEGER  V 
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COMMON  AC (51) 

COMMON  LU 
COMMON  MAC 
BO  1  Is* ll»l 
RU)*«.I 

1  CONTINUE 

DO  2  IM.NAC 
Rtn=BCtn 

2  CONTINUE 

DO  3  !sNAC+li6f  , 

R(I)*l.f 

3  CONTINUE 

C  REMEMBER  R(l)=l 

R(BI*1.B 
C 

C  BY  BARTLETT'S  APPROXIMATION)  EOU  2.1.11 

DO  5  K*ltll 

1=1.1 

DO  i  V:I.NAC 

M4R(V)**24R(V4K)*R(V-KM*R(I()*R(V)*R(V-K)42«R<V)«*2*R{|0* 

1*2 

6  CONTINUE 

VR(KI=I/FLQAT(LH) 

STOCK) sSORT  CVRCKI ) 

5  CONTINUE 

PRINT*.'  ' 

PRINT*.'  ' 

PRINT*.'  FIRST  LINE  VAR JRUUl. SECOND  LINE  STDIRIKM' 

■  PRINT*.'  » 

PRINT*.'  1  2  3  ». 

1'  4  5  i  '. 

2'  7  8  9  IB' 

PRINT  23. CVR(K) .K=l»lf 1 
PRINT  23. (STD (K). KM. If) 

23  FORMATUI.lflFU.in 

RETURN 
END 

SUBROUTINE  RAMDS4HXROU.HROU.NEC.KOL* IRHO.C.R.S* INK. IRS) 

DIMENSION  R(MXROU.t) .S(MXROU.l) .0(1) » IUKU) 

CHARACTER  INHAT 

DATA  1BLK/1H  /.1STR/1H*/.1ZR0/B/ 

K0LH1 3KOL- 1 
DO  5  ICM.KOL 
DO  5  IRM.HRON 
RIIR.IC)>«.« 

5  StIR.  10*1.1 
HN»HRQN-NEC 
DO  11  1*1. NN 
NEC1*NEC4I 
8(NEC1.1)*1.« 
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if  aucciiti'CM) 

BO  2«  1*1. NEC 
NECI*NEC*2-I 
Sll. I)*l.l 
21  R(l.l)*C(NECll 

IF  (MHO  .Efl.  •)  C0T04I 

bo  ai  im.hrou 

NEC1M-NEC-1 

31  R(I»l)*R(I»l)*(*l)**(NECI) 

41  J1*NR0H 
BO  71  IC>2*K0(. 

1011*10-1 
Jl*Jl*t 
00  41  1RM.JI 
IRP1MR+1 

IF  (AIS(R(lR»ICHl)t  .£0.  f.f)  COTOSf 

S(IR>IC)*S(IRPl.lCHl)*(R(IRPl.!CRl)/R(lRtICRl)*1.0) 

COTOM 

54  S(1R.1C)**S(1RP1»ICH1) 

41  CONTINUE 
Ji-Ji-l 
00  71  IR*tiJl 
IRPIMRM 

7f  R(IR>IC)*R(IRPl*ICHl)*(S(IRPl»IC)/S(lR.IC)*l.f) 

IF  (IRS.NE.f)  RETURN 
75  FQRMTM1I 

PRINT*.*  THE  4  AND  S  ARRAYS  HAVE  BEEN  CRUNCHED.' 

PRINT*. *  DO  YOU  ACTUALLY  HANT  TO  SEE  THEN?' 

READ15.75)  IHHAT 

IF  MNHAT.NE.'TM  RETURN 

IF  (IRHO  .£0.  II  COTOM 

HRITEU.2M) 

HRITE47.2MI 

COTOM 

M  URITE(4.21f) 

UR1TE(7.21I) 

M  NRITE(A.ZZf)  (1C.IC*1,K0L)  * 

HRITE(7.22f)  (IC.ICM.KOU 
DO  ttf  IRM.HROU 
IRN*IR-NEC-l 
DO  1M  1C*I.K0L 
IHKIICMIBU 

Iff  IFIIR.E0.NEC+2-IC)  IHKMOMSTR 

URITE(6.23f )  IRN. (R(IR.K) . IUK(K) >K*1 .KOI) 
HRITE(7»23f)  IRN.  (R(IR.K) . IHK(K)  .KM  »KOU 
ill  CONTINUE 

PRINT*.'  TTPE  ANY  CHARACTER  TO  CONTINUE**)' 

READ (5. 75)  IHHAT 

UNITE (4. 241)  IZR0.UC.ICM.K0UI1I 
HRITE(7.24f)  IZROt (1C. IC-t tKOLfll) 
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BO  13#  IRd.HRON 
IRN-IR-NCM 
DO  12#  XC<i>KOL 
IHK(IC)M8LK 

IF  IIR.EQ.NEC+2-IC)  IHKdCMISTR 
12#  IF  I1R.E0.MEC^3-1C>  IUKUO-1STR 

HRITE(A.23#)  IRNi(S(IRiKMUK(K)iK=1iK0L) 

HRITE(7*23f)  1RM* (S( IR.K)  •  IUK(K)  iK-1  .KQL) 

13#  CONTINUE 
RETURN 

2M  FORMAT (//2IX*27H(  FIR)  =  (<-l)»«H)*ACF(H>  )) 

21#  FORMAT  I //25X • 17H (  FIR)  *  ACF(fl)  t) 

22#  FORMAT I//3#X 1 7HR- ARRAY //3H  Kilt(9X»I2>> 

23#  FORMAT (2H  (.I3.1H)ill(lXiF9.4.Al)) 

24#  FORHAT I //3#X *7HS- ARRAY/ /3H  K. 11(91. 12)) 

END 

SUBROUTINE  OSTAT (HXROU>HROM»NEC>R>Si IP. IOtIPNT) 

C 

C  . 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  H1ROU  -  RON  DIHENSION  OF  R  AND  S  ARRAYS  IN  CALLING  PROGRAM 
C  ,  (INPUT) 

C  IKON  -  NUMBER  OF  RONS  TO  NH1CH  R  AND  S  HERE  CALCULATED. I INPUT) 

C  NEC  -  NUMBER  OF  RONS  OF  NEGATIVE  LAC  IN  R  AND  S  ARRAYS! INPUT) 

C  R  -  R  ARRAY.  DIMENSIONED  HXROH  BY  K  IN  CALLING  PROGRAM* 

C  HHERE  K  IS  GREATER  THAN  OR  EQUAL  TO  KOI. (INPUT) 

C  S  -  S  ARRAY.  DIMENSIONED  AS  R  ARRAY  IN  CALLING  PROCRAH. 

C  (INPUT) 

C  IP  -  (INPUT) MAXIMUM  ORDER  OF  AUTOREGRESSION  PERMITTED. 

.  C  (OUTPUT)  ORDER  OF  AUTOREGRESSION  SELECTED. 

C  10  -  < INPUT )HAK IHUM  ORDER  OF  MOVING  AVERAGE  PERMITTED. 

C  (OUTPUT)  ORDER  OF  AUTOREGRESSION  SELECTED. 

C  IPNT  •  PRINT  CONTROL  PARAMETER  (INPUT) 

C  IPNT  *  «  0  STAT  VALUES  ARE  PRINTED 

C  IPNT  NE  0  PRINTING  IS  SUPPRESSED. 

C  ,HITH  IPNT  *  f  THE  WRITE  STATEMENTS  REQUIRE  THE  INPUT  VALUE 

C  OF  10  TO  IE  LESS  THAN  TEN.' 

C 

C 

C  NOTE 

C  NEC  MUST  BE  GREATER  THAN  OR  EQUAL  TO  IP+IQ+3*  FOR  IP  AND 

C  10  THEIR  INPUT  VALUES.  FURTHERHOREHRON  RUST  BE  GREATER  THAN 
C  OR  EQUAL  TO  NEC* IP* 10+4*  FOR  IP  AND  IQ  EQUAL  TO  THEIR  INPUT 

C  VALUES.  IF  (!P*1)*(IQ*1)  IS  GREATER  THAN  S#(HHER£  IP  AND  IQ  ARE 

C  AT  THEIR  INPUT  VALUES)  THEN  THE  DIHENSION  OF  STAT  MUST  BE  CHANCED 
C  ACCORDINGLY. 

C 

C 

C . 
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c 


'DIMENSION  R<RXR00it)iS(KXRQUti)iSTAT{5f) 

DATA  IPR/AMZRO/f/ 

JPPUPU 
JQQMQU 
IOP=JPP«JOO 
DO  4#  tMiJPP 
JP=I- 1 

DO  41  i*UM 
JQ=J-1 

JK*J00«<M>4J 
ISUD2=NEC-JP*JO+2 
I SUil ‘ISOI2-1 
iptupn 
INI*I-I 

S0T*<SUSUB2>I)4S(ISUBl>IPtn««2 
IHI.CT.il  CO  TO  2f 
DO  If  Jl*l .3 
NECH*NEC-JP-JQ-JI*1 
MECP=M€C- JP+ JQ+ JI+1 

if  BOT’IOT+R (HECK* 1 I **2*R (NECPi l  )«2 

CO  TO  35 

2f  ISUB3MSUR2M 

BOT*BOT/ <SUSUB3i  INI )*S(  ISUB2i  1 1 1 «2 

DO  3f  JI«I.3 

NECN*NEC- JP- JO- Jl*  i 

NECP’NEC-JP+JO+JIM 

NECN1*NEC(H1 

NECP1*NECPM 

3f  D0T*B0T4<R(NECN.[WR(NECHi»INin««24(R(NECP.IWR{NECPi.INl)l«tZ 
35  NECf*NEC-JP-JQ 
NECI*NEC#M 

ARC>S(NECf . IP1 I / (BOT«S(NECi » II I 
4f  STAT  4JK)*MS(ARC) 

INXM 

XMI*STAT(il 
DO  5f  1*1, W 

IF (STATUI.LT. IM)  M  TO  54 
HUM 

im*STAT(ll 
Sf  CONTINUE 

IP-(INI-il/JQQ 
IOMNX-(IP*JM)-t 
IFUPNT.NE.fi  RETURN 
URITEUPRiiff) 

URITEI7* iff I 

JQQsJQQ-t 

URITE47*  1  if  I  IZROtUilUiJM) 

HRITEUPRiltf)  IZROi(IiI*l»JQfl) 

DO  4f  JM.JPP 
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JNWJ-I 

IST=JM1*(J08U)+I 

IENBMST+JM 

UftiTC(7*12fl  JHt . (STAT ( 1 1 » I - 1ST .  IEND) 
hi  WRITE (1PR. 128)  JM1  >  (STATU)  >lsIST»  IEND) 

WRITE (7. 1381  IP.  IQ 
WRITEUPR.13#)  IP.IO 
RETURN 

IN  FORNAT(///5II.*D  STATIST!C*//29It 'ORDER  Of  MAM 
111  FORMAT*’  ORDER  OF  AR  MK4M3.4DI 
121  F0RNAT(/41.i2i4I»lf<lIt£I8.4>) 

131  FORMAT*/’  ORDER  OF  AUTOREGRESSION  SELECTED  ’>14/ 

C’  ORDER  OF  MOVING  AVERAGE  SELECTED  ’>14) 

END 

SUBROUTINE  PCRAMtZ.LW) 

DIMENSION  Z(348) 

REAL  1NTNCTI 
SINT *8. 8 
PI *3. 141592454 
PRINT* i •  LH>  ' iLU 

PRINT*. »  MOW  MANY  INTENSITY  VALUES  00  YOU  WANT  TO  SEE?* 
PRINT*.*  •  *  DEFAULT  OF  451  N  .LT.  LU/2  *«>’ 

READtS.O  N 
IF  (N.EQ.fl  N-45 
IF  (N.CT.Ul-ll  N=45 
WRITE (4. 78) 

URITE(7.78i 
00  21  IM.N 
A*#. 

1*1. 

BO  II  ITM.LW 

YARC=2 .8*PI*flOAT  * 1 1 (FLOAT (IT! /FLOAT (LUI 
A*A  ♦  ZIITXCQSUARC) 

B=l  *  ZIITXSINUARC) 

If  CONTINUE 

A*2.*A/FL0AT(LW) 

B-2 .*B/FLOAT  <LU) 

INTMCTY  * (FLOAT  (Ui)  /2 . ) *  ( A* A  *  D*B) 

FO*l.f*I/LU 

P’l.l/FQ 

WRITE (4.77)  I .FQ.P. INTMCTY 
WRITE(7.77>  I.FO.P.INTNCTY 
SINT  «SINT-*  INTNCTY 
21  CONTINUE 

77  FORMATdX. I3t2l*FS.3t2X«F7.3.2I .FI6.3Y 

78  FORMAT  421# ' I • .51. *F8» .41 » ’P’ . 151 . * INTMCTY ’ I 
PRINT*.’  ’ 

WRITE  (4.2431)  SINT/N 
WRITE  (7.2431)  SINT/N 
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2431  FORMAT  (Mi* AVERAGE  INTENSITY*' «F2I. ID) 

RETURN 

END 

4 

SUBROUTINE  OURNlUZiLN) 

DIMENSION  1(348) 

INTEGER  IND(8)iI£R 

REAL  ARPS (41 iPNAS(4) >CR(Z4) iAA(74i) 

CHARACTER  INHAT 
1ND(I)*LU 
IND(5)-5M 
I NO  <  4  >  =  4 
INO (7) =0 
INS (81 >5 

PRINT*.'  INPUT  P>  ORDER  OF  Afi**>» 

READtS.*)  IP 
IN0(2) SIP 

PRINT* •'  INPUT  0.  ORDER  OF  HA**>' 

READ(S.«)  10 
IN0(3)*IQ 

PRINT*.'  INPUT  D.  ORDER  OF  DIFF**>» 

RCA04S.*}  IS 
IND(4)*1B 

PRINT*.'  INPUT  S.  NUMBER  OF  SIGNIFICANT  DIGITS  REQUIRED"-)' 

REAS(S.«I  IS 

IN0(41*IS 

PRINT*.'  INPUT  I. 'ITERATIONS  REQUIRED  HM  CHANGE  TO  OICITS**)' 
READ(S»«)  II 

PRINT*.'  DO  YOU  UANT  TO  INPUT  ESTIMATIONS  FOR  THE  COEFICIENTS?' 
REAS(S.IN)  INHAT 
IF ( INHAT. Efl.'Y'l  THEN 
1ND(7)*1 

IF(IP.EQ.f)  CO  TO  152 

PRINT*.'  INPUT  AR  ESTIMATES  INDIVIDUALLY' 

BO  151  1*1. IP 
PRINT*.'  V 
READ (5 i*)  ARPS(I) 

151  CONTINUE 

152  IF(IO.Efl.f)  GO  TO  154 

PRINT*.'  INPUT  HA  ESTIMATES  INDIVIOUAU!' 

BO  153  1*1.10 
PRINT*.'  ?» 

REA0(5.«I  PHASUl 

153  CONTINUE 

154  ENOIF 

IN  FORNAT(lAi) 

CALL  FTMlL(ZtlNOiASPS.PMAS.PHAC.UNViGR.AA.IER) 

INI  FORMAT (/.'  IND(S)*  '.13.'  IER*  >»I3./. 

*  NNV*  • .F8.2. '  PNAC*  '.FI2.4I 
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1M2  FORMAT </»•  ARPS  VALUES*) 

1N3  FORMAT C / 1 *  PNAS  VALUES*) 

URITEU.1M1)  IMO(S)  t  ICRtUNVtPHAC 
URITE(7>tMi)  INO(S) •  1ER*UMV»PHAC 
URITEU.IM2) 

URlTE<7tlN2) 

00  2f  1*1. IP 

UR1TEIA.*)  ASPS) I) 

URITE(7.f)  ARPS(I) 

21  CONTINUE 
URlTE(6»lfl3) 

URITE(7.tM3) 

00  41  1*1.10 

URIT£U,«)  PflAS(l) 

UNITE (7.*!  PNASd) 

41  CONTINUE 
RETURN 
END 
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Enlisted  into  military  service. 

Enlisted  into  military  service  and  shipped  to 
Basic  Training  location  or  duty  station. 
Shipped  to  Basic  Training  location  or  duty 
station. 

Qualified  but  not  yet  enlisted  into  military 
service. 


Entry  Status 


Enlistment  without  delayed  entry  status.  That 
is.  Enlistment  or  induction  into  the  active  force 
without  delay  status,  or  reservist  or  National 
Guard  (not  delayed  enlistment  program  enlistee) 
ordered  to  active  'duty  or  enlisted  for  active 
duty  or  reserve  enlistment. 

Enlistment  from  delay  program.  That  is,  enlist¬ 
ment  in  active  force  from  an  authorized  delay 
program  to  include  returning  delayed  reservists. 
Reservist  ordered  from  delay  program.  That  is, 
reservists  ordered  to  active  duty  from  a  delay 
status. 

Enlistment  into  delay  program.  That  is.  Enlist¬ 
ment  into  an  authorized  delay  program,  such  as 
delayed  enlistment  programs. 
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FINAL  HOBCL  0C  OCCION  HOKl  9 


fIMl  ftQIU  W  KSJO«  H09EL4 


FIML  NOBEL  SB  RESIO*  SOUL  4 


flMl  M*£l  tf€5T  Hiit*  Mill  4 


•AIQESI 


STATISTICS  ALL 


Accea 


ation 


APQfP 

Age 

AirP 

Army 

ArmyR 

Black 

BMA 

ColDg 

DODHS 

DODNP 

Hi  spa 

HSDip 

HSS 

Incom  * 

Marin 

Marri 

Month 

Navy 

NoDip 

Prior 

QMA 

Sex 

Unemp 

White 

#MHSS 


Multivariate  Analysis  Variables 
DODNP  &  DODHS  (See  definitions  below) 

Armed  Forces  Qualification  Test  Score 
Age  of  applicant 
Air  Force  applicant 
Army  applicant 

Army  Recruiters  as  percent  of  DOD  by  DRC 

Black  applicant 

Black  Military  Available 

Associates  Degree  or  higher 

DOD  High  School  Diploma  Male  Accessions 

DOD  Non-Prior  Service  Male  Accessions 

Hispanic  applicant 

High  School  Diploma,  less  than  College  Degree 

High  School  Senior 

Median  Disposable  Family  Income 

Marine  applicant 

Married  applicant 

Month  of  earliest  transaction  on  applicant 

Navy  applicant 

No  High  School  Diploma 

Prior  Service  applicant 

Qualified  Military  Available 

Male,  Female  applicant 

DRC  Overall  Unemployment 

White  applicant 

Number  of  Male  High  School  Seniors 
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SWRRC 


NERRC 


SERRC 


1A  Albany 

3A  Atlanta 

4A  Albuquerque 

IB  Baltimore 

3B  Beckley 

4C 

Dallas 

1C  Boston 

3C  Charlotte 

4D 

Denver 

ID  Concord 

3D  Columbia 

4E 

Houston 

IE  Harrisburg 

3E  Jacksonville 

4F 

J  ackson 

IF  Ft  Monmouth 

3F  Louisville 

4G 

Kansas  City 

1G  New  Haven 

3G  Miami 

4H 

Little  Rock 

1H  Newburgh 

3H  Montgomery 

41 

New  Orleans 

11  Long  Island 

31  Nashville 

4J 

Oklahoma  City 

1J  Niagara  Falls 

3J  Raleigh 

4K 

San  Antonio 

IK  Philadelphia 

3K  Richmond 

1L  Pittsburgh 

3L  San  Juan 

IN  Syracuse  ,  . 

NWRRC 

WRRC 

5A  Chicago 

6A  Salt  Lake  City 

5B  Cincinnati 

6E  San  Francisco 

5C  Cleveland 

6F  Honolulu 

5D  Columbus 

6G  Los  Angeles 

5E  Des  Moines 

6H  Phoenix 

5F  Detroit 

61  Portland 

5H  Indianapolis 

6J  Sacramento 

51  Lansing 

6K  Santa  Ana 

5J  Milwaukee 

6L  Seattle 

5K  Minneapolis 

5L  Omaha 

5M  Peoria 

5N  St  Louis 

ral  Codes  Used  in  Multivariate  Analysis 


Regular  Forces 


Reserve  Forces 


Collective  Code* 


PS,F 


National  level,  all  Regions,  all  Components 
all  Regions,  one  Component 
Prior  Service  and  Female  records  included 
in  analysis 


Lonal  Codes 


1  **  NE 


5  «  MW 


Northeast  Region  Recruiting  Command  (NERRC) 
Southeast  Region  Recruiting  Command  (SERRC) 
Southwest  Region  Recruiting  Command  (SWRRC) 
Midwest  Region  Recruiting  Command  (MWRRC) 


WEST  *  Western  Region  Recruiting  Command  (WRRC) 
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Factor  Analysis  Table  Codes 


the  DRC  (or  Region)  row  is  associated  with  the  factor 
column,  i.e. ,  the  DRC  or  Region  in  the  given  MX"  row 
possesses  the  particular  factor  in  the  column  above 
the  ”XH . 

,  or  —  =  a  given  variable  in  the  signed  row  is 

associated  with  the  factor  in  the  signed  column. 

50  to  99  percent  of  the  variable  is  like  the  positive 
pole  of  the  factor. 

25  to  50  percent  of  the  variable  is  like  the 
positive  pole  of  the  factor. 

50  to  90  percent  of  the  variable  is  like  the  nega¬ 
tive  pole  of  the  factor. 

a 

25  to  50  percent  of  the  variable  is  like  the  nega¬ 
tive  pole  of  the  factor. 


FACTORS 

2  3  .4  5  6  [709012 


DRC  6K  81 
DRC  6l  80 


FACTORS 


Discriminating  Power  of  Discriminant  Functions 
(No  Accession  Variables) 


a 


Discriminating  Power  of  Discriminant  Functions 
,  w  (No  Accession  Variables) 

Ja 


DRC 

Canonical 

Correlation 

Wilks' 

Lambda 

Degrees  of 
Chi -Squared  Freedom 

Significance 

a 

1A 

.536 

.713 

38.41 

5 

.0000 

L 

.  ■ 

IB 

.411 

H 

00 

• 

50.34 

7 

.0000 

¥ 

1C 

.365 

•866 

28.52 

6 

.0001 

1 

ID 

.459 

.789 

20.74 

7 

.0042 

» ' 

IE 

.345 

.881 

26.11 

5 

.0001 

IF 

.454 

.794 

32.84 

5 

.0000 

H 

1G 

.466  • 

.783 

53.22 

5 

.0000 

1H 

.464 

.785 

27.51 

3 

.0000 

V* 

11 

.529 

.720 

44.98 

8 

.0000 

i 

CD 

1J 

.307 

.906 

13.83 

3 

.0032 

IK 

.401 

.839 

37.85 

7 

.0000 

» 

X- 

1L 

.443 

.804 

31.38 

4 

.0000 

■ 

IN 

.468 

.781 

36.91 

8 

.0000 

>; 

3A 

.460 

.788 

54.90 

6 

.0000 

-■? 

30 

.616 

.620 

20.52 

6 

.0022 

3C 

.525 

.724 

46.60 

7 

.0000 

3D 

.524 

.726 

56.11 

4 

.0000 

3E 

.423 

.621 

61.19 

10 

.0000 

i 

3F 

.351 

.877 

15.44 

7 

.0307 

* 

3G 

.474 

.776 

38.34 

8 

.0000 

3H 

.455 

.793 

39.22 

6 

.0000 

€ 

31 

.324 

.895 

17.08 

6 

.0090 

3J 

.461 

.788 

29.14 

6 

.0001 

3K 

.355 

.874 

22.22 

3 

.0001 

N  177 


Canonical  Wilks ' 
RC  Correlation  Lambda 


Degrees  of 
Chi -Squared  Freedom 


Significance 


t 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


Canonical 
DRC  Correlation  Lambda 


Degrees  of 
Chi-Squared  Freedom 


Significance 


a 


Discriminating  Power  of  Discriminant  Functions 
(With  POD  Accession  Variables) 


Canonical 

Wilks' 

Degrees  of 

DRC  Correlation 

Lambda 

Chi -Squared  Freedom  Significance 

P 

IF 

.457 

.791 

32.34 

5 

.0000 

* 

11 

.551 

.697 

49.11 

10 

.0000 

\  V 

IN 

.485 

.765 

39.75 

9 

.0000 

■ 

3A 

.492 

.758 

63.8I 

8 

.  .0000 

3E 

.422 

.822 

60.94 

9 

.0000 

31 

.334 

.889 

18.11 

7 

.0115 

3L 

.553 

.694 

19.17 

5 

.0018 

•;> 

4C 

.399 

.841 

18.19 

6 

.0058 

4D 

.328 

.892 

15.54 

5 

.0083 

■ 

hS> 

TO 

.542 

.706 

34.66 

5 

.0000 

4F 

.447 

.800 

23.37 

6 

.0007 

4K 

.444 

.803 

20.23 

6 

.0025 

i 

5F 

.479 

.771 

46.29 

8 

.0000 

*  '  a 

5H 

.487 

.763 

55.03 

7 

.0000 

* 

5J 

.431 

.814 

26.25 

7 

.0005 

5K 

.392 

.847 

25.31 

10 

.0048 

»  /- 

5L 

.386 

.851 

19.48 

4 

.0006 

5M 

.395 

.844 

28.60 

4 

.0000 

«  . 

5N 

.500 

.750 

54.87 

8 

.0000 

’"V 

6a 

.503 

.747 

26.73 

7 

.0004 

6L 

.505 

.745 

29.90 

9 

.0005 

« 
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Dlscrlalnatli 


.Power  of  Discriminant  Functions  by  Region 
(No  Accession  Variables^ 


Canonical Wilks*  Degrees  of 

Correlation  Lambda  Chi-Squared  Freedom  Significance 


All  PS,F 

.352 

.876 

1903.5 

15 

0 

All 

.333 

.889 

1260.2 

14 

0 

REG  all  Pd,F 

.378 

.857 

1766.0 

16 

0 

REG  all 

.333 

.889 

1020.6 

13 

0 

RSV1  PS.F 

.493 

.757 

145.4? 

11 

.0000 

RSV1 

.468 

.781 

99.57 

9 

.0000 

RSV3  PS ,F 

.645 

.584 

147.32 

11 

.0000 

NSV3 

.600 

.639 

89.43 

8 

.00O0 

RSV4  PS  ,F 

.520 

.730 

63.86 

10 

.0000 

RSV4 

.459 

.789 

35.74 

8 

.0000 

RSV5  PS.F 

.577 

.667 

175.24 

10 

0 

RSV5 

.569 

.676 

132.53 

7 

.0000 

RSV6  PS.F 

.61? 

.619 

115.47 

15 

.0000 

RSV6 

.530 

.719 

65.10 

11 

.0000 

REG1  PS,F 

.380 

.855 

443.98 

13 

0 

REG1 

.332 

.890 

253.97 

10 

0 

REG3  PS.F 

.430 

.815 

525.72 

13 

0 

REG3 

.393 

.845 

318.80 

8 

0 

REG4  PS.F 

.380 

.856 

231.95 

12 

0 

REG4 

.373 

.861 

168.31 

8 

0 

REG5  PS.F 

.346 

.880 

373.20 

9 

0 

BGC5 

.297 

.912 

207.03 

11 

0 

HBG6  PS.F 

.337 

.866 

179.99 

12 

0 

REG6 

.26? 

.929 

83.14 

6 

.0000 

183 


Vi 

VI 


Discriminating  Power  of  Discriminant  Functions  by  Region 


(With  POD  Accession  Variables) 


Canonical  Wilks 


Degrees  of 


Correlation  Lambda  Chi-Squared  Freedom  Significance 


.380  .855  178?. 8  17  0 

.336  .887  1038.9  12  0 


150.70 

13 

.0000 

14?.32 

11 

.0000 

176.32 

11 

0 

263.00 

9 

0 

528.79 

11 

0 

323.21 

10 

0 

213.53 

11 

0 

87.66 

9 

.0000 

% 


s. 

STANDARDIZED  CANONICAL  DISCRIMINANT  FUNCTION  COEFFICIENTS 
(NO  ACCESSIONS  VARIABLES) 


DRC 

u  Hi 8 pa 

Marri  #iiHSS 

Unemp 

ArmyR 

Inc  ora 

Centroids 

Non-Army 

* 

Army 

1A 

-.75 

-.417 

.950 

IB 

-.18 

.27 

-.390 

.518 

1C 

-.35 

.30 

-.279 

.546 

ID 

-.34 

1.03 

-•353 

.741 

is 

-.49 

-.251 

.535 

IF 

.21 

.27 

.31 

-304 

.841 

1G 

-.20 

-.359 

.764 

1H 

.24 

-.291 

.926 

11 

-.39 

-.32 

-.440 

.871 

1J 

-.224 

.458 

IK 

-.38 

.43 

CM 

• 

• 

.740 

1L 

.28 

-.329 

.730 

IN 

.23 

-1.87 

-.396 

.700 

3A 

.15 

.16 

-.538 

.494 

3B 

* 

.60 

-.734 

.798 

3C 

.34 

-.73 

-.500 

.751 

3D 

.38 

-.575 

.650 

3B 

.51 

.33 

.62 

-.406 

.534 

3F 

.47 

.51 

-.40 

-.363 

.381 

3G 

-.30 

.26 

.66 

1.03 

-.403 

.708 

3H 

.20 

-.422 

.612 

31 

-.42 

-.258 

.449 

3J 

.28 

-495 

.536 

% 

-.331 

.431 

188 

DRC  Hi spa  Marrl  #MHSS  Unemp  ArmyR  Incom  Non-Army  Army 


31 

.29 

-1.358 

.289 

4A 

.46 

-.46 

-.425 

.552 

4C 

-.57 

-.85  1.78 

-.349 

.543 

40 

-.52 

i 

. 

ro 

& 

.374 

4G 

.36 

.42 

-.31 

-.531 

.784 

4F 

.56 

-.384 

.478 

4G 

-.365 

.5^8 

4H 

-.32 

-.571 

.690 

41 

.45 

-.362 

.686 

4J 

-.259 

.280 

4K 

-.29 

-.357 

.605 

54 

.27 

-.34 

VO 

.-1 

« 

1 

.424 

5B 

.55 

-.271 

.330 

5C 

-.299 

.455 

5D 

.32 

-.482 

.505 

5E 

.50 

-.52 

-.279 

.493 

5F 

-.48 

-.434 

.691 

5H 

.39 

-.33 

-.439 

.680 

51 

-.253 

.412 

5J 

-.25 

-.31 

-1.44 

-2.31  1.07 

-.410 

•  564 

% 

.42 

.36 

-.294 

.394 

5L 

1.29 

** 

-.382 

•  519 

5M 

*•221 

.544 

.30 

1.27 

-•526 

.613 

6A 

.27 

-.479 

.654 
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Standardized  Canonical  Discriminant  Function  Coefficients 

by  PRC 

(With  POD  Accession  Variables) 


i 

>  * 

i 


1  «r 

i 


4 


Standardized  Canonical  Discriminant  Function  Coefficients 
(With  POD  Accession  Variables) 


DBG 

Month 

AFQT 

NoDip 

HSS 

HSDip 

ColDg 

Age 

White 

Black 

HI  spa 

IF 

-.55 

-.55 

.21 

11 

-56 

-.32 

.17 

.46 

-.79 

-.57 

-.40 

IN 

-1.24 

-.67 

.22 

.26 

.46 

-.39 

3A 

-.73 

-.92 

-.28 

.17 

.39 

3E 

-.77 

-.39 

.27 

-.35 

.40 

-.8? 

-.4? 

31 

-.85 

-.35 

.40 

.48 

-.29 

3L 

-.41 

-.81 

.74 

-.28 

4C 

-1.01 

-.54 

.46 

4D  . 

-.37 

.58 

.49 

46 

-.93 

.45 

4F 

-.74 

.55 

.38 

4K 

-.73 

-.64 

.87 

-.32 

5F 

-.69 

.59 

.22 

.42 

.37 

-1.06 

-1.30 

5H 

-.33 

.71 

.32 

-.44 

5J 

-.31 

.26 

.81 

-.26 

5K 

-.31 

.76 

.39 

-.32 

.45 

.44 

.49 

51 

-.30 

.54 

-.27 

5M 

-.66 

-.48 

5N 

-1.16 

-.66  - 

-.27 

.25 

.59 

-.30 

6a 

-1.10 

.45 

.31 

.67 

.45 

-.39 

6l 

-.56 

.56 

-.59 

o 

• 

i 

ll 
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Standardized  Canonical  Discriminant  Function  Coefficients 
(With  POD  Accession  Variables) 


Marrl  #MHSS  Unerop  ArmyR  DODNP  DODHS  Incom  Non_j^ro^ 


IP 

.29 

.32 

-.31 

.85 

11 

-.32 

-.36 

.50 

-.4? 

.92 

IN 

.29 

-.57 

1.66 

-.41 

.73 

3A 

.17 

.86 

.69 

-.59 

.54 

3B 

1.38 

-.83 

-.41 

.53 

31 

-.43 

-.56 

-.27 

.46 

3L 

-.38 

-1.41 

.30 

4C 

-.77 

.38 

1.29 

-.35 

.54 

4D 

-1.04 

.74 

-.28 

.42 

4E 

.38 

.65 

-.48 

-.53 

.78 

4F 

-1.58 

-1.25 

1.42 

-.44 

.55 

4K 

-.29 

.32 

-.38 

.64 

5P 

-.4? 

-.43 

.68 

5H 

.37 

-.61 

.33 

-.45 

.69 

5J 

-.32 

-1.93 

2.19 

-.41 

.56 

5C 

.26 

.65 

-.80 

-.37 

.49 

5L 

.71 

-.36 

.48 

5M 

.58 

.95 

-.27 

.67 

5H 

1.62 

-.45 

-.53 

.62 

6A 

.36 

-.49 

.67 

6l 

.36 

1.77 

-.47 

1.10 

-1.03 

-.37 

.91 
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Standardized  Canonical  Discriminant  Function  Coefficients 


(Region.  No  Accession  Variables) 


Prior  Sex  Month  AFXJT  NoDip  HSS  HSDip  ColDg  Age  White  Black 


«r 


Standardized.  Canonical  Discriminant  Function  Coefficients 


Ho  Accession  Variables) 


Hi spa  Marri  QMA  #MHSS  Unemp  ArmyR  Incom  BMA 


Centroids 
Non-Army  Art 


All  PS,F 

.07 

-.03 

.12 

.05  -.18 

-.08 

-.36 

.39 

All 

.07 

-.06 

.12 

.05  -.18 

-.06 

-.30 

.41 

Rg  PS,F 

-.06 

-.05 

-.11 

.07 

.06  -.13 

-.36 

.46 

REG  all 

.03 

-.07 

.11 

.05  -.16 

-.28 

.44 

RSV1  PS,F 

.15 

-.10 

.30 

-.24 

-.83 

.39 

RSVl 

.16 

-.14 

-.13 

.31 

-.36 

-.66 

.42 

RSV3  PS,F 

-.12 

.23 

.35 

.14 

-.71 

» 

-1.01 

.70 

RSV3 

-.18 

.22 

— l.o6 

-.71 

.78 

RSV4  PS,F 

-.17 

.44 

.35 

.37 

-.21 

-.69 

.53 

RSV4 

.19 

.26 

.38 

-.54 

.48 

RSV5  PS.F 

-.22 

Cv- 

• 

1 

H 

tH 

• 

1 

.28 

-1.13 

.44 

RSV5 

-.28 

-.53 

.32 

-.93 

.51 

RSVPS.F 

.16 

.22 

.15 

-.50 

-.25 

.35 

.59 

-.80 

.76 

RSV6 

-.18 

.29 

-.71 

-.28 

.26 

.73 

-.51 

.76 

REG1  PS.F 

-.0? 

.10 

.07  -.13 

-.72 

.53 

REG1 

-.08 

.11 

-.09 

-.08 

-.24 

.51 

REG3  PS.F 

-.06 

.12 

-.13 

-.13 

-.16 

-.48 

.47 

RBG3 

-.07 

.14 

.28 

-.39 

.47 

REG4  PS.F 

-.09 

.14 

.18 

.09 

.18 

-.38 

.45 

RBG4 

.13 

.14 

-.32 

-.35 

.47 

RSG5  PS.F 

.07 

-.33 

.41 

RBG5 

.08 

.12 

.09 

H 

• 

1 

• 

1 

-.26 

.37 

RBG6  PS ,F 

-.08 

-.19 

-.36 

.21 

.29 

-.30 

.43 

SBG6 

-.16 

.27 

-.22 

.35 
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Standardized  Canonical  Discriminant  Function  Coefficients 
(Region,  With  POD  Accession  Variables) 


Prior 

Sex 

Month 

AFQ? 

NoDip 

HSS 

HSDip 

ColDg 

REG  all  PS,F 

.34 

.24 

.16 

-.66 

-.22 

-.27 

-.60 

REG  all 

.12 

-.76 

-.18 

-.25 

-.64 

RSV1  PS,F 

.28 

.33 

.21 

-.18 

.38 

-.46 

HSV3  PS,F 

.47 

.63 

-.27 

.46 

RSV5  PS,F 

.40 

.38 

.26 

-.53 

-.69 

REG1 

-.67 

-.45 

.09 

REG3  PS,F 

.33 

.19 

.17 

-.67 

-.30 

-.08 

R3G3 

-•78 

.09 

-.29 

HEG5 

-.72 

-.37 

-.69 

.08 

REG6 

-.73 

.27 

-.24 

Age 

White 

Black 

HI  spa 

Marrl 

QMA 

#MHSS 

Unemp 

REG  all  PS,F 

.35 

-.10 

.13 

-.05 

-.03 

-.15 

REG  all 

.46 

-.08 

.09 

-.14 

RSV1  PS.F 

.45 

-.14 

.12 

-.16 

.32 

RSV3  PS,F 

.37 

.31 

-.12 

.23 

.35 

RSV5  PS,F 

.30 

-.25 

REG1 

.46 

-.29 

.11 

-.21 

REG3  PS,F 

.33 

-.28 

REG3 

.40 

-.21 

-.06 

.24 

BEG5 

.40 

.09 

.13 

BBG6 

.64 

.15 

-.28 

-.16 

Standardized  Canonical  Discriminant  Function  Coefficients 
(Region,  With  DQD  Accession  Variables) 


ArnyR  Incom  BMA  DODNP  DODHS  gentrolds 
_ _ _ _  _ Non-Army  Army 


REG  all  PS.P  .10  -.10  -.15 


REG  all 


.13  -.09 


.16  -.37 

.19  -.29 


RSV1  PS,F 

-.27 

-.fco 

-.84 

.39 

RSV3  PS,P 

.14 

-.71 

.18 

-1.01 

.70 

RSV5  PS,P 

-.13 

-.35 

.32 

-.12 

-1.14 

.1*4 

REG1 

-.45 

.69 

-.25 

•52 

REG3  PS,F 

.36 

-.11 

.21 

-.1*8 

.48 

REG  3 

-.12 

.15 

-.39 

.47 

REG5 

-.11 

-.09 

-.52 

.56 

-.27 

•38 

REG6 

.26 

.25 

-.22 

.36 
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Percentages  of  Cases  Correctly  Classified.  Into  Army  arid  Non-Army 

(No  Accession  Variables) 


Analysis  Cases  Random  Sample 


DRC 

Non-Army 

Army 

Total 

Non-Army 

Army 

Total 

1A 

75.6 

55.6 

69.5 

72.8 

36.1 

62.5 

IB 

65.2 

70.6 

67.5 

69.5 

58.7 

65.O 

1C 

67.4 

62.3 

65.7 

63.9 

54.1 

61.1 

10 

76.2 

56.7 

69.9 

74.6 

43.9 

63.4 

IE 

72.0 

49.3 

64.8 

66.4 

50.0 

61.4 

IF 

72.2 

64.1 

70.1 

70.0 

54.1 

66.2 

1G 

74.8 

74.6 

74.8 

73.3 

56.9 

68.6 

1H 

74.2 

71.4 

73 .5 

71.9 

64.5 

70.1 

11 

81.1 

72.9 

78.3 

65.I 

59.2 

63.2 

1J 

69.8 

51.1 

63.6 

65.2 

51.2 

6I.3 

IK 

72.6 

56.1 

68.3 

76.9 

63.6 

73.0 

1L 

67.6 

65.2 

66.9 

75.4 

58.3 

71.5 

IN 

74.7 

67.9 

72.3 

63.9 

62.5 

63.4 

3A 

68.1 

78.9 

73.7 

66.9 

70.4 

68.5 

\  3B 

80.f0 

73.9 

77.1 

59.4 

73.3 

66.1 

3C 

83.3 

63.3 

75.3 

62.2 

64.5 

63.1 

3D 

69.5 

72.6 

71.0 

60.4 

72.7 

66.3 

3E 

68.3 

66.4 

67.5 

65.5 

59.0 

62.6 

3F 

73.0 

63.3 

68.3 

66.2 

61.7 

64.1 

30 

77.0 

64.9 

72.6 

61.7 

57.4 

60.3 

3H 

69.9 

67.6 

69.0 

79.4 

72.0 

76.0 

31 

65.3 

62.1 

64.2 

60.2 

61.0 

60.6 

3J 

68.2 

67.2 

67.7 

63.8 

64.2 

64.0 
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] 

1 

I 


h 


ft 


l 


M 


h 


M 


Analysis  Cases  Random  Sample 


DRC 

Non- Army 

Army 

Total 

Non-Army 

Army 

Total 

3* 

62.1 

67.1 

64.3 

69.1 

65.8 

67.6 

3L 

100.0 

78.7 

82.5 

84.2 

69.4 

72.8 

4a 

65.4 

72.5 

68.5 

46.2 

68.1 

55.4 

4C 

64.2 

69.8 

66.4 

69.3 

42.0 

58.4 

4D 

69.0 

59.6 

65.3 

67.0 

69.2 

68.0 

4E 

80.6 

78.6 

79.8 

67.5 

61.8 

65.3 

4P 

62.3 

73.5 

67.3 

73.0 

61.3 

67.2 

40 

65.6 

65.6 

65.6 

61.2 

62.5 

61.8 

4H 

73.4 

81.1 

76.9 

66.7 

65.5 

66.1 

41 

69.4 

68.4 

69.1 

63.O 

56.0 

60.6 

4J 

69.2 

52.8 

6I.3 

6  5.7 

51.5 

61.0 

4K 

75.4 

72.2 

74.2 

76.1 

57.4 

68.6 

5A 

72.4 

56.2 

65.5 

62.4 

64.3 

63.3 

5B 

61.1 

59.5 

60.4 

63.0 

62.3 

62.7 

5C 

87.6 

41.0 

69.2 

88.2 

43.9 

71.6 

5D 

66.3 

74.1 

70.1 

63.O 

72.2 

67.O 

5E 

71.7 

55.9 

66.0 

70.0 

45.5 

60.5 

5P 

74.3 

71.8 

73.4 

74.6 

45.9 

63.3 

5H 

78.7 

64.6 

73.2 

79.5 

52.1 

67.1 

51 

63.2 

55.6 

60.3 

68.6 

60.6 

65.6 

5J 

74.0 

67.9 

71.4 

65.5 

48.3 

58.5 

5K 

65.9 

61.8 

64.2 

64.8 

51.4 

58.8 

5L 

62.5 

64.2 

63.2 

51.3 

63.6 

56.3 

5M 

63.4 

60.0 

62.4 

60.0 

61.9 

60.7 

5N 

71.7 

67.0 

69.5 

68.9 

59.3 

64.9 

6A 

75.0 

70.7 

73.2 

64.9 

43.6 

57.8 
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Total 


Total 


Analysis  Gases 
Non-Army  Army 


Random  Sample 
Non-Army  Army 


6e 

63.0 

51.9 

57.4 

59.4 

63.2 

61.4 

6F 

68.8 

70.0 

69.2 

53.1 

65.I 

57.5 

6g 

76.5 

64.8 

71.3 

67.1 

34.0 

53.5 

6h 

57.4 

60.0 

58.2 

52.6 

52.4 

52.5 

6l 

72.7 

63.5 

69.8 

64.9 

54.9 

60.7 

6J 

94.3 

50.0 

83.1 

91.5 

34.8 

77.7 

6k 

74.2 

62.7 

70.1 

69.0 

56.4 

64.5 

6l 

70.2 

67.7 

75.9 

58.7 

49.0 

54.8 

1 

Percentages  of  Cases  Correctly  Classified  into  Army  and  Non-Army 
(With  DOD  Accession  Variables) 

9 

Analysis  Case3  Random  Sample 

DRC  Non-Army  Army  Total  Non-Army  Army  Total 


IP 

74.1 

71.8 

73.5 

67.5 

54.1 

64.3 

s 

11 

78*9 

75.0 

77.6 

61.3 

55.1 

59.4 

IN 

74.7 

71.4 

73.6 

63.9 

62.5 

63.4 

in 

3A 

69.9 

78.9 

?4.6 

64.6 

67.6 

66.0 

P 

3E 

70.0 

67.9 

69.1 

66.1 

57.6 

62.3 

► 

L- 

31 

68.3 

58.6 

64.8 

54.2 

59.8 

57.0 

r- 

L 

3L 

100.0 

76.6 

80.7 

78.9 

74.2 

75.3 

4C 

64.2 

69.8 

66.4 

68.0 

42.0 

57.6 

4D 

70.2 

56.1 

64.5 

67.0 

64.6 

66.0 

4E 

79.0 

81.0 

79.8 

70.0 

63.4 

67.8 

i  C' 

4P 

62  f  3 

71.4 

66.4 

74.6 

53.2 

64.0 

4K 

68.9 

69.4 

69.1 

74.6 

51.1 

65.3 

!■ 

a 

5P 

75.2 

74.6 

75.0 

74.6 

47.3 

63.8 

5H 

78.7 

64.6 

73.2 

76.9 

53.1 

66.2 

* 

5J 

74.0 
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HSV1  PS,F  72.8  73.7  73.4 

RSV3  PS,F  90.4  74.1  80.8 
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i 


I  F8  P 

1  .812  81.  HI 

2  .*2S  48.588 

3  .837  27. Ml 

♦  .849  28.258 

S  .882  18.288 

8  .874  13.588 

7  .888  11.571 

8  .899  18.125 

9  .111  9.888 

18  .123  8.188 

U  .138  7.384 

12  .148  8.758 

13  .188  8.231 

14  .173  5.788 

>3  .185  5.488 

18  .198  5.883 

17  .218  4.785 

18  .222  4.588 

19  .235  4.283 

28  .247  4.858 

21  .259  3.857 

22  .272  3.882 

23  .284  3.522 

24  .298  3.375 

25  .389  3.248 

28  .321  3.115 

27  .333  3.888 

28  .348  2.893 

29  .358  2.793 

38  .378  2.788 

31  .383  2.813 

32  .395  ‘  2.531 

33  .487  2,455 

34  .428  2.382 

35  .432  2.314 

34  .444  2.258 

37  .457  2,189 

38  .489  2.132 

39  .481  2.877 

.494  2.825 

AVERME  INTENSITY* 


1NTNCTT 
83278274.377 
5235338.593 
1719818.278 
3579549.389 
1719448.891 
2411587.427 
491489.172 
1818288.182 
338492.184 
584928.887 
439578.592 

139719.883 
2423898.524 

588393.978 

451279.883 
1371888.843 
1522421.178 
879873.988 
1218188.884 

118787.781 
8S231.885 
881248.958 
788999.175 
283138.414 
449778.928 
148898.881 
248841.291 
1714288.933 
354857.138 
558887.917 
182281.432 
28143.982 
288172.371 
89922.787 
272449.232 
488888.585 
298358.198 
33587.349 
95789.782 
385335.894 
2421741.8245555488 
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i  Fa  p 

1  >#13  aa.Mt 

2  .«25  48.811 

3  .#38  24.647 

4  .«SI  21. N# 

5  .143  16. N# 

6  .175  13.333 

7  .MS  11.429 

8  .IN  18.N8 

9  .113  8.889 

II  .125  8.HI 

11  .138  7.273 

12  .151  6.447 

13  .163  6.154 

14  .175  5.714 

15  .188  5.333 

16  .2N  5.8N 

17  .213  4.786 

18  .225  4.444 

19  .238  4.211 

21  .258  4.8N 

21  .263  3.818 

22  .275  3.436 

23  .2N  3.478 

24  .3N  3.333 

25  .313  3.2N 

24  .325  3.177 

27  .338  2.963 

28  .358  2.857 

29  .363  2.759 

38  .375  2.647 

31  .388  2.581 

32  .4N  2.588 

33  .413  2.424 

34  .425  2.353 

35  .438  2.286 

36  .458  2.222 

37  .463  2.162 

38  .475  2.185 

39  .488  2.851 

48  .5N  2.N8 

average  intensity 


INTNCTT 
425263.264 
95272.234 
142867.746 
438195.258 
389961.518 
554123.886 
117826.221 
221731.274 
276989.756 
768311.686 
368433.641 
15813.418 
2414888.482 
148289.517 
316635.128 
2424474.818 
1781752.256 
1355883.652 
1184398.989 

285641.225 

284164.837 
2536959.417 
1258784.786 
669837.154 
1846371.276 

769922.838 
278799.387 

5259793.421 
2436559.857 
1521658.649 
171317.897 
267971.656 
1334378.454 
367933.398 
1499523.795 
748633.271 
678635.832 
626438.953 
979614.676 

238584.225 
918483.2147895396 
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Appendix  D-5 


AGG13AH  -  Bi'tDP  Output  for  ARI  (2,1)  Model 


MIMA 


VARIABLE  IS  ACC13AM. 
Bf ORDER  IS  1. 

ARORDERS  ARE  Ml. 2)'. 
CENTERED./ 


ESTIMATION  RESIDUAL  >  RACC ./ 


ESTIMATION  BY  CONDITIONAL  LEAST  SQUARES  METHOD 

RELATIVE  CHANCE  IN  EACH  ESTIMATE  LESS  THAN  .1888E-83 

VARIABLE  VM.  TYPE  MEAN  TIME  DIFFERENCES 

I 

ACC 13AM  RANDOM  REMOVED  1-  81  Il-B  I 

PARAMETER  VMIASLE  TYPE  FACTOR  ORDER  ESTIMATE  ST.  ERR.  T-RATIO 

1  ACC 13AM  M  1  I  -.2814  .1892  -2.57 

2  ACC  13AM  M  1  ,  2  -.278,4  .1117  -2.74 

RESIDUAL  SUN  OF  SQUARES  *  27219117.115255 
DECREES  OF  FREEDOM  *  74 

RESIDUAL  MEM  SOUARE  *  358144.277832 


ESTIMATION  BY  BACKCASTINC  METHOD 

RELATIVE  CHANCE  IN  RE61DUAL  SUM  OF  SOUARES  LESS  THAN  .1888E-84 

VARIABLE  VM.  TYPE  MEAN  TIME  DIFFERENCES 

l 

ACC13AH  RANDOM  REMOVES  1-  81  It'D  I 


PARAMETER  VMIABLE  TYPE  FACTOR  ORDER  ESTIMATE  ST.  ERR.  T-RATIO 

1  ACC 13AM  M  1  I  -.2743  .1938  -2.95 

2  ACC13AH  M  1  2  -.3249  .8937  -3.49 

RESIDUAL  SUM  OF  SOUARES  *  27385853.812457  (BACKCASTS  EXCLUDED! 
DEGREES  OF  FREEDOM  *  74 

RESIDUAL  MEM  SOUARE  *  359287.558144 
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VARIABLE  IS  RACC.  HAILAC  IS  24./ 


ACF 


NUHBER  Of  OBSERVATIONS  >  81 
MEAN  OF  THE  (DIFFERENCEB)  SERIES  *  -11.4824 
STANDARD  ERROR  OF  THE  MEAN  49.819# 
T-VALUE  OF  MEAN  (AGAINST  ZERO)  >  -.1445 


AUTOCORRELATIONS 

1-  12  -.#3  -.#3  -.#9  -.#5  .13  .14  -.#7  -.#4  .14  -.#7  .13  .12 

ST.E.  .11  .11  .11  .11  .11  .11  .12  .12  .12  .12  .12  .12 

13-  24  #.f  .22  -.17  -.14  .#2  .11  .#3  .85  -.1#  #.#  .14  .87 

ST.E.  .12  .12  .12  .13  .13  .13  .13  .13  .13  .13  .13  .13 


PLOT  OF  SERIAL  CORRELATION 

-1.8  -.8  -.4  -.4  -.2  .8  .2  .4  .4  .8  1.8 


LAC 

CORN.  4~ 

1 

1 

-.832 

4 

11  4 

2 

-.829 

* 

11  4 

3 

*.895 

♦ 

111  4 

4 

-.847 

♦ 

11  * 

5 

.125 

4 

1111  ♦ 

4 

.139 

4 

111!  4 

7 

-.849 

4 

111  4 

8 

-.843 

4 

II  * 

9 

.845 

4 

III  * 

18 

-.847 

4 

III  4 

11 

.133 

4 

till  4 

12 

.822 

4 

II  4 

13 

-.888 

4 

1  ♦ 

14 

.214 

4 

1 1 1 1 1 1 4 

15 

-.871 

4 

III  4 

14 

-.837 

4 

11  * 

17 

,.821 

4 

11  4 

19 

.US 

4 

IIII  * 

19 

.829 

4 

II  4 

28 

.852 

4 

II  4 

21 

-.184 

4 

IIII  4 

22 

-.884 

4 

1  4 

23 

.837 

4 

11  4 

24 

.874 

4 

III  <4 
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VARIADLE  IS  RACC.  HAILAC  IS  24./ 


PACF 


NUHSER  Of  OISCRVATIOMS  *  81 
MEAN  Of  THE  (DIFFERENCED)  SERIES  *  -11.4824 
STANDARD  ERROR  OF  THE  KEAN  49.6191 
T-VALUE  OF  (CAN  (AGAINST  ZERO)  *  -.1445 


PARTIAL  AUTOCORRELATIONS 

1-  12  -.«3  -.03  -.11  -.84  .12  .14  -.14  -.12  .18  -.86  .89  .84 

ST.E.  .11  .11  .11  .11  .11  .it  .11  .11  .11  .U  .It  .11 

13-  24  .83  .23  -.85  -.83  .  83  .13  -.81  -.81  -.81  .81  -.84  .87 

ST.E.  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11 


PLOT  OF  SERIAL  CORRELATION 


|Af 

-1.8  -.8  -.4 

-.4  .2 

.8  .2 

Ufv 

1 

11  ♦ 

i 

-.832 

♦ 

2 

.-.838 

* 

11  ♦ 

3 

-.897 

■* 

III  * 

4 

-.855 

■* 

11  ♦ 

5 

.117 

♦ 

1111  4 

4 

.139 

♦ 

1111  ♦ 

7 

-.843 

♦ 

III  ♦ 

8 

-.824 

♦ 

II  ♦ 

9 

.181 

♦ 

mi  ♦ 

18 

-.878 

♦ 

in  ♦ 

11 

.887 

♦ 

nx  ♦ 

12 

.848 

.*■  n  ♦ 

13 

.824 

♦ 

n  ♦ 

14 

.232 

♦ 

niu+i 

15 

-.854 

♦ 

u  ♦ 

14 

-.832 

♦ 

ii  ♦ 

17 

.828 

n  ♦ 

18 

.128 

♦ 

mi  * 

19 

-.889 

♦ 

i  ♦ 

28 

-.887 

♦ 

i  ♦ 

21 

-.812 

♦ 

i  ♦ 

22 

.889 

-.848  * 

♦ 

i  ♦ 

23 

♦ 

ii  ♦ 

24 

.871 

♦ 

in  ♦ 
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MUM 


VARIABLE  IS  ACC13AM. 
DFORBER  IS  t. 

ARORDERS  ARE  '(1.2tl4l\ 
CENTERED./ 


ESTIMATION  RESIDUA.  «  RACC  ./ 

EST1MT10N  RT  CONDITIONAL  LEAST  SQUARES  METHOD 
RELATIVE  CHANCE  IN  EACH  ESTIMATE  LESS  THAN  .1989E-93 


VARIABLE  VAR.  TIRE  MEAN  TINE  DIFFERENCES 

I 

ACC  13AM  RANSOM  REMOVED  I*  SI  U-l  ) 


PARAMETER  VARIABLE 

TTPE 

FACTOR 

ORDER  ESTIMTE 

ST.  ERR. 

T-RATIO 

1  ACC  13AM 

AR 

1 

*  1 

-.359* 

.0983 

-3.45 

Z  ACC 13AM 

AR 

1 

2 

-.2042 

.9945 

-2.94 

3  ACC 13AM 

AR 

1 

14 

.1031 

.185* 

2.15 

RESIDUAL  SUM  OF  SQUARES  «  USD 1944.398 23 4 
DECREES  OF  FREEDOM  *  43 

RESIDUAL  MEAN  SQUARE  «  234191.498305 


ESTIMATION  Bf  BACKCASTINC  METHOD 


VARIABLE  VAR.  TTPE  MEAN  TIME  DIFFERENCES 

l 


ACC13AN 

RANDOM 

REMOVED  1- 

81  ll-B  I 

PARAMETER  VARIABLE 

TTPE 

FACTOR 

ORDER  ESTIMTE 

ST.  ERR. 

T-RATIO 

1 

ACC13AH 

AR 

1 

1  -.1835 

.9893 

-2.29 

•  2 

ACC 13AM 

AR 

l 

2  -.2429 

.9784 

-3.99 

3 

ACC13AN 

AR 

1 

14  .2479 

.9724 

3.48 

RESIDUAL  SUM  OF  SQUARES  *  15743199.443337  (BACKCASTS  EXCLUDED! 
DECREES  OF  FREEDOM  *  43 

RESIDUAL  MEAN  SQUARE  *  2592*7.927989 
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ACf 


VARIABLE  IS  MCC.  MILM  IS  24./ 


HUNKS  OF  OBSERVATIONS  *  01 
MEAN  OF  THE  (DIFFERENCES)  SERIES  *  *11.1227 
STANDARD  ERROR  Of  THE  MEAN  *  40.1453 
T- VALUE  OF  HE AN  (ACAINST  ZERO)  «  -.1432 


AUTOCORRELATIONS 

1-  12  -.05  -.03  -.10  -.04  .10  .19  -.12  *.09  .14  -.07  .13  -.03 

ST.E.  .11  .11  .11  .11  .11  .11  .12  .12  .12  .12  .12  .12 

13-  24  -.11  -.02  .13  .01  .13  .11  .04  .02  -.09  .02  0.0  .09 

ST.E.  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12 


PLOT  OF  SERIAL  CORRELATION 

-1.0  -.0  -^4  -.4  -.2  .0  .2  .4  .4  .0  1.0 


LAC 

CORR.  4~ 

1 

1 

-.048 

* 

11  4 

2 

-.027 

♦ 

II  4 

3 

-.105 

♦  1111  ♦ 

4 

-.040 

♦ 

II  ♦ 

S 

.104 

4 

III!  4 

4 

.193 

4 

1X1X114 

7 

-.022 

4 

11  4 

0 

-.092 

4 

III  4 

9 

.057 

4 

II  4 

10 

-.073 

4 

III  4 

11 

.132 

4 

IIII  4 

12 

-.025 

,  ♦ 

II  4 

13 

-.004 

4 

I  4 

14 

-.020 

4 

II  4 

15 

.030 

4 

II  4 

14 

.009 

4 

I  4 

17 

.025 

4 

II  4 

10 

.112 

4 

IIII  4 

19 

.034 

4 

II  4 

20 

.023 

4 

II  4 

21 

-.084 

4 

III  4 

22 

.019 

4 

I  4 

23 

.002 

4 

1  4 

24 

.007 

4 

III  4 
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PACF  VARIABLE  IS  RACC.  HAILAC  IS  24./ 


MUN9ER  OF  OtSESVATIOMS  »  01 
MEAN  OF  THE  (BIFFERENCEOI  SERIES  *  '11.1227 
STANMRO  ERROR  OF  THE  MEAN  *  40.1453 
T-VALUE  OF  MEAN  (AGAINST  ZERO!  *  -.1432 


PARTIAL  AUTOCORRELATIONS 

1-  12  -.95  -.93  -.11  -.95  .19  .19  9.9  -.1 

ST.E.  .11  .11  .11  .11  .11  .11  .11  .11 

13-  24  .91  .91  .92  .91  9.9  .13  .  90  .1 

ST.E.  .11  .11  .11  .11  .11  .11  .11  .11 


PLOT  OF  SERIAL  CORRELATION 


LAC 

-1.9  -.0 

rnoo  : 

-.4  -.4  - 

.2  .9 

Uv 

l*UNK«  »  ▼ 

I 

l 

-.940 

♦ 

II  * 

2 

-.929 

♦ 

II  * 

3 

-.190 

+  mi  ♦ 

4 

-.953 

* 

II  ♦ 

5 

.994 

♦ 

III  ♦ 

4 

.195 

♦ 

min 

7 

-.992 

♦ 

i  ♦ 

0 

-.971 

♦ 

in  ♦ 

9 

.997 

•*  ♦ 

iii  ♦ 

(9 

-.947 

♦ 

in  ♦ 

11 

.975 

♦ 

in  ♦ 

12 

-.945 

♦ 

a  ♦ 

13 

.990 

♦ 

i  ♦ 

14 

.997 

t 

i  ♦ 

15 

.91S 

♦ 

i  ♦ 

14 

.913 

♦ 

i  ♦ 

17 

.992 

♦ 

i  ♦ 

10 

.129 

♦ 

IIKI  * 

19 

.995 

♦ 

ill  *' 

29 

.919 

♦ 

I  ♦ 

21 

-.959 

♦ 

II  ♦ 

22 

.911 

♦ 

I  ♦ 

23 

-.997 

♦ 

I  ♦ 

24 

.917 

♦ 

I  ♦ 

.19  -.97  .97  -.95 
.11  .11  .11  .11 

-.94  .91  -.91  .92 
.11  .11  .11  .11 


.4  .0  1.9 

♦ - ♦ - ♦ 
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BMDP  Output  for  Optimal  Unemployment  Rate  Model  and 
a  Comparison  of  Cumulative  Perlodogram  Plots 


MUM  VMIAKE  |8  UNEMP. 

ARMIERS  ME  Ml»2>». 
ARVALUES  •  1.13»i  -.1477. 
CCMTEftCO./ 


THE  COMPONENT  MS  KEN  ABBED  TO  THE  MODEL 

* 

THE  CUMEMT  NOBEL  HAS 
OUTPUT  VARIABLE  >  UNEMP 
INPUT  VM1ABLE  «  NOISE 


ESTIMTION  RESIDUALS  «  Rl. 

NETHOI  IS  CIS./ 


ESTIMTION  IT  COMITIONAL  LEAST  SOUARES  METHOD 
RELATIVE  CHANCE  IN  EACN  ESTIMATE  LESS  THAN  .1M0E-I3 


SUNMNT  OF  THE  MOREL 

OUTPUT  VARIABLE  -  UNEMP 
INPUT  VARIABLES  -  NOISE 

VARIABLE  VAR.  TYPE  MEAN  TIME  DIFFERENCES 
UNEMP  RANDOM  REMOVED  1-  SI 


PARAMETER  VARIABLE  TYPE  FACTOR  ORDER  ESTIMTE  ST.  ERR.  T-RAT1D 

1  UNEMP  M  1  1  1.1301  .1127  10.03 

2  UNEMP  M  1  2  -.1477  .1154  -1.21 

RESIDUAL  SUM  OF  SQUARES  «  5.045450 

KCREES  OF  FREEDOM  '  77 

RESIMML  KAN  SOUARE  >  .045742 


VARIAHE  IS  RI.  MM  LAC  IS  29./ 


j 


\ 


ACF 


MIRIER  OF  OBSERVATIONS  * 

74 

REAR  Of  THE  (DIFFERENCES!  SERIES  « 

.9191 

STANDARD  ERROR  OF  THE  HEAR 

.9284 

T-VALUE  OF  REAM  (ACAINST  ZERO!  * 

.3527 

AUTOCORRELATIONS 

-.12  .12  9.9  .14  -.12  .12  .13  -.12  .11  -.11  .If  .14 

.11  .11  .11  .12  .12  .12  .12  .12  .12  .12  .12 


1- 

ST.E 


12 


t. 


13-  21  .94  -.14  -.15  .94  .14  -.91  .98  -.91 

ST.E.  .12  .12  .12  .12  .12  .13  .13  .13 


PLOT  9f  SERIAL  CORRELATION 

-l.f  -.9  -.4  -.4  -.2  .9  .2  .4  .4  .9  1.9 


LAC 

CORN.  ♦ - ♦ - ♦- 

I 

1 

-.924 

*  11 

2 

.129 

4  IIII 

3 

.993 

♦  I 

4 

.141 

4  mu 

5 

-.923 

4  11 

4 

.921 

4  11 

7 

.939 

4  II 

8 

-.923 

4  11 

4 

.974 

4  III 

19 

-.lit 

4  IIII 

11 

.984 

4  III 

12 

.941 

4  111 

13 

.941 

4  HI 

14 

-.141 

4  mu 

15 

-.944 

4  II 

14 

.934 

4  II 

17 

.184 

4  iiiiii 

19 

-.919 

4  1 

14 

.977 

4  III 

21 

-.994 

4  1 

VARIAKE  IS  RI.  MILAS  IS  21./ 


PACF 


NUH8ER  OF  OBSERVATIONS  >  •  79 
(CAN  OF  THE  (B.IFFERENCEB)  SERIES;  *  .fill 
STANDARD  ERROR  OF  THE  MEAN  *  .8286 
T-VALUE  OF  DEAN  (ASA INST  ZERO)  *  .3527 


PARTIAL  AUTOCORRELATIONS 

1-  12  -.82  .12  .81  .15  -.82  -.81  .83  -.85  .80  -.11  .86  .18 

ST.E.  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11  .11 

13-  28  .84  -.13  -.18  .84  .22  .82  .86  -.86 

ST.E.  .11  .11  .11  .11  .11  .11  .11  .11 


PLOT  OF  SERIAL  CORRELATION 

-1.8  -.8  -.6  -.4  -.2  .8  .2  .4  .6  .8  1.8 


LAS 

CQRR. 

1 

1 

-.824 

..  .  ♦ 

II  * 

2 

.119 

♦ 

11(1  « 

3 

.889 

♦ 

1  ♦ 

4 

.149 

* 

mu  * 

5 

-.818 

♦ 

i  ♦ 

6 

-.814 

♦ 

i  * 

.833 

♦ 

II  4 

0 

-.848  < 

♦ 

II  4 

4 

.888 

♦ 

IK  4 

18 

-.188 

♦ 

HU  4 

11 

.862 

* 

III  4 

12 

.181 

♦ 

IIII  4 

13 

.856 

♦ 

II  4 

14 

-.126 

♦ 

IIII  4 

15 

-.183 

♦ 

(III  4 

16 

.842 

♦ 

II  4 

17 

.223 

♦ 

mini 

18 

.821 

* 

II  4 

19 

.86# 

♦ 

III  4 

21 

-.864 

♦ 

HI  4 
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UNEMP  RATE 
ARIMA  (2.  0,  0) 

t 

-  5%  KS  BAND 

-  25%  KS  BAND 


AD-A124  84S  FACTORS  INFLUENCING  ARHV  ACCESSIONS(U)  AIR  FORCE  INST  4/4 
OF  TECH  URIGHT-PATTERSON  RFB  OH  SCHOOL  OF  ENGINEERING 
K  H  KALINICH  ET  AL.  DEC  82  AFIT/G0R/0S/82D-? 

F/G  12/1  NL 


UNCLASSIFIED 


MINA 


VMIMLE  IS  ACC13M. 

ARORKRS  ME  *11*41  •  • 

ARVALUES  «  .9200i  .4409. 

CEMTEREI./ 

TME  COMPONENT  MIS  KEN  Mil  TO  THE  NOKL 

TIC  CURRENT  NOKL  MM 
OUTPUT  VMIMLE  •  ACC13AH . 

INPUT  VMIMLE  ■  HOIK 

INKP  VMIMLE  IS  UNEHP. 

UPORKRS  *  *  (I  i7*8)  * . 

UPVALUES  >  743.»i  -441.0.  379.1./ 

THE  COMPONENT  MA6  KEN  ANEB  TO  THE  MODEL 

THE  CUMENT  NOKL  HAS 
OUTPUT  VMIAILE  ■  ACC13AH 
INPUT  VMIAILE  •  NOISE  INEMP 

ESTIMATION  RiSIBUALS  IS  RVI. 

NETHOI  IS  CLS./ 

ESTIMATION  If  COMIITIONAL  LEAST  S0UANE6  NETHOI 

NELATIVC  CHANCE  IN  KSIOUAL  SUN  OF  SOIMKS  LESS  THAN  .100K-04 

SUNNMT  OF  THE  NOKL 

OUTPUT  VMIAILE  -  ACCI3AN 
INPUT  VMIMLES  -  NOISE  UNEHP 

VMIAILE  VM.  TTPE  NEM  TIME  DIFFERENCES 

ACC13M  RAMON  REMOVE!  I-  II 


imp  UHIM 

1- 

It 

KRAKTER  VMIAILE 

TTPE 

FACTOR 

ORDER  ESTIMATE 

ST.  ERR. 

T-RATIO 

1  ACC19M 

M 

1 

I  .S2M 

.0007 

9.94 

2  ACCUM 

M 

1 

4  .4419 

.0077 

9.29 

3  INEMP 

UP 

1 

1  712.9327 

143.9109 

4.00 

4  UNEMP 

UP 

1 

7  -441.4121 

231.3032 

-1.91 

9  UNEMP 

UP 

1  370.0971 

231.3290 

1.44 

KSIOUAL  SUN  OF  SOUMES  *  14401993.01999* 
KOREES  OF  FREEMN  ■  42 

AESIOUAL  NEM  MUAK  •  247744.904129 
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VARIAOLE  16  MI.  NAILAC  IS  21./ 


ACF 


mm  OF  OOSERVATIONS  >  47 
KM  OF  INC  (DIFFERENCED)  SCRIES  >  -0.4231 
STMBMI ERAM  OF  TK  KM  >  41.2432 
T-9ALK  OF  KM  (AGAINST  ZERO)  »  -.I4M 


AUTOCORRELATIONS 

1-  12  -.17  .04  .11  -.01  .10  -.13  -.05  .07  -.02  -.03  .11  -.09 

ST.E.  .12  .13  .13  .13  .13  .13  .13  .13  .13  .13  .13  .13 

13-  20  .03  -.11  .07  -.12  -.09  .01  -.09  -.04 

ST.E.  .13  .13  .14  .14  .14  .14  .14  .14 


PLOT  OF  SERIAL  CORRELATION 


VARIAKE  IS  Rll.  AAILAC  IS  20./ 


HCf 


NUMER  Of  QISERVATIQNS  *  47 
REM  OF  THE  (BIFFERENCE9)  SERIES  *  -8.4Z3I 
STMMRS  ERROR  OF  Ttt  REM  *  41 .2432 
T-VN.UE  Of  HEM  (RM1MST  ZERO!  •  -.1403 


PARTIAL  AUTOCORRELATIONS 

■  I*  12  -.17  .13  .12  .13  .19  -.12  -.11  .14  .04  -.12  .12  -.04 

ST.E.  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12 

13-  20  -.03  -.13  .17  -.11  -.07  -.12  -.17  -.09 

ST.E.  .12  .12  .12  .12  .12  .12  .12  .12 


PLOT  Of  SERIAL  CORRELATION 


lie 

-1.0  -.0 
MM  4,-1  ,, 

-.4  -.4  -.2  .0 

UHr 

0 

Ml.  1 - V" 

I 

1 

-.171 

♦  1(111 

2  ' 

.029 

♦  II 

3 

.129 

♦  IIII 

4 

.029 

♦  11 

9 

.093 

♦  III 

4 

-.119 

♦  IIII 

7 

-.114 

♦  IIII 

0 

.037 

♦  II 

9 

.037 

*  II 

10 

-.017 

♦  1 

11 

.120 

♦  ’  IIII 

12 

-.049 

♦  III 

13 

-.134 

♦  II 

14 

-.130 

♦  IIII 

19 

.040 

'  ♦  III 

-.113 

*  IIII 

17 

-.049 

*  III 

10 

-.014 

♦  I 

19 

-.049 

♦  III 

2# 

-.007 

♦  III 
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MINA 


VMIAOLE  1$  UHCMP. 

MfiUlltS  ABC  *11.21*. 
ABVAUICS  *  1.13#.  '.1477. 
CENTERS!./ 

THE  COMPONENT  HAS  SEEN  AMES  TO  THE  HODCL 

THE  CUMENT  RQOEL  HAS 
OUTPUT  VMIAOLE  •  UNEHP 
INPUT  VMIAOLE  •  NOISE 


FORECAST  CASES  ABE  24.  JOIN./ 


FORECAST  OH  VMIAOLE  UNEHP  FROM  TINE  PEN  1  OB  02 


PER 101 

FORECASTS 

ST.  ERR 

02 

4.43750 

.25*54 

03 

1.40447 

.30704 

04 

4.34140 

.40342 

H 

4.31434 

.5*041 

I* 

4.2*01* 

.42*33 

07 

4.22244 

.40320 

N 

4.170*7 

.73354 

Of 

4.13531 

.77005 

ff 

4.0420* 

.01441 

fl 

4.05130 

.05742 

42 

4.010*0 

.04141 

43 

0.47075 

.4237* 

44 

0.43173 

.45331 

45 

0.04352 

.4H74 

4* 

0.05*11 

1.00*44 

47 

0.01440 

1.03044 

40 

0.703*1 

1 .05243 

44  . 

0.74044 

1.07405 

IN 

0.71411 

1.04341 

101 

0.40043 

1.112*2 

112 

0.44747 

1.13027 

IN 

0.41510 

1.14*44 

IN 

0.50350 

1.1*270 

IN 

0.552*3 

1.177*0 

STANOAU  ERROR  «  .25*444  III  CONOITIONAL  METHOD  I 

PSIUE1CNT  NA1PSI  >  40./ 

41  PSl-NEICHTS  ARE  STORE! 


UNIVARIATE  TINE  SCRIES  MEL  ERASES 


ARINA  VARIAILE  IS  ACC13AH. 

ARMBERS  ARE  MliAI*. 
ARVALUES  *  .5218,  .4MS. 
CENTEIEI,/ 


THE  CQNfQNENT  HAS'BEEH  ASHER  TO  THE  HOOEL 

THE  CURRENT  NOBEL  HAS 
OUTPUT  VARIAILE  •  ACCI3AN 
INPUT  VARIAILE  •  WISE 


tNBCP  VARIAILE  IS  UNENP. 

WORSENS  *  1  <1 17*8>  * . 

UPVALUES  •  7*3.1,  37M. 

CENTEREI./ 


THE  CONPOHEHT  HAS  BEEN  ASKS  TO  THE  MODEL 

THE  CURRENT  NOIEL  NAS 
OUTPUT  VARIAILE  >  ACC13AN 
INPUT  VARIAILE  «  WISE  UNENP 


FORECAST  CASES  '  41/ 


F0RECA6T  CASES  «  4#./ 

FORECAST  ON  VNUMLE  ACC13AM  FROM  TINE  PERIOO  02 


PER 101 

FORECASTS 

ST.  ERR 

02 

4777.45045 

519.85844 

03 

5844.97587 

424.29125 

14 

4948.743N 

408.43980 

15 

4882. 42211 

735.11545 

84 

4494.10734 

774.89383 

87 

4731.33875 

018.54484 

88 

4725.19448 

079.44444 

89 

4495.89901 

921.17142 

98 

4754.05254 

941.44798 

91 

4490.90499 

991.404N 

92 

4577.10283 

1814.83174 

93 

4575.24554 

1834.27479 

94 

4554.00842 

1841. 18N2 

95 

4527.49499 

1894.43149 

94 

4539.52844 

U24.441N 

97 

4521.05505 

1152.87122 

98 

4457.83394 

1173.18457 

99 

4424.57848 

1189.19910 

IN 

4399.92849 

1284.44709 

111 

4374.74348 

1224.02741 

182 

4348.31424 

1240.88142 

IN 

4358.18739 

1249.77744 

104 

4324.52374 

1207.45433 

IN 

4292.7I4N 

1382.14350 

IN 

4245.79387 

1315.79885 

MSO  ON  THE  AVAILABLE  MTRi  OH!  25  FORECASTS  CAN  K  MADE 
STANBARI  ERROR  •  519.051  (Of  CONDITIONAL  METHOD  > 


Appendix  E-6 

Enhanced  Transfer  Function  and  Forecasts 


MUM 


VMIMLE  18  ACCI3M. 

AMRIERS  ME  Ml.41'. 
ARVALUES  *  .5200*  .4415. 
IMOMERS  ME  •III*. 
CENTERED./ 

TIC  CONPONENT  MS  KEN  ADBEB  TO  THE  KOOEL 

THE  CURRENT  NOBEL  MS 
OUTPUT  VMIMLE  «  ACCOM 
INPUT  VMIMLE  ■  NOISE 


> 

1N0EP  ,  VMIMLE  IS  UNEHP. 

UPORBERS  •  Mli7t8)'. 

UPVALUES  «  57l.li  -S18.0i  434.0./ 

TIC  CONPONENT  MS  KEN  ADKB  TO  THE  HOKL 

THE  CURRENT  NOKL  HAS  ‘ 

OUTPUT  VMIAOLE  •  ACCOM 
INPUT  VMIABLE  >  HOIK  UNEHP 

ESTIMT1QN  RESIDUALS  IS  RU. 

NETHOO  IS  CIS./ 


fj'  EST1MT10N  II  COMITIQML  LEAST  S8UARES  HETHOD 

RELATIVE  CHANCE  IN  RESIIUAL  SIM  OF  SHIRKS  LESS  THAN  .1MK-I4 

SUNMRf  OF  THE  NOKL 

OUTPUT  VMIABLE  --  ACC13AH 
INPUT  VMIABLES  -  HOIK  UNEHP 


VMIAOLE  VM.  TTPE  KEAN  TIK  DIFFERENCES 
ACCOM  RAMON  REMOVES  1'  01 


lUgrUP  ftAUtUYM 

1* 

01 

MAKTEI  VMIAOLE 

TTPE 

FACTOR 

ORDER  ESTIMTE 

ST.  ERR. 

T-RATIO 

I  ACC13M 

M 

1 

1 

.3149 

.1443 

1.91 

2  ACCOM 

M 

l 

t 

.4032 

.1029 

4.44 

3  ACCOM 

M 

1 

4 

.3114 

.1*0 

3.17 

4  UNEHP 

UP 

I 

i 

7*1.919* 

03.3072 

5.24 

5  UNEHP 

UP 

1 

7 

•493.0114 

24I.73M 

*2.04 

4  UNEHP 

UP 

1 

» 

442.9091 

240.2444 

1.93 

RESIDUAL  SUN  V  SONRES  ■  154*3750.5479*9 
DECREES  K  FREEOQN  41 

KSIIUAL  DEAN  SRHAK  •  Z57274.73042I 
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V 


ACF  VARUBLE  IS  RYI.  MltAG  IS  21./ 


HUMES  OF  OBSERVATIONS  *  47 
HEM  OF  THE  (DIFFERENCED)  SERIES  >  4.8272 
STMMRI  ERROR  OF  THE  HEM  >  59.5717 
T-VNJJE  OF  HEM  (RCRIHST  ZERO)  *  .Ntf 


AUTOCORRELATIONS 

1*  12  -.12  .13  .09  -.11  .18  -.11  -.IS  .14  -.12  -.14  .17  -.14 

ST.E.  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .12  .13 

13-  2«  .13  -.11  .11  -.14  Sli  -.11  -.18  -.14 

ST.E.  '  .13  .13  .13  .13  .13  .13  .13  .13 


I 


PUT  OF  SERIAL  CORRELATION  t 

-1.1  -.8  -.4  -.4  -.2  .1  .2  .4  .4  .8  1.1 


LAC 

CORR.  ♦ - ♦- — 

i 

1 

-.818 

♦ 

i 

2 

.13# 

* 

ii 

3 

.187 

♦ 

m 

4 

-.HA 

* 

i 

S 

.177 

♦ 

in 

4 

-.111 

♦ 

i 

7 

-.#44 

♦ 

ii 

1 

.134 

♦ 

ii 

9 

-.#19 

♦ 

i 

If 

-.#42 

♦ 

in 

11 

.173 

♦ 

in 

12 

-.144 

♦ 

ii 

13 

.132 

♦ 

u 

14 

-.1«S 

♦ 

IIII 

IS 

.#•9 

♦ 

i 

14 

-.154 

*  uni 

17 

-.#97 

♦ 

in 

-.#13 

* 

i 

19 

-.#79 

♦ 

in 

21 

-.#43 

♦ 

m 

« 
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pacf 


VARIABLE  IS  RU.  HAILAC  IS  21./ 


NUMBER  OF  OBSERVATIONS  >  47 
HEM  Of  THE  (DIFFERENCES)  SERIES  *  4.8272 
STANSARI  ERROR  OF  THE  MEAN  *  59.57*7 
T-VAUIE  OF  NEAR  (ACAINST  ZERO)  *  .Nil 


PARTIAL  AUTOCORRELATIONS 

1-  12  -.12  .03  .19  «.f  .17  -.12  -.IS  .02  -.01  -.IS  .07  -.13 

ST.E.  .12  *12  .12  .12  .12  .12  .!2  .12  .12  .12  .12  .12 

13-  21  .03  -.12  .02  -.18  -.08  -.01  -.14  -.04 

ST.E.  .12  .12  .12  .12  .12  .12  .12  .12 


PLOT  OF  SERIAL  CORRELATION 

-1.0  -.8  -.4  -.4  -.2  .0  .2  .4  .4  .0  1.1 

LAC  CQRR.  ♦ - ♦ - * - ♦- — — ♦- — ♦ 

I 


-.118 

♦ 

1 

.029 

♦ 

11 

.008 

♦ 

III 

-.004 

* 

I 

.072 

* 

III 

-.014 

♦ 

I 

-.050 

♦ 

11 

.022 

4 

II 

-.012 

♦ 

I 

-.043 

III 

.070 

4 

III 

-.029 

II 

.032 

4 

II 

-.114 

4 

IIII 

.024 

4 

II 

14 

-.178 

4  mu 

17 

-.082 

4 

III 

18 

-.011 

4 

I 

If 

-.042 

4 

II 

21 

-.058 

4 

II 

MUM 


*4 


VMIMLE  IS  IMENP. 

ARORBERS  ME  Mlt2l\ 

M VALUES  «  1.131,  -.1477. 
CENTEREI./ 


THE  COMMENT  IMS  KEN  MOEO  TO  THE  HOOEL 

THE  CURRENT  N08EL  HAS 
OUTPUT  VMIMLE  ■  UHENP 
INPUT  VMIMLE  *  HOIK 


FORECAST  CASES  ME  24.  JOIN./ 


FORECAST  M  VMIMLE  IMENP 

PER1M 

FORECASTS 

82 

9.4S7S* 

83 

9.4*947 

84 

9.34148 

85 

9.31434 

84 

9.24814 

87 

9.22294 

88 

9.17847 

89 

9.13531 

98 

9.19284 

91 

9.«SI3* 

92 

9. *1*4* 

93 

8.97*75 

94 

8.93173 

95 

8.89352 

94 

8.8S4II 

97 

8.81941 

98 

8.78341 

99 

8.74849 

IN 

8.71411 

1*1 

8.48*43 

1*2 

8.44747 

1*3 

8.41518 

1*4 

8.58358 

1*5 

1.55243 

FRM  TINE  PERIOO  82 

ST.  ERR.  ACTUAL 
.25451 
.387*4 
.48342 
.54191 
.42433 
.48321 
.73359 
.77885 
.81991 
.85742 
.89191 
.92374 
.95331 
.98179 
1.M444 
1.13*44 
I .*5293 
1 . *74*5 
1 .*9391 
1.11242 

1.13*27  , 

1.14494 

1.1427* 

1.1774* 


STAMARI  ERROR  «  .254499  (II  CONOITIOHAL  NETHOI  1 
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P$ HEIGHT  NAIPSI  «  «./ 


4#  PSI-MCIGMTS  ME  STOKED 


ERASE  MODEL./ 


UNIVARIATE  TIME  SERIES  MODEL  ERASED 


MINA  VARIABLE  IS  ACC13AH. 
ARORIERS  ARE 
ARVALUES  *  .S2S8i  .4415. 
NAMIERS  ARE  MU’. 
CENTERED./ 


THE  COMPONENT  HAS  KEN  ADDED  TO  THE  MODEL 

THE  CURRENT  MODEL  HAS  ' 

OUTPUT  VARIABLE  «  ACC13AH 
INPUT  VMIAILE  «  NOISE 


INBEP  VMIAILE  IS  UNEAP. 

UPORBERS  •  Ml»7ill'. 

UPVALUES  •  571. I.  -511.11  434.1. 
CENTERED./ 


THE  COMPONENT  HAS  KEN  ADDED  TO  THE  MODEL 

THE  CUMENT  NQKL  HAS 
OUTPUT  VMIAILE  «  ACC  13AM 
INPUT  VMIAILE  «  NOISE  UNCAP 
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f 

FORECAST  ON  VM1ASLE  ACC 1 3AM  FROM  TIME  PERIOD  82 


PERIOD 

FORECASTS 

ST.  ERR. 

02 

4411. 70357 

523.48040 

03 

4011.71812 

504.48025 

04 

4021.84704 

441.12405 

85 

4405.40440 

473.43440 

04 

4410.14441 

782.15512 

07 

4717.42421 

720.13057 

08 

4442.47102 

702.41033 

00 

44N.57373 

027.27401 

Of 

4450.40804  , 

057.83*48 

01 

4410.14000 

870.45010 

02 

4520.14543 

004.37084 

03 

4548.01222 

087.41050 

04 

4530.04447 

020.35123 

05 

4404.23172 

057.45814 

04 

4405.30230 

003.4*741 

07 

4470.73134 

1N4. 14821 

08 

4420.54704 

1810.35482 

00 

4410.20444 

1838.0501# 

IN 

4415.22137 

1844.27030 

lfl 

4377.82844 

1841.11580 

112 

4344.08N7 

1070.S45N 

IN 

4351.48223 

1804.73203 

lf4 

4322.48074 

1111.81501 

IN 

4383.41324 

1122.10305 

1#4 

4200.83427 

1132.08848 

MSEI  ON  THE  AOAILAILE  DATA*  ONU  25  FORECASTS  CAN  K  MADE 

STANDARD  ERROR  •  523.41#  III  CONDITIONAL  METHOD  ) 
l PACE  If  ACC13AH  WITH  UNEMP 
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Appendix  bI-7 

MDP  Output  for  Leading  Indicator  Transfer  Function 
Models  (with  Forecasts)  with  Producer  Price  Index 

as  Input  Variable 
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•PIT 


VARIABLE  ACCISAH  IS  FILTERED.  RESULTS  ARE  STORED  IN  VARIABLE  RV 
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EffRCft  «  S7B7S8  1ST  COtCITIOML  ICTMOO 


VARIABLE  ACC  UAH  IS  FILTERED.  RESULTS  ARE  STORED  IN  VARIABLE  RT 


ut  si  noiiM 
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irr  conditional  ictmoo 


ARINA 


VARIABLE  IS  AC LI  JAM. 
DFORDER  IS  1. 

ARORDERS  ARE  '(1,2)'. 
CENTERED./ 

FORECAST  CASE  IS  10.  START  IS  50./ 


FORECAST  ON  VARIABLE  ACC13AH  f  RON  TIME  PERIOD  50 

i 


PERIOD  FORECASTS 


50 

1748.97683 

51 

1678.35866 

52 

1655.71074 

53 

1685.05384 

54 

1684.34892 

55 

1674.95158 

56 

1677.77880 

57 

1680.06950 

58 

1678.51230 

59 

1678.19378 

ST.  ERR.  ACTUAL 
713.64942  1742.00000 
880.91520  1650.00000 
943.40842  24B8. 00000 
1045.89281  2504.00000 
1150.15525  2340.00000 
1229.25173  2203.00000 
1304.73549  2172.00000 
1380.47766  2398.00000 
1450.75478  2915.00000 
1516. 87421  2653.00000 


STANDARD  ERROR  =  713.649 


<BY  CONDITIONAL  METHOD  ) 


AR I  HA 


VARIABLE  IS  AC LI 3AM. 

DF ORDER  IS  1. 

ARORDERS  ARE  '(1,2)'. 
CENTERED./  t 

FORECAST  CASE  IS  10.  START  IS  60./ 


FORECAST  ON  VARIABLE  ACC13AH  FROM  TIME  PERIOD  60 


PERIOD  FORECASTS 


60 

2556.39389 

61 

2668.73581 

62 

2669.27230 

63 

2632.3999 9 

64 

2642.41359 

65 

2651.69991 

66 

2645.86041 

67 

2644.43839 

68 

2646.74024 

69  ' 

2646.56902 

STANDARD 

ERROR  -  661.074 

ST.  ERR.  ACTUAL 
661.07372  2375.00000 
816.01676  2495.00000 
873.90600  2323.00000 
968.84021  2500.00000 
1065.42146  3113.00000 
1138.69078  3029.00000 
1208.61353  2965.00000 
1278.77565  2649.00000 
1343.87534  2322.00000 
1405.12365  2998.00000 

(BY  CONDITIONAL  METHOD  ) 


ARINA  VARIABLE  IS  ACC13AN. 

DFORDER  IS  1. 

ARORDERS  ARE  '(1,2)'. 

FORECAST  CASE  IS  10.  START  IS  75./ 


FORECAST  ON  VARIABLE  ACC  1 3 Art  FROM  TINE  PERIOD  75 


PERIOD  FORECASTS 


75 

3194.09557 

76 

3155.86233 

77 

3178.81816 

78 

3184.97302 

79 

3175.76808 

80 

3176.29970 

81 

3179.16185 

82 

3178.19716 

83 

3177.52811 

84 

3178.02834 

ST.  ERR.  ACTUAL 

602.63919  3630.00000 

743.88630  4455.00000 

796.65851  4454.00000 

883.20116  4933.00000 

971.24528  4396.00000 

’1038. 03807  4455.00000 

1101.78012  5097.00000 

1165.74037 
1 225 .013568 
1280.92005 


STANDARD  ERROR  -  602.639 


(BY  CONDITIONAL  METHOD  ) 


FORECAST  CASE  IS  10.  START  IS  50./ 


FORECAST  ON  VARIABLE  ACC13AA  FROfl  TIME  PERIOD  50 


PERIOD 

FORECASTS 

ST.  ERR. 

ACTUAL 

50 

1631.13598 

608.05057 

1742.00000 

51 

1624.99382 

784.98051 

1650.00000 

52 

1646.37702 

867.48427 

2488.00000 

53 

1662.90076 

963.67506 

2504.00000 

54 

1854.64815 

1059.92722 

2340.00000 

55 

1701.97338 

1141 .99861 

2203.00000 

56 

1521.88813 

1217.76167  , 

2172.00000 

5? 

1653.96637 

1290.51636  ‘ 

2398^00000 

58 

1860.62967 

1359.30526 

2915.00000 

59 

1766.05508 

1424.47338 

2653.00000 

STANDARD  ERROR  =  608.051  <BT  CONDITIONAL  METHOD  ) 

1  PAGE  11  ACCESSIONS 


FORECAST  CASE  IS  10.  START  IS  60./ 


FORECAST  ON  VARIABLE  ACC13AH  FROM  TINE  PERIOD  60 


PERIOD 

FORECASTS 

ST.  ERR. 

ACTUAL 

60 

2549.33471 

549.21368 

2375.00000 

61 

2632.00274 

709.02331 

2495.00000 

62 

2512.52297 

783.54375 

2323.00000 

63 

2542.93764 

870.42681 

2500.00000 

'  64 

2598.95116 

-  957.36531 

3113.00000 

65 

2556.72016 

1031.49521 

3029.00000 

66 

2774.60056 

1099.92720 

2965.00000 

67 

2749.14360 

1165.64191 

2649.00000 

68 

2657 . 1*0  415 

1227.77459 

2322.00008 

69 

2643.60233 

1286.63683 

2998.00000 

352 


ARINA 

VARIABLE  IS 

ACC13AM. 

DFORDER  IS  1 

• 

ARORDERS  ARE 
CENTERED./ 

'(1,2,14)'. 

FORECAST 

CASE  IS  10. 

START  IS  75./ 

FORECAST  ON 

VARIABLE  ACC  13 AH  FRON  TIHE 

PERIOD  75 

PERIOD 

FORECASTS 

ST.  ERR. 

ACTUAL 

75 

3239.12204 

496.13972 

3630.00000 

76 

3155.69930 

640.50595 

4455.00000 

77 

3216.53659 

707.82501 

4454.00000 

78 

3389.29928 

786.31202 

4933.00000 

79 

3320.38748 

864.84911 

4396.00000 

80 

3273.98160 

931.81537 

4455.00000 

8! 

3244.86950 

VV3. 63435 

50V7.00000 

82 

3149.68556 

1052.99864 

83 

3354.49160 

1109.12704  ■ 

84 

3470.80177 

1162.30106 

STANDARD  ERROR  =  496.140 

1 P AGE  9  ACCESSIONS 


(BY  CONDITIONAL  METHOD  ) 


ARINA 


VARIABLE  IS  ACC13AM. 
DFQRDER  IS  t. 

ARORDERS  ARE  '(1,2,14)'. 
CENTERED./ 


FORECAST  CASE  IS  10.  START  IS  80./ 


FORECAST  ON  VARIABLE  ACC  1 3 AN  FRQN  TIME  PERIOD  BO 


PERIOD  FORECASTS 


80 

4361.11308 

81 

4413.59049 

82 

4325.12868 

83 

4509.10059 

84 

4633.25590 

85 

4461.12185 

86 

4456.14689 

87 

4496.20323 

88 

4529.84149 

89 

4620.19912 

ST.  ERR.  ACTUAL 

502.22676  4455.00000 

648.36419  5097.00000 

716.50918 
795.95913 
875.45978 
943.24763 
1005.82505 
1065.91768 
1122.73471  ■ 

1176.56110 


STANDARD  ERROR  -  502.227  (BY  CONDITIONAL  METHOD  ) 

1PA6E  9  ACCESSIONS 
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VITA 

Kenneth  Michael  Kalinich  was  bom  on  30  November  1951  in  Endi- 

cott,  New  York.  He  enlisted  in  the  Army  in  July  1969  to  attend  the 

United  States  Military  Academy  Preparatory  School  at  Ft.  Bel voir, 

Virginia.  He  attended  and  graduated  from  the  United  States  Military 
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Permanent  address i  2473  Felters  Rd. 

Binghamton,  New  York  13903 
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VITA 


Dennis  Charles  Wenzel  was  born  on  1  August  1949  in  San 
Antonio,  Texas.  He  enlisted  in  the  Army  in  February  1972 
and  served  in  Korea  as  a  computer  programmer  from  September 
1972  to  October  1974.  Under  an  Army  enlisted  student  status, 
he  received  a  Bachelor  of  Arts  in  Mathematics  and  a  Bachelor 
of  Business  Administration  in  Personnel  Management  from 
Saint  Mary's  University  of  San  Antonio  in  May  1976.  He  was 
commissioned  in  October  1976  and  served  in  Korea  for  one 
year  as  a  computer  data  base  manager.  In  January  1978  he 
performed  duty  as  a  personnel  management  officer  in  th|S  6th 
A^r  Cavalry  Brigade  and  as  a  systems  analyst  in  III  Corps 
at  Fort  Hood,  Texas.  In  1980  he  was  assigned  six  months 
Temporary  Duty  to  Fort  Harrison,  Indiana,  as  a  Data  Processing 
Officer  Student  and  Instructor.  He  then  attended  the 
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Engineering,  Air  Force  Institute  of  Technology,  in  May  1981. 

Permanent  address:  4236  Chestnut  Hill 

San  Antonio,  Texas  78218 
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